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Abstract 



We employ the theoretical framework of positive operator valued measures, to study 
Markovian open quantum systems. In particular, we discuss how a quantum system 
influences its environment, or parts thereof. Using the theory of indirect measurements, 
we then draw conclusions about the information we could hypothetically obtain about 
the system by observing the environment. Although the environment is not actually 
observed, we can use these results to describe the change of the quantum system due to 
its interaction with the environment. We apply this technique to two different problems. 

In the first part, we study the coherently driven dynamics of a particle on a rail of 
quantum dots. This particle can tunnel between adjacent quantum dots and the tunnel 
rates can be controlled externally. We employ an adiabatic scheme similar to stimulated 
Raman adiabatic passage in quantum optics, to transfer the particle between different 
quantum dots. 

We compare two fundamentally different sources of decoherence. The first is the 
frequent but weak measurements of the position of the particle by nearby quantum 
point contacts. This Markovian effect destroys the position coherence required for a 
good particle transport fidelity. Second, we couple the particle to two-level fluctuators, 
which are randomly located near the quantum dot rail. These result in non-Markovian 
dephasing, which reduces the transport fidelity in a quite different way. 

In the second and larger part of this thesis, we study the dynamics of a free quantum 
particle, which experiences random collisions with gas particles. Previous studies on this 
topic, which applied scattering theory to momentum cigenstatcs, found controversial 
results and were not conclusive. Therefore, we present a supplementary approach, where 
we start by solving the time dependent Schrodinger equation for two colliding particles, 
each described by a Gaussian wave function. 

Next, wc use our results about a single collision, to develop a rigorous measurement 
interpretation of the collision process, in which the colliding gas particle performs a 
simultaneous momentum and position measurement on the tracer particle. This in turn 
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leads to the collisional transformation of the tracer particle's density operator. We then 
derive a master equation by using the correct collision statistics. Finally, we study the 
collisional decoherence process in terms of the Wigner function. 

Because our approach is more complex than the typical scattering theory approach, 
we restrict ourselves to one spatial dimension. Nevertheless, we find some interesting 
new insight, including that the previously celebrated quantum contribution to position 
diffusion is not real, but a consequence of the Markovian approximation. Further, we 
discover that the leading decoherence process is due to phase averaging, rather than 
induced by the information transfer between the Brownian particle and the gas. 
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Acronyms and Notations 



Throughout the thesis, we used the following acronyms: 



Table 1: Acronyms 



CTAP 


coherent tunneling by adiabatic passage 


POVM 


positive operator valued measure 


QBM 


quantum Brownian motion 


QD 


quantum dot 


QLBE 


quantum linear Boltzmann equation 


QPC 


quantum point contact 


QPD 


quantum contribution to position diffusion 


STIRAP 


stimulated Raman adiabatic passage 


TLS 


two level system 



We used hats for all operators acting in Hilbert space, except for the density operator, 
the free evolution operator, and the Hamiltonian. We further used calligraphic fonts for 
super operators and tildes for continuous measurement outcomes. The index g refers 
to a gas particle. Below is a list of the notations, which where used in more than one 
chapter. 



XV 



XVI 



ACRONYMS AND NOTATIONS 



Table 2: Notations 



a 


distance between QD rail and 




tunneling rate times h 




QPC rail 




p 


momentum 


a 


measurement sensitivity (part 


p 


momentum after collision 




I) 






p 


momentum operator 


a 


mass ratio mg/m (part II) 


p 


momentum measurement re- 


A 


measurement Kraus operator 




sult 


Xn 


coupling 


constant to n-th TLS 




effect operator 


d 


distance between neighboring 


m/ m 


J. J. J. £ ■ 1 / 

state vector of a smgle / com- 




QDs 








posite system 


d 


differential 




q 


momentum transfer 


D 


decoherence rate 




distance between i-th QD and 




momentum diffusion constant 






Dxx 


position diffusion constant 


p{x,p) 


classical phase space probabil- 


D{x,p) 


Glauber- 


Sudarshan displace- 




ity distribution 




ment operator 


P / Q 


density operator of a single / 


r 




coarse graining time 




composite sysiem 


E 


energy 






K 


measurement rate 


7 


friction constan 


K, K[Pg), 


collision rate operators 


phase space projection opera- 


RsiXgjPg) 




tor 








total scattering cross section 


h 


reduced 


Planck constant 


a 


Tt{x = 0,p = 0) 




/i/(2vr) 








Pauli z-operator for Ti-th TLS 


H 


Hamiltonian 


s 


scattering operator 


n 


Hilbert space 


tc 


collision time 




imaginary unit 


T 


temperature 




Boltzmann constant 


f 


transition operator 


K 


see Eq. ( 


5.7 


) 


U{t) 


free evolution operator 


c 


Liouville super operator 


V 


velocity 


m 


mass 






V 


volume 




Maxwell-Boltzmann distribu- 


w 


position variance of Gaussian 




tion 








wave packet 


Hg 


gas particle density 




Wigner function for the state 


6 


general operator 


P 





general observable 


ry* nf* rv* /y* 


analogous to p's 



Chapter 1 

General Introduction 



All physical systems, whether governed by the laws of classical or quantum physics, are 
open systems, that is, they interact to some extent with their surrounding environment. 
For quantum systems, this interaction leads to the process known as decoherence, that 
has the effect of destroying the peculiar quantum feature of such systems, which is that 
quantum systems can behave as if they are simultaneously in a number of different, 
distinct states. This is both, a benefit and a curse: decoherence is believed to be 
responsible for the 'emergence' of the classical laws of physics from their underlying 
quantTim form, hnt it also makes it difficult to exploit these same quantum features in, 
for instance, the development of quantum computers. 

One of the most fundamental differences between classical and quantum physics are 
their respective descriptions and interpretations of a measurement process. The classical 
one is rather trivial, in that a measurement changes our knowledge about the state of 
the measured system, but does not necessarily affect the state of the system itself. In 
fact, a measurement just means using our senses to gain information about a system, as 
e.g. our eyes for its position. If we can not use our senses directly, cither because they 
lack in precision, or because we do not have the required sense for a physical observable 
(i.e. magnetic field), we employ a classical measurement apparatus. This magnifies the 
desired observable and / or translates it into another observable for which we have a 
sense, as the position of a pointer, or the click of a Geiger-Miiller counter. We say, the 
position of the pointer becomes correlated to the observable. 

A quantum mechanical measurement process also translates and magnifies a physical 
observable, but in general, it also irreversibly changes the state of the measured system. 
The reason is that in quantum mechanics, the building up of correlations between two 
systems (observed system and measurement apparatus), necessarily changes the states 



1 



2 



CHAPTER 1. GENERAL INTRODUCTION 



of both. 

If we are interested in the evohition of the state of a quantum system, we should 
therefore know when a measurement on it is performed. Of course, from the systems 
point of view, a measurement is performed each time it becomes correlated with any 
other system, regardless whether or not the latter can be observed by one of the human 
senses. We conclude that a quantum system is also measured by its environment [1]. 

This poses the following questions: Can we use the very well developed theory of 
quantum measurements, to describe a general open quantum system? If so, is this new 
approach advantageous to the many other means of studying open quantum systems? 
Rather than trying to answer these very general questions, in this thesis, we will apply 
the above reasoning to two special problems, which are both of general interest to the 
physics community. In both instances, we find that our measurement approach presents 
us with the necessary insight for a good interpretation of our results. 

A challenge in our approach is certainly to work out the details of the measurements, 
which are carried out by the environment. But once that is achieved, in some instances, 
it is straightforward to write down a master equation, which governs the dynamics of 
the quantum system of interest. A further nice property of our approach is, that in 
addition to the mathematical master equation, it naturally presents a physical interpre- 
tation thereof. Especially in the second part of our work, this leads us to a number of 
conclusions, which were not found in previous studies. 

We wish to note that the idea of a measurement interpretation of the influence of 
an environment is not new. Let us quote Joos and Zeh in their influential paper about 
collisional decoherence [2j: The destruction of interference terms is often considered 
as caused in a classical way by an 'uncontrollable influence ' of the environment on the 
system of interest. In fact, this interpretation seems to date back to Heisenberg J^. 
But the opposite is true: the system disturbs the environment, thereby dislocalizing the 
phases. If the system is initially in a superposition Yin '^n \i^n) , it may influence the 
environment as if being measured by it. 

As this point of view is now widespread, especially in regard to collisional decoherence 
and quantum Brownian motion, it is surprising that there is no rigorous attempt in the 
literature, to work out the precise details of such an environmental measurement. In 
the second part of this thesis, we present these details for the very topic of collisional 
decoherence. Interestingly, these measurements turn out to be too "weak" to account 
for a leading contribution to decoherence. Instead, we find that the main cause of 
decoherence is indeed the 'uncontrolled influence' of the random gas particle momentum. 



Chapter 2 

General Concepts in Quantum 
Mechanics 



In this chapter we introduce the basic concepts used throughout this thesis as well as 
establish the notation. There are many sources for the material presented here: the 
lecture notes of Hornberger 0] provide a good introduction, with a more thorough 
development to be found in Breuer and Petruccione [S]. 

2.1 Density operators and composite quantum systems 

One of the basic tenets of quantum mechanics is that a closed quantum mechanical 
system S is described by a state vector, which is a normalized element of a Hilbert space 
Us- We will use Dirac's bra-ket notation when talking about state vectors \^)s, where 
the subscript refers to the quantum mechanical system of interest. The generally time 
dependent state vector describing the system S evolves according to the Schrodinger 
equation 

^^\m)s = ^Hs{t)\m)s, (2.1) 

where Hs is the Hamiltonian of the system S. 

An alternative approach to quantum mechanics is the use of a state operator ps{t) 
to describe S. This operator is known as the density operator, or density matrix when 
talking about its components with respect to some basis of Us- In the most simple 
case, when a ket \ip{t))s can be assigned to the system, then the corresponding density 
operator is the projector onto this ket, i.e. ps(t) = \tp{t))s{tp(t)\. It follows from 



Eq. (2.1 ), that the density operator then evolves according to the von Neumann equation 



^^ps{t) = ^[Hs{t),ps{t)], (2.2) 



4 



CHAPTER 2. GENERAL CONCEPTS IN QUANTUM MECHANICS 



where [•, •] denotes the commutator. 

One of the major advantages of the density operator formahsm over the state vector 
formahsm is that it includes the description of statistical mixtures. Such mixtures are 
commonly needed if one is not certain about the state of the quantum system, but only 
knows the probability pj for the system being in the state \(pj)- Then, one can use the 
density operator 

PS = ^Pj\<Pj)s{(pj\- (2.3) 
j 



which will again evolve according to the von Neumann equation (2.2), if the system S 
is a closed one. 

Some basic properties of the density operator follow directly from the probability 



interpretation of pj in Eq. (2.3), in particular its positivity, i.e. s {(Pl Ps \ (t>)s ^ for all 
\(p)s £ 'Hs, and its normalization Trps = 1. If the state of the system is known with 
certainty, that is, if the density operator is a projector, then the system is said to be in 
a pure state, whereas otherwise we say the system is in a mixed state. 

It is important to note that for a given density operator ps, the decomposition 



Eq. (2.3) is not unique. For example, a system being in the state \(t)i)s or \4>2)s with 
respective probabilities 1 /2 has the same density operator as a system being in the state 
{\4>i)s + \(t'2)s)/V^ or {\4>i)s — \4'2)s)/V^ with the same probabilities. This imposes 
the question whether the density operator formalism is sufficient to describe all physical 
situations. The answer is a definite yes if one is interested in the prediction of probabil- 
ities of measurement outcomes only, as is usually the case in any physical theory. More 
precisely, it can be shown that no measurement can distinguish between two physical 
situations which are described by the same density operator. 

The most general density operator can always be written in terms of an orthonormal 
basis \i)s as 

PS = '^Cii'\i)s{i'\, Cii>0, ^Cn = l, (2.4) 

where cn' = si^l Ps Ws- Ps is positive and therefore Hermitian, one can choose the 
eigenvectors of ps as orthonormal basis to diagonalize the density matrix 

PS = ^Pi\i)s{i\, Pi>0, ^pi = l, (2.5) 

i i 

where pi is the eigenvalue of ps corresponding to the eigenvector \i)s- 
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The concept of density operators is most appreciated when the system S of interest 
is part of a larger, composite quantum systems S + T. If \i)s and \j)T are respective 
orthonormal bases of Hs and Ht, then \i)s <8) |i)r is a basis for the product Hilbert 
space H = 'Hs'Si'Ht of the combined system. The symbol (g) denotes the direct product, 
and will frequently be omitted for shorter notation. Any pure state of the combined 
system can now be written either as state vector \^)s+t = Ylij lij \^)s b)r or as density 
matrix qs+t = Ylu'jj' '^u'jj' \j)T{j'\, and the two descriptions are related by 

Cii'jj' = lijlitji- If the state vector of the composite system can be written as a direct 
product |*)5+T = \4')s ® W)t, then the composite system is said to be in a product 
state. If this is not the case, then the subsystems S and T are called entangled and this 
is the general situation if the subsystems are interacting with each other (or did interact 
in the past). 

If there is entanglement between S and T, then it is not possible to assign a state 
vector to any of the subsystems. On the other hand, the density operator formalism 
tells us how the reduced density operators for each subsystem is obtained in terms of the 
partial trace of the density operator of the composite system: 



with Cjj/ = '^ii'jj- Note that if the two systems are entangled, the reduced density 
operator of each subsystem is never in a pure state even if the composite system is. 

It is this property of the possibility to describe subsystems of larger systems, which 
makes the density operator formalism so suitable for the study of open quantum sys- 
tems, where the system of interest is only a small part of a miich larger system often 
incorporating an environment with infinite degrees of freedom. 

An important property of density matrices is Trp|. < 1 where the equal sign holds 
if and only if ps describes a pure state. Therefore Trp|. is often used as a measure of 
how close the system is to a pure state and is commonly referred to as the purity. If 
S* is a subsystem of S" + T and the combined state is pure, then Trp|. can also serve as 
a measure of entanglement between S and T. Whether a composite system described 
by a mixed state qs+t is entangled or not is not as easily seen. One can show that a 
system S" + T is not entangled if and only if its density operator is separable, that is if it 
can be written as Qs+T = J2m lmPs,m ® PT,m where ps,m and pr,m are density matrices 
of the respective subsystem and > 0. 

Unitary operations 1^)^ — >■ \'4^')s = Us\'4^)s which are well known from the state 



PS 




(2.6) 
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vector formalism of quantum mechanics are easily generalized to density operator trans- 
formations by ps — )• p'g = UspsUg. But there is a larger class of "allowed" linear 
transformations of density operators. Any linear map which maps operators on oper- 
ators is called a super operator and will be denoted with calligraphic fonts. A super 
operator : Ps ^ p's — Tsips) is called positive if Ts{ps) is positive for all positive 
Ps- It is called completely positive if (8> 1 is positive for all identity operators 1 in 
any dimension. An "allowed" linear map Ts is trace preserving and completely positive. 
Trace preserving and positivity are required for the transformed density operator to 
be a valid density operator. Complete positivity is less intuitive at first glance, but is 
required because there is always another system T which does not interact with S (e.g. 
a molecule on the moon). The combined system is then transformed by (8) It which 
again should be a positive operation therefore demanding complete positivity for Ts- 
Any such transformation can be written in terms of Kraus operators 

Tsips) = ^KmPsKl (2.7) 

m 

with KmKli = 

Remember that any linear operator O is uniquely specified by its action on a complete 
basis i.e. knowing O \i) for all i is equivalent to knowing O- Similarly, Ts{\i){i'\) for 
all i,i' uniquely determines the super operator Ts- Note that the action of Ts on all 
diagonals \i){i\ alone is not sufficient to uniquely determine Ts- 

2.2 Decoherence 

A diagonal density operator ps = YliPi is sometimes called "classical", because 

it can be interpreted as having a state \i)s with the classical probability pi. If there are 
"quantum" superpositions between states \i)s and \i' ^ i)s involved, then the density 
operator will also include off-diagonals, which are termed coherences- The coherences 
therefore reflect quantum behavior. 

The reduction of coherences is called decoherenc^ and leads to a classical like be- 
havior of a quantum system. This process typically results from the coupling to an- 
other system, which can be either classical or quantum in nature. If studying quantum 
properties of a system, one aims to reduce decoherence effects as much as experimen- 
tally feasible. For this reason, decoherence is a major obstacle in quantum information 
science, where superpositions of many "quantum bits" have to be maintained over a 
^The nomenclature to this topic varies vastly in the literature. A good overview is given in [3]. 
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sufficiently long period of time. On the other hand, decoherence is necessary for any 
measurement process in quantum mechanics, and is believed to be responsible for a 
classical macroscopic world. 

Of course, whether a density operator is diagonal or not, depends on the basis 
\i)s. Therefore, when studying decoherence, one has to single out a basis of the Hilbert 
space. Often, the microscopic decoherence mechanism points out an appropriate bases in 
which decoherence is most pronounced. In this pointer basis, all coherences often vanish 
completely over time. One should mention that many problems regarding pointer bases, 
are not resolved yet py, and the concept remains somewhat vague, except for some 
special cases. 

In the most simple situation, when the diagonals are constants in time, the off- 
diagonals often follow exponential decay 

s{i\psit)\iys = e-^-'*5(i|p5(0)|i>'5. (2.8) 

The decoherence rate Dn' is often called dephasing rate in processes which do only affect 
the off-diagonals. If dephasing is due to multiplication of coherences by a random and 
unknown phase, then we will also call it phase averaging. 

A somewhat different definition of decoherence is the loss of purity. This notion 
is independent of any designated basis, but instead depends on the initial state of the 
system. 

2.3 The von Neumann-Liiders measurement 

The von Neumann-Liiders measurement [HI E] is easily generalized to density operators. 
As in the state vector formalism, one assigns an Hermitian operator 

d = Y,rm\m)s{m\ (2.9) 

m 

to the observable to be measured. Here, are the possible measurement outcomes 
which will be found with the probability 

P{rm) = s{m\ps\m)s. (2.10) 

The expectation value of an observable measured in a von Neumann-Liiders measure- 
ment is easily obtained via the trace operation 

(O) = ^P(r„)r™ 

m 

= Tv{dps). (2.11) 
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After a measurement result is found, the system is known to be in the measured 
state \m)s{m\. That is, if a second measurement of the same observable is performed, 
the measurement outcome will be with certainty. This change of the system's state 
into the measured one is known as the projection postulate of quantum mechanics. The 
set of states {l^n-)} are called the measurement basis of the measurement. 

2.4 Generalized measurements 

For both parts of this work the concept of generalized measurements IH O |8] will be 
of importance. We will refer to this as the theory of positive operator valued measures 
(POVM). It generalizes the von Neumann-Liiders measurement to a larger class of mea- 
surements. The von Neumann-Liiders measurements can be considered as the limit of 
a strong measurement in that it ensures that the measured system is in the measured 
state after the measurement. POVMs also include measurements where the measured 
system will generally not end up in the measured state, but might more or less differ 
from this state. Consecutive measurements on the same system will then generally give 
different results, but there will be some correlation between probabilities of consecutive 
measurement results (if there were no correlation at all, the name "measurement" would 
barely be justified). If the correlation between consecutive measurements are small, the 
measurement is called a weak one. These typically disturb the system only slightly, but 
do not reveal much information about the measured system. 

If a measurement is performed on a system described by a general density operator 
Ps, then for every possible measurement outcome there exists a positive operator, 
the so-called effect operator TTm such that the probability of this measurement outcome 
is 

P(r„) = Tr(psvrm). (2.12) 

The normalization of probabilities requires Ylm^"^ ~ ^s- The state of the system 
changes due to the measurement with outcome according to 

= Tr(p5-^) ^'-''^ 
where J2n Amn-Amn = TTm. It follows that the density operator after a non-readout mea- 
surement, that is after a measurement is performed but without knowing the outcome, 
is 

P'S = J^^mnPsiL- (2-14) 
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The operators Amn will be called Kraus operators according to their more general use 
outside of measurement theory (see previous section). A special case which is often 
encountered is if for every effect operator there is only one Kraus operator — 
Am Am- In this case, the measurement is called efficient and the transformation ps — >■ 
AmPsAm/TT{AmPsAm) maps pure states to pure states. In some sense, these measure- 
ments can be considered as optimal, because no unnecessary decoherence is introduced 
during the measurement process. 

The generalized measurements include the von Neumann-Liiders measurements as 
the special case TTm = Am = Pm where Pm are projection operators. Interestingly, one 
can also derive POVMs from indirect von Neumann-Liiders measurements. For this, 
one couples the system S to some ancilla system M which is often called measurement 
apparatus or meter in the literature (although it can only be considered as the "first step" 
of a measurement apparatus). After some coupling time, which is usually assumed to 
be short compared to all other time scales, the ancilla's state depends on the initial 
state of S. Performing a von Neumann-Liiders measurement on the ancilla will then 
reveal information about the state of S. The effect and Kraus operators acting on S to 
describe this indirect measurement may depend on the initial state of M, the coupling 
Hamiltonian Hj^is{t)^ and on the type of measurement performed on M. 

Because it is possible to show that all POVMs can in principle be derived in such 
a way [S], one is tempted to consider the von Neumann-Liiders measurement as the 
fundamental one. This would however be misleading as often strong von Neumann- 
Liiders measurements are actually the result of a large number of weak POVMs. 

One should mention here that there is still no satisfactory microscopic derivation of 
a quantum mechanical measurement process. Especially the projection postulate or its 



generalization Eq. (2.13) has to be used without further justification [8]. 

There are two distinct causes for the possibility of an indirect measurement to reveal 
only little information about the system. First the coupling of M to the system S might 
be a weak one in a sense that the combined state of M and S is only little influenced by 
the coupling. Then M and S are only weakly correlated and even the best measurement 
on M gives only little information about S. Such a measurement influences the system 
only little, and is therefore called a weak measurement. In this case the effect operator 
is far away from being a projection operator, as are the Kraus operators. Such a weak 
measurement can also be an efficient one. 

If the coupling between M and S is optimal in that M and S are perfectly cor- 
related in the desired measurement basis of S, it is still possible that a measurement 



10 



CHAPTER 2. GENERAL CONCEPTS IN QUANTUM MECHANICS 



on M gives only little information about S. This can happen either because the von 
Neumann-Liiders measurement on M has a different measurement basis which is not 
the basis of the correlations, or because the measurement on M is a weak POVM itself. 
In both cases, the effect operators acting on S indicate little information gain, whereas 
the Kraus operators indicate strong disturbance on S and can even be projectors. Such 
a measurement will influence the system considerably despite giving only limited infor- 
mation about it. These sort of measurements will never be efficient ones. In some sense, 
on a microscopic description a strong measurement is performed, but then followed by 
an imprecise readout of the meter. 



2.5 Open quantum systems 

If the quantum system S of interest is coupled to an environment E, then we call it 
an open quantum system. The environment is usually a quantum system with a large 
number of degrees of freedom, such as the electromagnetic field in a vacuum or the 
phonons in a solid state matrix. The combined system is then a closed system (see 



Fig. 2.1) evolving according to the von Neumann equation 



^PS+E{t) = \[Hs+Eit),QS+Eit)]. (2.15) 

at ih 

Once this equation is solved one could trace out the environment to find the density 
operator ps{t)- However, because of the many degrees of freedom of the environment, 
the solution of the von Neumann equation is generally not feasible. Instead one usually 
tries to find a first order differential equation for the reduced density operator ps{t) = 
TieQS+e^ which is called the master equation. In open systems (classical as well as 
quantum), the change of the state during an infinitesimal time interval (t,t + dt) does 
not only depend on the system's state at time i, but also on its past. That is because the 
state of the environment at time t depends on the history of the state of the interacting 
system 5*, and the evolution of S in turn depends on the state of E. This process is 
called hack action, and the open system is then said to behave non-Markovian, or the 
environment is said to exhibit memory. The master equation can be shown, by, for 
instance, the Zwanzig-Nakajima projection operator technique, to take the following 
form 



^^ps{t) = ^ dt'}Cs{t-t')[ps{t% (2.16) 
where the super operator ICs{t — t') is the memory kernel. 
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Figure 2.1: The system S interacts with the environment. The environment's state then changes 
depending on the system's state. Therefore, at some later time the system will see an environment 
which depends on its own past. This process is called back action, and leads to memory effects. 

The Markovian approximation ICs{t — t')[ps{t')] = 5{t — t')Cs[psit')] results in a 
master equation which is local in time. It is generally valid on a coarse grained time 
scale on which correlations between environment and system vanish. It can be shown 
that under quite general conditions the Liouvillian Cs has to take the following Lindblad 
form 191 



in order to ensure completely positive and trace preserving dynamics. The Hamiltonian 
consists of the system Hamiltonian Hs plus a Lamb-shift term H^s which is a coherent 
modification due to the coupling to the environment. The non-coherent effects of the 
environment are described by the Lindblad operators L^- 

It should be noted here that it is often very difficult to derive a Markovian master 
equation from the underlying microscopic physics. But once this is achieved it provides a 
clear physical picture in terms of quantum trajectories j5J, and there are many methods of 
solving such equations. One should however not forget, that any Markovian description 
is an approximate one. It involves one or several of the following assumptions [5]: 

1. weak coupling to the environment, 

2. high temperature of the environment, 

3. low density of the environment, 

4. instantaneous change of state at random times (e.g. collisions, measurements, ...). 




m 
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Furthermore, any Markovian master equation is only valid on a coarse grained time scale, 
not resolving times small compared to correlation time between system and environment. 

An example which will be used in this thesis is if the environment can be described 
as a measurement apparatus, which performs measurements at random times at a rate 
R [To]. In this case the Lindblad operators are the Kraus operators from the theory of 
POVM, multiplied with \/^- It is then clear that the time scale on which the master 
equation is valid is the time a single measurement takes. 

Most work about open quantum systems in the literature is concerned with Marko- 
vian systems described by a master equation of Lindblad form. This is on the one hand 
because in many physical situations the aforementioned assumptions are good approx- 
imations, and on the other hand because general non-Markovian master equations are 
much harder to solve. We will encounter non-Markovian behavior in the first part of 
this work, but we limited ourselves to a very small "environment" which enables us to 
solve the full dynamics of S + E. 



Part I 



Coherently Driven Adiabatic 

Transport 
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Chapter 3 

Introduction 



In the last two decades, the field of quantum computing has drawn the attention not 
only of a large number of physicists, but also mathematicians, computer scientists, 
engineers and to a smaller extent even the general public. At first glance this might 
be surprising, as we are not anywhere near of being able to build a useful machine for 
quantum information processing, but in fact there are many reasons for the excitement 
about the topic of quantum computing. 

Quantum mechanics these days is very well understood and was used for the devel- 
opment of a broad range of items which influence our lives. Just to name a few, there 
are lasers and all its applications; improvements in semiconductors which are used in 
LEDs, solar cells or classical computers; nuclear magnetic resonance essential for mag- 
netic imaging in medicine; nuclear fission as an energy resource as well as for radiation 
needed for medical purposes; and many more. Yet it is very hard to imagine what 
consequences a quantum computer would have in our daily life. The situation might 
be compared with the introduction of electronics during the last, say six decades. Elec- 
tromagnetism was well understood at that time and its applications already changed 
the world. But yet a new revolution was to come with electronics, the processing of 
information by means of electromagnetism. 

Quantum physicists are very excited about the idea of quantum computing because 
it is believed that only such a machine could be capable of simulating many particle 
quantum systems. Mathematicians are fascinated by the challenge a quantum computer 
imposes on the Church- Turing thesis [llj, which was long believed to be true without 
ever being proven. For computer scientists a whole new range of possible algorithms, 
so-called quantum algorithms open up. Finally, engineers are closest to the technical 
challenge of building such a machine. 
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One of the main challenges to a quantum computer is decoherence, a process which 
leads to classical like behavior of a quantum system. It is either caused by imperfect 
control of physical parameters, or by coupling to an environment. As an environment 
can be considered to change physical parameters in an non-controllable way, these two 
causes of decoherence behave very similarly and can generally be described within the 
same language. A problem lies in the fact, that systems which naturally couple only 
very weakly to any environment (e.g. nuclear spins), also couple weakly to measurement 
apparatuses, creating a new challenge in the readout of a computational result. 

In this part of our work, we are concerned with transfer of quantum information as 
it would be necessary to transport information in any computational device. Especially, 
if a two-qubit operation between distant qubits is needed (as will be in any realistic 
computation), the two qubits have to get close together to interact. One could also 
imagine that during the processing, quantum information has to move from some storage 
area to a processing area. 

In particular, we are concerned with information stored on a single quantum system. 
To be specific, we will use an electron and discuss the encoding of information into its 
spin (spin qubit), as well as into its location (charge qubit), and it will turn out that 
the encoding into the position is somewhat disadvantaged over the spin qubit. But our 
work equivalently applies to ions, where the information can be encoded in its energy 
levels, or other systems. This electron is allowed to tunnel between discrete positions, 
which we assume to be quantum dots (QD) fixed in space, but again one could use ion 
traps or similar. 




Figure 3.1: The electron carrying the quantum information can be located on any of the 
quantum dots (blue circles) of the quantum dot rail. It can move from one dot to another by 
tunneling, which can be controlled externally by applying voltages to the gates between them. 



We suppose that the tunneling rates between the QDs can be controlled externally, 
in our specific example by voltages applied to external gates |12j as is shown in Fig 3.1 
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We will use a particular transport concept called Coherent Tunneling by Adiabatic 
Passage (CTAP) [I2j, which in a different application is well know in quantum optics as 
Stimulated Raman Adiabatic Passage (STIRAP) [13]. Indeed, it was recently shown \14:\ 
[T5] by a Schrodinger-wave description with experimentally feasible parameters, that 
CTAP can be applied to a system where the QDs are replaced by phosphorus ions 
embedded in a silicon matrix. This latter system is a promising candidate for the 
implementation of quantum computing [16] . 

As an adiabatic process, CTAP does not require a very precise control over the 
tunneling rates. All what is required is the possibility to suppress tunneling and to 
switch it on slowly, but the precise time when tunneling is switched on as well as the 
tunneling rate need only qualitative control. This advantage however comes with a 
downside in that adiabatic processes are typically slow, therefore giving the environment 
more time to destroy quantum behavior. 

We include two fundamentally different types of environmental decoherence in our 
study, and investigate their influence on the transfer fidelity. First, we assume that 
there are measurement apparatuses near the chain of QDs to measure the position of 
an electron. This could be a realistic situation as we might need these apparatuses for 
readout of quantum information, but more generally there are many types of Markovian 
environments (including dephasing), which have very much the same effect as such 
a measurement apparatus. Therefore our model for Markovian decoherence is quite 
realistic in many practical circumstances. 

Second, we include a non-Mar kovian decoherence effect by coupling the electron's 
position to a two level system (TLS). This could model a two level fluctuator often found 
in solid state physics [17]. However, our non-Markovian model is very much simplified 
because we do not take into account the coupling of the TLSs to their environment. The 
reason is that full non-Markovian open quantum systems are very complicated. Further- 
more, we are more interested in qualitative differences of Markovian and non-Markovian 
effects, and already with our simple model we find these to be quite pronounced. 

One should mention that there are many other types of decoherence sources, de- 
pending on the physical system CTAP is applied to. E.g. for spin qubits in a solid 
state matrix, the nuclear spins of nearby atoms lead to decoherence |l8j. Although a 
^^Si matrix does not have nuclear spins, a small fraction of ^^Si is unavoidable. Other 
defects might also have nuclear spins. For charge qubits, fluctuating electrical fields are 
a prominent decoherence source |19| . Furthermore, l//-noise is of importance in most 
solid state applications |2Uj . Of course, this work can not take into account all possi- 
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ble decoherence sources, but our conclusions from the case studies of a Markovian and 
a non-Markovian decoherence process should be vahd at least qualitatively for similar 
processes. 

3.1 Quantum dots and quantum point contacts 

Although our work should be usable for many physical systems as indicated in the 
introduction, we used the specific example of electrons on quantum dots. Similarly, 
we could model any Markovian dephasing source, but as a specific example we model 
dephasing due to measurements performed with quantum point contacts. We therefore 
devote this section of the introduction to say a few words about quantum dots and 
quantum point contacts (QPC), as well as two-level fluctuators, which also lead to 
decoherence. 

A quantum dot |21j is mostly made of semiconducting material, but is so small 
(between one and a hundred nanometer), that it is often called a zero dimensional 
quantum system. An electron in a quantum dot is so much confined, that it can move 
in no direction (therefore zero dimensional). The quantum dot acts as a potential 
well for the electron, which can only inhabit the discrete energy levels obtained by 
solving the time independent Schrodinger equation. The smaller the dot, the larger is 
the difference between neighboring energy levels. The spacing between energy levels is 
further increased by semiconductor materials which result in a low effective mass of the 
electrons, typically about ten percent of the free electron mass. 

The quantum properties of a quantum dot are most apparent when the difference 
between energies of different levels are comparable (or even large) to the thermal energy 
of the electron. Therefore, in quantum dot physics, one uses either small temperatures or 
extremely small quantum dots. As the optical properties of any semiconductor depend 
on the energy level structure, one can influence the optical properties of quantum dots 
by engineering their size. 

It is also possible to add an additional electron to a quantum dot. The energy levels 
of this additional electron is also determined by the potential well of the quantum dot, 
that is, by the dots size. If the quantum dot is small enough and if the temperature is 
sufficiently low, then it is possible to ensure that the electron only occupies the ground 
state. This is the realm we will use. 

If two QDs are located close to each other, they can act as a double well, and if the 
barrier between the wells is small enough, electrons can tunnel between the two QDs. 
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The tunneling rate depends crucially on the distance between the quantum dots, but 
as electrons couple strongly to electric fields, the tunneling rate can also be controlled 
externally. For this purpose one can place electrodes between the two QDs and by 
adjusting voltages at these electrodes (often called gates), the electric field and hence 
the tunneling rate can be varied. 

A quantum point contact [22] is a narrow constriction between two wide electri- 
cally conducting regions, of a width comparable to the electronic wavelength (nano- to 
micrometer). At low temperatures and voltages, the conductance of a QPC is quantized 
according to G = NGq. Here, is a natural number which depends on the properties 
of the QPC, and Gq = 2e^//i is the conductance quantum. 

At finite temperatures, there is a steep, but continuous transition between the quan- 
tized conductance plateaus, at which G becomes extremely sensitive to its electrostatic 
environment \22\ I23| . Therefore, if the QPC conductance is set between such a transi- 
tion, it can act as a very sensitive charge detector, able to detect single electrons on a 
QD in its vicinity. 

Two-level fluctuators [T7] are often found in solid state systems and can be a 
major source of decoherence. It is not exactly clear what these are, and there might 
be different physical reasons for their existence for different systems. What is known is 
that they behave like a two level system, randomly fluctuating between the two states. 
Their energy splitting is thought to have a broad distribution and the coupling between 
the TLS and the quantum system of interest also seems to be quite random and might 
differ between consecutive experiments. 
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Chapter 4 

Adiabatic Dark State Population 
Transfer 

The coherent transport of quantum information is an essential element in any scalable 
architecture for a quantum information processor. Much attention has been focused on 
dark state adiabatic passage for coherent state transport. Originally studied in the con- 
text of quantum optics [Ml where it is called stimulated Raman adiabatic passage (STI- 
RAP), dark state transport uses the existence of a "dark state", which is a zero-energy 
eigenstate of a driven quantum system. By manipulating the driving of the system, one 
can sculpt this dark state to coherently transport quantum states using STIRAP-like 
procedures. This intra-atomic dark-state transport has been demonstrated experimen- 
tally [25j. However, the method has more recently been applied to spatial transport 
of quantum information (which we specifically denote CTAP - coherent tunneling by 
adiabatic passage following [12]) in a variety of physical systems, including chains of 
neutral atoms |2B] , quantum dots [121 123 [2H] , superconductors pS] , and photons in 
nearby waveguides |30] . It has also been proposed as a crucial element in the scale 
up to large quantum processors [31]. The method possesses two very crucial benefits 
over other quantum transport methods: since the transport is via a zero energy state 
the quantum state acquires no dynamical phase, and due to the adiabatic theorem, the 
process is very robust to a wide range of system variations. 

4.1 Adiabatic theorem 

This section is devoted to the adiabatic theorem, because the adiabatic state trans- 
fer described in the following section relies heavily on it. This theorem (often called 
adiabatic approximation) [22j is an approximation used in many fields of quantum me- 
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chanics, typically if the coherent dynamics induced by a time dependent Hamiltonian 
involves two different time scales. Its application requires that the time dependence of 
the Hamiltonian H{t) is slow compared to internal time scales, defined by frequencies 
^■mn{t) = {En{t) — Em{t))/h, where Ej{t) are the eigenvalues of the Hamiltonian. The 
adiabatic theorem then states, that the time dependance of the Hamiltonian will not 
induce transitions between different instantaneous eigenstates of the Hamiltonian. 

This situation is commonly achieved when the Hamiltonian can be controlled exter- 
nally. For example one can slowly change the phase or magnitude of a laser field, which 
determines the Hamiltonian of an atom; or slowly change the direction of a magnetic 
field to control particles with a magnetic moment. There are also important applica- 
tions to systems which are not controlled externally, as for example in molecular physics, 
where the Hamiltonian describing the electronic structure depends on the configuration 
of the nucleuses within the molecule. Because the dynamics of the heavy nucleuses 
are slow compared to the electronic dynamics, the adiabatic theorem can be applied 
for the latter. In particular, the adiabatic theorem is a requirement for the use of the 
Born-Oppenheimer approximation in molecular physics |32j . 

The theorem can be stated in a more precise way: 
Adiabatic Theorem: Let H{t) be non-degenerate with H{t) \n{t))= En{t) \n{t)). 

Assume \{ra{t)\ n{t))\ < \Ern{t) - En{t)\/h Vm / n. 
If the initial state of the system is an eigenstate |^(0)) = |n(0)), 
then it will approximately stay in the corresponding eigenstate 
e**^^*^ 1^(0) during the entire evolution. 

Here |ri(t)) denotes the time derivative of the time dependent energy eigenstate \n{t)). 
We will not go into the proof of the adiabatic theorem as it is quite technical and can 
be found in standard text books [_32j. Instead, we will discuss some properties which 
are important for this chapter in more detail. First the reader should note, that the 
assumption in the second line of the theorem is not required for all n, but only for the 
one corresponding to the initial state. 

We also want to say a few words about the scaling of errors due to the adiabatic 
approximation. Under quite general conditions, the undesired populations Pm{t) = 
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\{mit)\tp{t))\ of states |m(t))/ |n(t)) satisfjQ 

This implies that if we want to change an initial state |V'(0)) = \n{0)) adiabatically to 
some final state \ip{tT)) oc |n(tT)), the undesired populations pm for m 7^ n decrease 
with the square of the transfer time tT (because |n(t))cx t^^) 



There are very rare situations when (4.1) and (4.2) do not hold [33], in particular if 



|(m(t)| n{t)) \ oscillates with the frequency But we will not encounter this situation 
in the following. 



More important for our work is that the relation (4.2) can be significantly improved, 
if some higher derivations of |n(t))are small. In our case, all higher derivations are small 
as our tunable parameters are described by Gaussian curves. Then, the final populations 
Pm{tT) decrease exponentially with It [31], which means that exceptional fidelities can 
be achieved with reasonable short transfer times. However, this statement is only correct 
once the change of the state is completed (and H{t > tx) = 0), while the populations 
Pm during the process (for the same reason as in footnote [T]) still follow relation (4.1). 



We will see such a behavior later in our studies, when small populations in undesired 
states |m(t)) appear during the state change, which almost magically disappear once the 
change of state is completed. 



4.2 Coherent population transfer 

We consider a system A described by an A^-dimensional Hilbert space where is an 
odd number, and we assume that we can tune the coupling between neighboring states. 



^This is easily seen from the proof of the adiabatic theorem as given in [32] and is essentially due 
to the fact that the integral of a product of a fast oscillating function (e''^"*"') with a slowly varying 
function (|{m(t)| h{t))\) vanishes in first order. 
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The Hamiltonian for such a system is 



Ha 






Qi 


• 
























































^N-2 














^N~2 













• 











(4.3) 
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and by appropriate tuning of the couphngs we want to adiabatically change the 
state of the system from to \N)a- The couphngs Qi = ^Ip and ^^n-i = ^.re 
often called pump and Stokes pulse, respectively. 

In the quantum optical process of STIRAP (often three-level atoms are used), the 
different levels are energy levels of an atom, and r2„ is the Rabi coupling induced by a 
laser field with frequency One can derive this form of the Hamiltonian [13j if the 
laser field has frequency Un = |e„ — e„_i|/^, where are atomic energy levels; that is, 
if there is no detuning. To do so, one has to use the interaction picture and perform 
the rotating wave approximation. If there is detuning, one has to include entries on the 
diagonals [35], and we will return to this situation in chapter [6| 

In this work, the different levels correspond to the discrete positions of an electron, 
which is allowed to tunnel between quantum dots. We aim to transport the elec- 
tron coherently from the first quantum dot to the last one. The spin of the electron is 
disregarded, which is justified if the tunneling rates Qn do not depend on it. Further- 
more, we assume that the electron populates only the ground states of each quantum 
dot. Quite realistically, the electron can only tunnel between neighboring quantum dots 
leading us to the structure of the Hamiltonian above. We assume that we can tune the 
tunneling rates within < < 0,max by external gates, as shown in Fig. 3.1 In fact, 
during the transport process we will only change the tunneling rates Jli(t) and $77v_i(f) 
between the first two dots and the last two dots, respectively, and keep Qn = ^max for 
n = 2,3, ■ ■ ■ , N — 2. The ket |n)yi denotes the state of the electron when it is found with 
certainty on the n-th quantum dot. 

Of importance for this process is the zero energy eigenstate (un-normalized) [12] of 



the Hamiltonian Eq. (|4.3|) 

|V'o)a 



(JV-l)/2 

cose|l)^ + (-l)(^-i)/2ginG|A^)A-A ^ (-l)i |2j - 1)^ (4.4) 

3=2 
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with 

= arctan — — —, X = 1 N i ^ (4.5) 

i'Af-l o /o2_|_n2 

It is easily seen that |'0o)a = if f^i = and SIat-i / 0, and |V'o)a = \^)a if / 
and VLn-i = 0. Therefore we can use this eigenstate for the adiabatic electron transport 
by performing the following steps: 

1. While keeping ili = one switches — ^max 

for ah n = 2, • • • ,iV - 1. This 

ensures that the Hamiltonian is non-degenerate 

2. Then Vti is increased and ^n-i is decreased until it vanishes. This process has 
to be done slowly to satisfy adiabaticity. It is this step where the electron moves 
from the first to the last quantum dot of the rail. 

3. Finally all couplings can be set to zero. 

In the following, we will refer to these steps as step one, step two and step three. 
We emphasize that the electron moves exclusively in the second step, and therefore, the 
first and the third steps can be done arbitrarily fast (up to experimental limitations). 
Hence, we set t = a.t the beginning of step two and t = ti at the end of step two, and 
call tx = ti — to the transfer time. 

During the entire process, the electron will be in the state |V'o(^)) (up to non-adiabatic 
corrections). It is interesting to note that this state does not populate the even numbered 
QDs, which is very counter intuitive in the following sense: Despite the fact that the 
electron can only tunnel between neighboring QDs, it can move from one end of the 
chain to the other end, without ever populating the even numbered QDs between them. 

This paradox can be understood if one has a look at the non-adiabatic corrections. 
These give amplitudes on the even numbered QDs which are proportional to the in- 
verse of the transfer time tT, therefore validating the intuitive picture that the wave 
function passes through all QDs. The reason we say that even numbered QDs do not 
get populated in the adiabatic limit is that population is the squared absolute value of 
the amplitude, therefore scaling like t^^. This gives population x transfer time oc t^^ , 
resulting in negligible population of the even numbered QDs in the limit of long transfer 
times. 

This described paradox has important implications for STIRAP, as even numbered 

levels are usually excited atomic levels exposed to spontaneous decay. This decay is 

^The Hamiltonian is then non-degenerate only if TV is an odd number. This is the reason for assuming 
odd A'^ in the first place. 
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successfully avoided in STIRAP experiments because these levels do not get populated 
significantly. For our work, this paradox is interesting to mention, but has no impli- 
cations as all QDs are equally subject to decoherence and there is no advantage of 
populating some more than others. 

The real advantage of CTAP compared to non-adiabatic population transfers is that 
the tunneling rates as well as the timing of their switching do not need to be very 
precise (see [26] for a detailed study). This is crucial in many experiments because 
these parameters are often hard to control. The only requirement is that tunneling 
can be suppressed completely when needed, and an approximate control of tunneling 
rates and their timing. As mentioned earlier, the drawback of CTAP is a relatively long 
transport time, which is limited by the adiabatic theorem and usually is about an order 
of magnitude longer compared to diabatic population transfer schemes. 

We conclude this section with an example of CTAP on five quantum dots. The 
tunnel rates 

(i - o.9r)2 



17 1 — ^max GXp 



1^4 — ^max GXp 



(T/4)2 

{t - o.5r)2 



(T/4)2 

^2,3 = ^max (4.6) 



are plotted over time in Fig. 4.1 (a), where the red and black marks on the time axis 
correspond to T = AOh/Qmax and T = QOH/Qmax, respectively. Note that due to the 
Gaussian shape of these rates, there is no clear distinction between step one, two, and 
three. As is the time when all pulses are switched on, we have approximately = 20 
for the red curves and tx = 30 for the black ones. The energy levels are shown in 
Fig. |4.1| (b) and show the required energy spacing during step two. The populations on 



the respective dots are shown in Fig. 4.1 (c) - (g). Comparing the red with the black 
lines, it is clearly seen how a relatively small increase in the transfer time results in a 
dramatic increase of the fidelity (from 0.973 to 0.998). In the limit of slow populations 
transfer there would be zero populations on the second and fourth quantum dot. For 
finite transfer times we find populations on these dots due to non-adiabatic corrections. 



During the transfer these scale according to proportionality (4.2), but most of it vanishes 
towards the end of the transfer (scaling exponentially in t^). This effect, which was 
explained at the end of the previous section, is responsible for the for extremely good 
fidelities (with reasonable times) in the coherent transport of an electron. 
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Figure 4.1: (a): Pump pulse ili/ftmax (dashed), Stokes pulse n4/flmax (dotted), and interme- 
diate couplings fl2,3 /^max (solid). (b): Eigenenergies of the Hamiltonian in units of flmax- (c) 
- (g): Populations on the different quantum dots, starting with the first dot. Red lines corre- 
spond to faster transport, black lines to slow transport. Time is in units of h/Vlmax- Note the 
exponential increase in fidelity with transfer time. 
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Chapter 5 

Modeling of Decoherence 



An important question, particularly with regard to the use of CTAP within large scale 
quantum computer architectures, is to determine the effects of decoherence on the trans- 
port. The effects of dephasing and spontaneous emission has previously been examined 
in the case of STIRAP in a three level atom in a A configuration j36[ I37j. In that work, 
a master equation was postulated and its effects on the population transfer studied. 

We examine the effects of two types of physically-motivated decoherence sources 
effecting the CTAP transport in a quantum dot (QD) chain. We first study the effects 
of delocalised non-readout measurements on the systems making up the QD chain. In 
particular we imagine quantum point contacts (QPC) close to each QD to measure the 
electric charge on the respective QD (see Fig. |5.1| ). These QPCs however, are non-local 
measurement devices in that their charge sensitivity falls off continuously with distance. 
Such devices or similar will be required to either modulate or readout the quantum 
information in a real device. In large scale quantum processors one will routinely wish 
to have quantum information in a superposition of two "distant" spatial location^ 
It is known that from numerous studies of cat-states in quantum Brownian motion 
- a single harmonic oscillator coupled to a bath of harmonic oscillators - the rate of 
decoherence suffered by the cat grows quadratically with the spatial separation of the 
two superposition states of the "cat" We find that such an effect is also present in 
our case, i.e. the decoherence rate of a "cat-state" on the CTAP chain increases with 
cat-separation, but surprisingly we find that this decoherence rate saturates beyond a 
critical cat spatial separation. This is a positive result for the CTAP transport protocol, 
and is essentially due to the rapid spatial fall off of the measurement sensitivity of the 

^If the qubit is encoded in the location itself, then such superpositions arise naturally. But even 
if the qubit is encoded on a different degree of freedom (i.e. spin), such superpositions appear during 
quantum information transport. 



29 



30 



CHAPTER 5. MODELING OF DECOHERENCE 



QPCs 
QDs 



d 



Figure 5.1: Realization of non-local measurement as proposed in [23]. The electric current 
through a QPC is influenced by the appearance of an electronic charge in its near environment. 
Therefore, the electric current on the QPC measures the position of the electron. 

QPCs. 

This first model is an example of Markovian decoherence. We also include a second 
non-Markovian dephasing source, and consider that each quantum dot interacts with 
a nearby two level system (TLS). This could model two-level fluctuators in a solid 
state CTAP scheme. Interestingly, we find that qubit transport still seems relatively 
robust in the presence of these combined Markovian and non-Markovian decoherence 
sources. More worryingly, however, we find that this non-Markovian dephasing slightly 
entangles the two level systems with the transported qubit. Surprisingly, this effect does 
not seriously detract from the transport of an electron isolated on one QD, but causes 
serious degradation if the qubit transported is in a superposition state. As the latter 
situation will be the typical case in a large scale quantum processor, the present analysis 
might indicate a much lower density of TLSs will be required when using CTAP in a 
large scale quantum computer, if operated with charge qubits. The problem is avoided 
by using an internal degree of freedom, like the spin of the electron, as then the transport 
of spatial superpositions is not required. 

5.1 Measurements as decoherence 

As we mentioned above, we first consider each QD of the QD-chain to be continuously 
measured in a manner that gives rise to Markovian decoherence. We consider the mea- 



surements to be made by QPC situated close to the quantum dots as shown in Fig. 5.1 In 
particular, we imagine that each time an electron travels through the QPC, this electron 
weakly measures whether there is an electron situated on a nearby quantum dot. The 
measurement executed by a QPC is caused by the modulation of the conductivity of the 



QPC due to the presence of a nearby electron (see section 3.1). The QPC conductivity 
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is modulated by a factor 1 — a/r, where r is the distance between electron and QPC and 
a is a constant reflecting the properties of the QPCs (see ^3]). This modulation results 
in an indirect position measurement of the electron's spatial position on the rail of QDs. 
However, it is a non-local measurement because even an electron on the neighboring QD 
influences the transmission through a QPC. The localness is parameterized by a/d, i.e. 
the distance between QPC rail and QD rail over the distance between two neighboring 
QDs, with small values representing more local measurements. Furthermore a/ a pa- 
rameterizes the sensitivity (signal over noise) of the measurements and is typically small 
and hence the measurements are weak ones. Such measurements are properly described 
in the language of positive operator valued measurements (POVM) jHE]. 

The purpose of a measurement apparatus is the readout of quantum information 
for which one would like to use strong local measurements. This can be approximated 
by using a large number of weak, non-local measurements of the type described. A 
large number of measurements in a reasonably short time is achieved by using a high 
measurement rate which in turn can be realized, for instance, by applying a voltage to 
the QPCs. In this case there will be a small but macroscopic current through each QPC. 
By measuring this current, an observer learns the precise distance between the electron 
and the QPC, and essentially performs a perfect von Neumann-Liiders measurement. 
However, such measurements also act as a strong source of decoherence. During quantum 
unitary operations such as transportation by CTAP it is preferred that this decoherence 
is absent. However, switching off the measurements might not be completely achievable 
in practice and one may be left with a small current through the QPCs, either caused by 
non-zero voltage (possibly fluctuating) , or by electrons which travel thermally across the 
QPC. Even if this current might be too small to be detected (no read out possible), each 
electron which travels across the QPC causes decoherence which is described by non- 
readout measurement^ Thus, non-readout measurements as a source of decoherence 
are included in the analysis presented here. 

As in the previous chapter, we restrict our treatment to the case of having only 
one electron in the rail of QDs. We also assume that the electron can only occupy the 
ground state of the QDs, and we take \i)A, to be the quantum state of the electron in the 
i*^ QD. Furthermore we neglect all interactions depending on the spin of the electron. 
Then, {\i)A, i = ^, --i A^}, form a basis for the Hilbert space Ha of this electron on the 
QD rail with dots. In the following we take the limit of a long rail, A^ — )• oo. We 

^We are not interested in a possible readout during the transport anyway, as in a realistic setup we 
could not undo the decoherence effects of a measurement. 
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denote the distance between the i-th QPC and the j-th QD by 

nj = Va2 + (|i-j|d)2 (5.1) 



where a and d are defined in Fig 5.1 

The probabihty of the j^^ QPC detecting the presence of an electron on the QD rail 
can be written as [1] 

P,(p^) = Tr(^,/) (5.2) 

where is the state of the electron on the QD rail and itj is the effect operator 
corresponding to the QPC measurement at site j. If the electron is spatially localised 



to be only on the i^^ QD, i.e. in the state \i)A, Eq- (5.2) reduces to 



Pji\^)A{i\) = A{'i'\^j\^)A- (5.3) 

As we noted above the measurement sensitivity of the QPC depends on the distance 
rij. The presence of an electron a distance Vij away from the QPC decreases the current 
flowing through the QPC by a factor 1 — a/rij and this leads to a reduced detection 
probability, 

P,(|i)^(i|)ocl-^ . (5.4) 



Fulfillment of Eq. (5.4) is certainly achieved with the effect operators 



N 



N 

i=l 



which can be checked by substitution into Eq. (5.3). The constant 



N 1 

k{N) = -. ^ = ^ ; (5.6) 



N 

a 



n = lim RiN) = 1 + lim V — — (5.7) 

j=i J 

is chosen to satisfy Xlj^i ~ ^- Note that each effect operator is almost proportional 
to the unit operator (remember that a /rij is small) which reflects the weakness of the 
measurements being performed. 

We will assume that the measurements performed by the QPCs are efficienlj^ (see sec- 
tion 2.4), that is they only introduce a minimum amount of decoherence. Measurement 



^This is quite realistic as we describe the measurement from a very microscopic viewpoint 
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theory states that the transformation of the density operator due to such a measurement 
with result j is described by [l] 



(5.8) 



with Aj = Uj ^J TTj for some arbitrary unitary operators Uj . This unitary depends on the 
interaction of the quantum system and measurement apparatus during the measurement 
process. In our case, if the electron is at site i, its position will not be changed by a 
detection event. For simplicity we also assume that a measurement does not introduce 
relative phases within the rail of QDs, leading to Uj = 1 and A = A^ , which means the 
measurements influence the electron state as little as possible: 



A, 



N 

E 



1 



(5.9) 



1=1 V 

To derive the master equation we now assume that detection events in the QPCs 
occur uniformly at random and at a constant rate R. Furthermore, we ignore the 
detection result and average over all possibilities (non-readout measurements). Following 
|10j . we can write down the master equation describing the evolution of the density 
operator as 



d/ 
dt 



N 



^ [Ha, /] -Rp^ + rY, A.p^Aj . 



(5.10) 



This equation is in Lindblad form with Lindblad operators Aj and thus the evolution 
is Markovian. The two incoherent terms can be interpreted in the following way: Each 
time a measurement occurs, the state p"^ is substituted by ^jP^^j corresponding to 
a non-readout measurement. This process happens with rate R. We note that R should 
scale proportional to A'^ as each QPC contributes equally to the overall measurement 
rate {R/N is the average measurement rate of a single QPC). For the case when Ha 
involves no electron transport along the QD rail, Eqn. Eq. ( 5.10| ) possesses a stationary 
state which is diagonal in and thus the evolution corresponds to a pure dephasing 
type of decoherence. 

For now we take Ha = 0, but later we introduce a time dependent Hamiltonian to 



induce CTAP. Expressing Eq. (5.10) in the basis \k)A and using Eq. (5.9), we obtain 



Kk ~ 0' ^"-"^ diagonal entries and p"^i 
rate is 



-Dkipti for off-diagonal ones. The dephasing 



Dki = R 



1 



oo 



N ^ 

j=~co 



a 

Tkj 



a 

no 



(5.11) 
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where rij is defined in Eq. (5.1). 



5.1.1 Properties of the dephasing rate 



In this brief subsection we study the dephasing rate Eq. (5.11) in more detail. First we 



show that it saturates for large distances \k — l\d ^ a. To this end we assume / > k and 
split the sum into two parts 



AT 2^ 



N 



1 



a 
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a 



N ^ 

j=-oo 



1 



+ 1^ 
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i=(fc+0/2 



1 



a 



(5.12) 



where we also used that r/j is large in the first sum on the right hand side of the equation, 
and rkj is large in the second sum. With the same argument we can extend the sums 
to infinity 



(fc+0/2 
K X - 



N 



a 

Tkj 



N 2^ 



N 
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Tkj 
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N 



E 1 

j=(k+l)/2 



(5.13) 



to find 



K, 

N 



1 



a 



a 



1 



(5.14) 



j=— oo ' ' J=— oc 

Due to the periodic distribution of the QDs and QPCs, this is clearly independent of k 
and / and is therefore the limiting dephasing rate for large separations |A; — ^|. 

Next we consider the limit of weak measurements 
be well justified in experiments 



<C 1 in Eq. (5.11), which should 



This weak measurement limit corresponds with the 
inability of the measurements to give detailed information on the position of the electron 
in the QD rail. Weak measurements of this nature feature in many models of continuous 
monitoring of a quantum particle's position [38] . In the limit of weak measurements we 



can expand the square roots in Eq. (5.11) to second order in a /Rij. Exploiting the 



D 



kl 



Ra^ ^ ' ^ ^ ^2 
8N 



periodicity of the setup we then use Y2j ^/'''kj = l/^j and applying Eq. (5.7) we find 

(5.15) 



E 



1 

rkj 



1 



where it is again apparent that R has to scale linear with as already discussed from 
a physical point of view. For spatial separations larger than the threshold, \k — l\d 3> a, 
we can again split the sum into two parts, extend them to infinity and find 
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2 ^ 



{a/dY + {k-jf 



2 ^ 1 



E 



(i2 [a/df+p 

j=-oo 

2tt coth {na/d) 



(5.16) 



d2 a/d 

The last equality is a standard formula which can be checked i.e. with Mathematica. 
Finally we get the large distance saturation value of the dephasing rate 

„ irRa^ coth (it a/ d) 

lim Dki = -1--^ \ / ' , (5.17) 

|fc-/|->oo ANd^ a/d ' ^ ^ 

valid in the weak measurement limit a/ a <^ 1. 

We also note that we can model local measurements, where the measurement result 
of any QPC only depends on the charge of the nearest QD, in the limit = 6ij, while 



keeping a finite, e.g. in Eq. (5.11). In this case we find from Eq. (5.11) 



N 

which is of course independent of k and / 



D = ^{2 + a-2VlT^), (5.18) 



The saturation of the decoherence rate in Eq. (5.17), is somewhat surprising when one 
compares this with the similar situation for a free particle in a spatial superposition cat- 
state, experiencing continuous position measurements [2]. In that case the decoherence 
rate suffered by the particle increases without bound according to the spatial separation 
of the cat, i.e. Dx^^x2 ^ ixi—X2)'^. This however depends very much on the measurement 
model. 

5.2 Coupling to two level systems 

We now consider a further source of decoherence. It is highly likely that in any physical 
device there will be unknown accidental two level fluctuators nearby to the quantum dot 



rail (see Fig. 5.2 ). In fact, experiments in solid state physics [T7] (and references therein) 
often find decoherence due to coupling to two level fluctuators. Although not much is 
known about these two level fluctuators, they certainly have to be taken into account 
in many solid state devices. In fact, from the study of superconducting qubits [39j, it 
seems likely that these TLF are within any layer of amorphous silicon oxide. Such layers 
are inevitable in most silicon devices, and would also be needed in the implementation 
of CTAP using phosphorus atoms embedded in a silicon matrix |14| . 
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Figure 5.2: As in Fig 5.1 but with additional local coupling to two-level systems. 



As a rail of QD's would most likely be embedded in a solid state matrix in any 
technical device, we include these mysterious systems in our studies. If these unknown 
two level systems (TLS) can couple to the electron on the rail, then these systems act 
as a source of decoherence which exhibits memory, i.e. is non-Markovian. That is, 
quantum coherences on the quantum dot rail can be transferred to the nearby TLSs, 
where they can remain for a period, before being transferred back. Typically the analysis 
of these types of non-Markovian effects are complex but in the following we are able to 
derive analytic solutions of the resulting reduced dynamics of the quantum dot rail. 
For simplicity we assume these fluctuators have no internal dynamics other than their 
coupling to the quantum dot raiQ and that this coupling is local to the nearest QD 
only. Furthermore we assume that these TLSs do not experience significant decoherence 
on the time scales of the transport. We are aware that these assumptions may not be 
completely satisfied in all realistic situations. However this model will serve to highlight 
the striking difference between the Markovian and non-Markovian evolutions in, for 
instance, the transport of a spatial superposition. 

We use and |0)j as basis of the Hilbert space Tij of the j— th TLS such that the 
interaction Hamiltonian is diagonal. If the electron is on the j-th QD, it is assumed to 
induce a phase shift on the j— th TLS, so that the Hamiltonian acting in the product 



Hilbert space Ti = Ha ^jLi T^j of electron and TLSs reads 



int 



Xn \n)A{n\ (S> (T^,„ 



(5.19) 



The coupling constants Xj ^re considered to be constant in time, and a^.j and Ij are 



*A possible internal Hamiltonian Htls does not influence the QD rail if [Htls, Hi„t] = 0, because 
it can be removed by the unitary operator e*^^^^''''', without changing Hint- 
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the Pauli Z-matrix and the identity operator, respectively, acting in Tij. 

As typically assumed, we now take the initial state to be in product form Q{t = 0) = 
Po ^ Po^^^ ^ where Pq^^^ does not have to be a product state of the different TLSs. We 
can now express the master equation describing the time evolution of the density matrix 
of the combined QD rail and TLSs under the effects of the above measurements (see 
Eq. dslol )), to be 



d ^ 

-£ = -i[Hinu q]-Rq + RY, ^oQ^o (5-20) 

with Aj = Y^(8) i^^- 

After some effort, one can trace out the TLSs to find the non-Markovian master 
equation for the reduced density matrix of the QD rail for an arbitrary initial product 
state of the TLSs, as 



d/ 
dt 



N 



-Rp^ + [R^nP^An - ^A„(^) [\n)A{n\ , /] 

n=l 

+7„(t) [\n)A{n\,[p^,\n)A{n\]] }, (5.21) 



where the last two terms describe the effects of the TLSs on the QD rail. Here the 
definition 

/,x . s'mXnt-tUJnCOSXnt , . 

-fn[t) - ^A„ t) = Xn — -. 7 5.22 

COS Xnt + iUJn Sm Xnt 

is used and a;„ =Tr[p„(5"2.„], is the inversion of the n-th TLS. Decoherence due to the 
TLSs is described by 7„ whereas represents the Lamb-shift. If ujn = ±1, then 
A„ = ±Xn and 7„ = 0. In this case the decoherence due to the coupling to TLSs 
vanishes and the coupling to the TLSs results only in a change of the energy of states 
|n)yi by ±x„. However, a more interesting case is when uJn = 0, since it includes the 
case where the TLSs may initially be in a complete mixture Pn{0) = ^l- In this case 
we find A„ = 0, and 7„ = Xntanxnt. The presence of singularities in 7n(t), may cause 
difficulties in studying the time evolution of p"^, using normal methods. To avoid this 



we do not work directly with (5.21). Instead we return to solve the complete dynamics 
of the coupled QDs and TLSs and then trace out the latter to obtain p^{t). 



To this end, we transform Eq. (5.20) with 

N 



\n)A{n\ (cos{xnt)tn + ^sin(xr^^)(5■2,n) (^ij- (5.23) 

n=l j^n 
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The combined density operator in this picture is given by Q{t) = e*^™'*^(t)e and 
is governed by the master equation 

d^ ^ 



dt 



-R0 + R'^AjgAj. 
j=l 



(5.24) 



Note that e*^*"**Ae ^^i^tt _ i Taking the components of Eq. (5.24) in the QD Hilbert 



space, Qf:i{t) = A{k\ g\l)A (which is stiU an operator in the TLSs Hilbert space), we can 
repeat the calculation of the previous section to find 



Qkiit) = = e-^^'*a/(0) 



Dkit . 



(5.25) 



where D^i is given in Eq. (5.11). After transforming back to the Schrodinger picture 



and tracing over the TLSs, we find for the components of the reduced density operator 



put) 



-Dkit 



cos(xfct) COS {xit)Pki{^), 



(5.26) 
(5.27) 



where we have assumed that the initial states of the TLSs is a completely mixed state 

1. 



-1. 



(5.28) 



Hence we see that information lost to the TLSs can return to the rail via the oscillatory 



terms in Eq. (5.27), which is in contrast to Markovian decoherence induced by the 



measurements. Note that the non-Markovian behavior of Eq. (5.27) results from tracing 
out the environmental TLSs. 



Chapter 6 

Coherent Transport 



Here we employ the technique of CTAP described in chapter [4] to the system of QDs 
coupled to QPCs as well as to TLSs as outlined in chapter [5] A minor difference to 
chapter |4] is that the system here consists of an extended rail of QDs with an infinite 
dimensional Hilbert space %a- The aim is to transport an electron coherently from the 
m-th QD to the n-th one, where the number (n — m + 1) of involved QDs is required to 



be odd because of reasons discussed in section 4.2 Although we can restrict the system 
Hamiltonian Ha as well as the local coupling to Hamiltonian Htls to the subspace 
corresponding to these QDs, we have to take into account measurements performed by 
QPCs on sites other than m, • • • , n. This is because of the non-local nature of these 
measurements . 

The system Hamiltonian reads 



/ 



Ha 



m 




m 








V 



\ 








0„_2 
^n~l / 



(6.1) 



where following section 4.2 we call VL 

0. 



Stokes pulse, and set il,- 



Vlp{t) the pump pulse and Q.n-i = ^s{t) the 
for all j = m + 1, • • • , n — 2. Also, we will use the 

to at 



division of the process into three steps as introduced in section 4.2, and set t 
the beginning of step two, just before the population transfer starts, as well as t = ti at 
the end of step two, when the population transfer is finished. 
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6.1 Measurements 



To introduce our technique for solving a master equation in the adiabatic approximation, 



we first neglect the coupling to the TLSs. The master equation to be solved is Eq. (5.10) 



dt 



1 



[Ha,p^]-Rp^ + R ^jP^^j- 



(6.2) 



with the system Hamiltonian Eq. (6.1 ). At t = to the (n — m+l) eigenstates of Ha which 
are superpositions of \m)A, • • • , \n)A, are denoted by \ipj{to)) with j = ^^^^y^, • • • , 
ordered by their energy Ej{to). These eigenstates evolve continuously to \ipj{t)) which 
depend on the choice of the pump and Stokes pulses Qp{t) and 0,s{t), respectively. The 



un-normalized adiabatic state responsible for the transport is according to Eq. (|4.4|) 

\Mt)) 



n — m 
2 



COS0 \m)A + (— 1) 



^ sine |n)A - X (-l)J' |2j - 2 + m)A 
i=2 



(6.3) 



with 



= arctan 



X 



pi '5 



nr. 



p ~\~ g 



Note that \tpo{to)) = \m)A and |'0o(*i)) = \n-)A, i-e. the states in which the electron 
is on the m— th and n— th QD, respectively. Furthermore during the entire process the 
energy EQ{t) = of \ipQ{t)) vanishes, which ensures that no dynamic phase appears for 



the state to be transported. We recall from section 4.1 that if 

d 



\Eo{t) - Ej{t)\ > 



dt 



(6.4) 



holds and if the system at t = to is in |V'o(io))) then the adiabatic theorem states that 
the system will stay in \ipo{t)), provided it is a closed system. Therefore an electron 
starting in \n)A will end up in \m)A- 

To generalize this concept to the open system described here, we follow [ID] and 



transform Eq. (5.10) with the unitary operators defined by 

UHt)\i;,{t)) = |^,(to)) forj = ^,. 

U\t)\j)A = \j)A forj<mandi>n 



n—m 
2 



(6.5) 



to get 

dp^ 
dt 



—I 



n — m 
2 



Ej{t)\i;j{to)){i;,{to) 



rrf'^U ^A 



N 



Rp^ + RY,^iP^Ai (6.6) 



1=1 
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with O = U^OU for any operator O. If, in addition to adiabaticity Eq. (6.4) we also 
assume weak coupling to the environment one can neglect the term ^f^^^ in Eq. (6.6) 



as is shown in |40] . This is the generalization of the adiabatic theorem to systems 
described by a master equation of Lindblad form. Hence we have achieved the time 
independence of the eigenspaces of the transformed Hamiltonian H: 



dt 



n — m 
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E 

m — n 

~ 2 



Ej{t)\ijj{tomjito)\,p' 



N 

Rp^ + R'^Aip^Ai 

i=l 



(6.7) 



For R = we get the von Neumann equation 
dp^ 



dt 



n — m 
2 



Ej\i^,{to)){i^,{to) 



(6.^ 



and P(x)(t) := {4'o{to)\ P^it) |V'o(*o)) = A{fn\ P^{to) \'>t^)a is then a constant of motion, 
which ensures perfect transport of the electron. At t = ti we use U{ti) to transform 



back to the Schrodinger picture. From Eq. (6.5) we find again that an electron initially 
in state \m)A = \ipo{tQ)) will be in state \ipo{ti))= at the end of step two. 



For R we use Eq. (6.7) to calculate the loss from perfect transport. To this 
end we note that the initial state is |'i/'o(^o))) and in the transformed picture p^ we aim 
to stay in this state. Therefore the loss of transfer fidelity due to the measurements 



increases in time according to the projection of Eq. (6.7) onto |^o(^o)): 
d 



dt 



4o{t) = -R{Mto)\p^(.t)\Mto))+R {Mto)\Ajp^{t)Aj\Mto)) (6.9) 



As we argued before, without measurements p^{t) is a constant in time in the adiabatic 
approximation. To first order in the measurements]^ we can therefore substitute p^{t) — )■ 
p^(to) = \Mto)){Mio)\ and p^{t) |^o(t))(V'o(i)| on the right side to find 



±7A 
dt 



pUt) 
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-R + R 



oo 

E 

j=-oo 
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(^o(to)MjlV'o(io))' 



-R 




-\mt)\i)A\'\ 



(6.10) 



^This is usually justified in systems described by a Markovian master equation. 

^That is for small measurement rates R/N and/or for weak measurements, i.e. small a/a. 
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where we used the definition Eq. (5.9) of the measurement operators Aj. Now, \tpQ{t)) 



can be substituted from Eq. (6.3) to obtain numerical values and the error probabihty 



of the electron transport is obtained by integrating Eq. (6.10) 

Error Probability = - [ dt^^^^. 



(6.11) 



One might also ask what happens to coherences Pkoi^) •= ^(^1 P'^it) iV'o(^o)) (for any 
k < m or k > n) during the transport. This question arises if one is interested in the 
transport of one part of a super position state, e.g. (ci \k)A + C2 \m)A) — >■ (ci \k)A + 



C2\i^)a)- In the same manner as Eq. (6.9) and Eq. (6.10) we arrive at 
d 



dt 



-R{k\p^{t)\Mto))+R J2 {k\Ajp^{t)Aj\,Po{to)) 

j=-oo 

oo 

-Rcicl + R cic*2{k\Aj\k){Mto)\Aj\Mto)) 

j=-oo 
oo 

-Rcicl + R ci4{k\Aj\k){Mt)\Aj\Mt)) 

-RciC2 



(6.12) 



J = -CO 



j=— oo ' t=m ' 



'o{t)\i)A\ 



(6.13) 



where again Eq. (6.3) can be substituted. In the line with the approximation sign, we 
substituted p^{t) with p^{to) on the right hand side of the equation which is again valid 
to first order in the measurements. Furthermore, we only kept cid^ |^)('0o(^o)| as the 
other terms of /?^(to) vanish once the Ajs are substituted in the last line. 



The terms in the square brackets of Eq. (6.10) and Eq. (6.13) become small when the 



measurements are sufficiently non-local, such that they can not distinguish well between 



the QDs involved in the transport (see also Fig 6.1). Since the loss of information 
increases linearly in the time ti required for transport, it is crucial that the couplings 
between the dots are as large as experimentally possible to give a big energy splitting 



which in turn allows a fast transport (see Eq. (6.4)) 



Also note that in the case of local measurements Eq. (|6.13|) reduces to ^P^q 



-Dp^ 



kO 



where the decoherence rate D = ^ (2 + a — 2\/l + a) is the same as the one obtained 



in Eq. (5.18) without transport. Therefore, the decoherence rate of a charge qubit on a 



quantum dot rail is the same during storage as during transport by CTAP if subject to 
local measurements. 
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Figure 6.1: Transfer loss during the transport of an electron along three (black), five (blue), 
seven (red), nine (green), and eleven (purple) QDs as a function of a/d, i.e. the localness of 
the measurements. The crosses show the numerical solution of the exact equation Eq. (6.2) (i.e. 



before the adiabatic approximation is performed) for local measurments and the circles show 
the empirical cross-over between local and non-local measurements (see text). Pump and Stokes 
pulses are as in Eq. (6.17) (see also Fig |6.2[ a)) with TQ„iax/fi chosen such that the additional 
loss due to non-adiabaticity is 24 x lO^^for all curves, to allow reasonable comparison. This 
is T = 150, 196, 225, 242, 249 for three, five, seven, nine, and eleven QDs, respectively. The 
intermediate couplings are equal one. System parameters are R = N, a = 0.04a. 
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The probability of not finding the electron in the desired state \n)A after the trans- 



port, calculated from Eq. (6.10) and Eq. (6.11), is shown in Fig 6.1 as a function of the 
parameter a/d for transport along 3, 5, 7, 9, and 11 dots. Remember that a/d is the 
distance of the QPCs to their nearest QD compared to the distance between neighboring 
QDs, and therefore is a measure of the non-localness of the measurements performed 
by the QPCs. For any given a this parameter can be changed with d and we fix the 
sensitivity vTj of the measurements by setting a = 0.04a. Also the rail of dots 
and measurement apparatuses extends on both sides {N ^ m — n). In the special case of 



local measurements we also numerically integrated the exact master equation Eq. (6.2) 



(i.e. not applying the adiabatic approximation) and the crosses in Fig 6.1 are from 
this solution. The great agreement with the curves indicate a good performance of the 
extension of the adiabatic approximation to open systems. 



One can distinguish two regimes in Fig 6.1 First, for ^ > " (right of the circles) 



the transport fidelity increases with the non-localness of the measurements. This is 
easily understood since sufficiently non-local measurements can not distinguish between 
the QDs involved in the transport, therefore not leading to decoherence. 

Second, for ^ < "4"^ (left of the circles) we find the reverse and the transport 
fidelity increases with the localness of the measurements (except for transport along 
three QDs). This needs some explanation. In this regime the measurements are local 
enough to distinguish well between the states |m)^ and |n)^ (during CTAP most of 
the population is found on these two states) and hence one can not expect a further 
decrease of the fidelity with the localness of the measurements. To account for the 
apparent increase of the fidelity with localness we need another argument. For very 
local measurements, almost exclusively measurement results from QPCs at sites n and 
m contribute to decoherence, whereas if the measurements are slightly non-local, QPCs 
in the near neighborhood of n and m also contribute, therefore increasing the overall 
decoherence. 

The small decrease of the fidelity for extremely local measurements might be due to 
the small populations on QDs other than m and n. 

It is also evident from Fig |6.1[ that if local dephasing is the main source of deco- 
herence, it is best to transport long distances in one step as is seen e.g. for the transfer 



|1)a — >• |11)a- For the parameters given in Fig 6.1 , the chance of not finding the electron 
on the desired state after the transport is then about 11 * 10~^ (see left end of Fig 6.1). 
This number increases to about 5* (4* 10~^) when we first transfer to |3)a, then to |5)a, 
|7)a5 |9)a; and finally to |11)a- On the contrary, if dephasing tends to be more global. 
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a better transfer rate is achieved by breaking a long distance into several smaller ones 



(see right end of Fig 6.1). 



6.2 Coupling to two level systems 



We now apply CTAP to QDs coupled to TLSs as shown in Fig 5.2, Hence the driving 
field Ha from Eq. (6.1) has to be added to Hint from Eq. (5.19), responsible for the 



QD-TLS-coupling. We first neglect the coupling to the measurement apparatus and 
point out its inclusion at the end of the subsection. Then, all couplings behave local 
and we can restrict ourself to study only the QDs of the rail which are involved in the 
electron transport, i.e. N = n — m + 1. 

To understand how CTAP works for this system, we note that the Hamiltonian 
H = Hint + Ha can be written in a block diagonal forrrj^ with 2^ blocks, each block 
having dimension N. Each block gives the Hamiltonian of the electron on the rail of QDs, 
conditioned on the state of the TLSs. As an example we show part of the Hamiltonian 
for = 3 explicitly 



H 



-XI 




-X2 
^2 





-X3 



-Xi 




-X2 
VL2 





X3 



-Xi 




X2 
^2 





-X3 



V 



where the first block is for all TLSs in the |0) state, the second for the first and second 
TLSs in |0) and the third TLS in |1), and so on. 

We can now study each block individually because there is no coupling between them. 
If we want to transport the electron from the first QD to the third QD, regardless in 
which state or superposition the TLSs are, we have to make sure that CTAP works in 
each block. That is the state |l)yi |0)|0)|0) should adiabatically evolve to |3)a |0)|0)|0) 

Hamiltonian H — Ha + Htls + Hint with [Htls, Hint] = always allows such a block diagonal 
form, because then Htls does not induce transitions between different eigenstates of Hint- 
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as well as |0)|0)|1) should evolve to |3)a |0)|0)|1) and so on. If this is achieved, 
then an electron initially found in \1)a will evolve in the pure state \3)a. During the 
population transfer the reduced state =Tvtlss[q] is not only in a superposition of 
|2)a and |3)a, but in a real mixture of these states (see Fig |6.3[ a)). This is 
because depending on the state of the TLSs, the electron will start moving earlier or 
later, which generally results in entaglement between the rail of QDs and the TLSs. 
However, once CTAP is completed the electron will be found in the desired state \3)a, 
and no entanglement will be left. This reduction of entanglement at certain times is a 
common feature of non-Mar kovian environments, where transfer of information between 
system and environment is possible in both directions [41j. 

Hence, for CTAP to be applicable despite coupling to TLSs, we only have to ensure 
that it works for a Hamiltonian of the form (where we return to an arbitrary long rail 
of dots) 



H 



np ±x 



\ 








±Xm+2 






\ 










±Xn-l 



(6.14) 



instead of the simpler form Eq. (6.1 ). CTAP should work for all possible permutations of 
+ and — , since each permutation represents one block of the Hamiltonian corresponding 
to a certain combination of environmental TLS states. The crucial requirement for this 



is again relation (6.4), but the energy eigenstates \iljj{t))as well as the energy eigenvalues 
Ej(t) depend now on the couplings Xm, ■ ■ ■ ,Xn and the state of the TLSs. An important 
feature which is due to the non-zero values on the diagonals is that this relation can not 
always be fulfilled by just performing the transport sufficiently slowly. Additionally one 
has to ensure that no level crossings occur during step two, i.e. Eo^t) — Ej{t) ^ 0. As 
before, -E'o(^) is the energy of the state \^o{t)) which we want to evolve adiabatically from 
|V'o(^o))= \it^)a to |V'o(*i))= |'^)a- But unlike in the previous subsection, EQ{t) ^ and 
depends on the state of the TLSs. The situation becomes even more complicated when 
anti-crossings are considered which appear quite naturally (see blue and red curves in 



Fig 6.2). At an anti crossing two energy levels can get so close that adiabaticity can not 



be achieved with any reasonable transfer times. 
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Furthermore it might happen that the adiabatic state |V'o(*)) does not connect \m)A 
to \n)A, but instead to another state |'0o(ti)) = \j)A with j ^ n. An example of this is 
shown in Fig |6.2| (green curves) where the population returns to the original QD. 

To calculate the energies is generally only possible numerically and an exact treat- 
ment when level crossings or anti-crossings appear is highly non-trivial. Some work on 
this issue is done in |35| . but the possibility of anti-crossings is not considered there (see 
also [l2])- However, some qualitative statements can be made. If the energy level EQ{t) 
of the adiabatic state is at all times during step two in the center of all energy levels 
involved in the transport, i.e. 

Tfl — Tl 

Eo{t)>Ej{t) for j = __,..., -1 

Eo{t)<Ej{t) for j = l,...,^— , (6.15) 
then we can be sure that no level crossings occur. In this case the process during step 



two is qualitatively the same as without coupling to the TLSs. Since condition (6.15) 
is valid with Xj = (previous subsection), it can always be achieved with sufficiently 
high driving fields ^Ij , as then the comparatively small couplings Xj do not influence the 
qualitative eigenvalue structure. 

Lets apply this condition to the block of the Hamiltonian in which we have minus 
for Xm and pluses for all other coupling constants. Clearly we need sufficiently strong 
tunneling rates Umax and Qs{to) to ensure that £"0(^0) = —Xm is not the lowest energy 



at the beginning of the transport, but at the centre of all energies (see Fig 6.2 c) where 
£"0 gets to the centre of the energy levels at t ^ 3 before the transport starts at to 20). 
The key feature to realize is that level crossings or anti-crossings are actually needed to 
get Eq to be at the centre of all energies. But it is crucial that this happens as long as 
rjp ~ 0, i. e. before the population transfer starts. 

For n—m = 2 we find ^2^(to) > (xi~X3)(Xi~X2) for j / P. For transport along many 
QDs {n — m ^ 2) we get the approximate inequality ^j{to) ^ cmax(xfc)(n — m) for some 
constant c, which states that for given maximal Qj (due to experimental limitations) 
and given Xj^ there exists a maximum number n — m for which we can use CTAP. For 
transport along more QDs one has to break up the transport into several smaller sections 
and apply CTAP successively to each section. The shorter the rail gets (compared to 
the maximal length), the further we can stay away from any level crossing and the faster 



we can transport the electron without violating the adiabaticity condition (6.4). 



As example we show in Fig 6.2 the eigenenergies and populations as a function of 
time for a system with five QDs and constant Xi = X for different values of x- We 
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consider the block of the Hamiltonian acting in the subspace where the first TLS is in 
the state |0) and the other four TLSs are in the state |1), which reads 



Ha 



( 


-X 


rip 








\ 




np 


+x 


^max 













^max 


+x 


^max 













^max 


+x 




\ 













+x ) 



(6.16) 



Of course, if we assume the TLSs to be in a statistical mixture one would have to 
calculate the populations for all blocks of the Hamiltonian and average over them. Since 
all qualitative features can be seen with this particular block we restrict the discussion 
to this one for pedagogic reasons. We use Gaussian pump and Stokes pulses 



— ^max exp 



{t - 3T/4)^ 
(r/4)2 



^5 = ^maxexpj-^^^^^^l (6.17) 
with the parameter T = 150. Because of the Gaussian pulse shapes the distinction of 



steps one, two and three is not sharp. Looking at Fig 6.2 (a) we can approximately say 
that step one is for —50 < t < 20, step two is for 20 < t < 130 when all couplings are 
switched on, and step three is for 130 < t < 200. T is chosen long enough to achieve 
very good transfer fidelity (0.9996) for X = 0) without coupling to TLSs. 

For X = 0.15 there appears a level anti-crossing at which the energy gap (0.0002) is so 
small, that in can be treated as a crossing on the time scales set by T, the Gaussian pulse 
width. If one were to enlarge T extremely, adiabaticity would be fulfilled even for this 
tiny energy gap which would result in no transport, similar to the case x = 0.45 described 
below. However, as fip ~ at the time of the anti-crossing, the state will remain in 
its initial state |m)with Em = —0.15 during this anti-crossing for any reasonable pulse 
width T. That is, the state crosses the levels as indicated by the dashed line in the 



magnified part of Fig. 6.2 (c). The smaller population transfer (fidelity 0.975) is because 
during step two, the energy level spacing is decreased compared to the case with x = 0- 
Therefore, the transfer fidelity can be increased arbitrarily close to one by choosing 
larger 

Enlarging x = 0.3 moves the anti-crossing into step two and adiabaticity can not 
be achieved any more. The result is wild oscillations of all populations and a low 



Non-adiabatic corrections scale with exp {~cT\Eo — Ej\) for some constant c (see also section 4.1 1. 
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Figure 6.2: (a): Tunneling rates (dasiied: fli, dotted: ^2,3, solid: ^4), (b)-(e): energy levels of 



Hamiltonian Eq. (6.16) and x = 0; 0.15, 0.3, 0.45 (see text), and (f)-(j): populations of \1)a 
- \5)a obtained by the numerical integration of the Schrodinger equation. The colors (see (j) 
for grey scale) black, blue, red, and green correspond to x = 0, 0.15, 0.3, 0.45, respectively. 
Time is measured in units of f^.^ix- What looks like level crossings in (c) and (d) are actually 
anti-crossings as can be seen in the magnified figures. 
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transfer fidelity. Finally at x = 0-45 the anti-crossing disappears and adiabaticity is 
restored which is evident in the lack of oscillations in the populations. But because the 
eigenenergy -0.45 at t = — oo connects to the same energy at t = c«, the population 
returns to the original state \m)A and the population transfer failed. 

This result should be quite surprising. Up to relatively large couplings to TLSs, 
— < 0.2 (for m — n = 4) adiabatic transfer can be achieved and the transfer fidelity 
very quickly approaches one with sufficiently large transfer times. This is contrary to 
the Markovian dephasing studied in the previous subsection, where the transfer loss 
increases with transfer time once the transfer time is long enough to ensure adiabaticity. 

There is however a disturbing effect due to the coupling to TLSs is that £'o / and 
hence we get a dynamical phase 

ip= f dtEo{t) (6.18) 



and worse even, this phase depends highly on the state of the TLSs. If the electron starts 
in a superposition :^(|^)a + \'>tT')a) it will not end up in ^ (|^)^ + but in a 

real mixture of \k)A and |n)^. That is, during the transport, the electron gets entangeled 
with the TLSs. This should not be surprising since we saw a similar behavior already 



in subsection 5.2 without transport (see Eq. (5.27)). However, transport prevents the 
lost information from returning back to the electron periodically. This can be seen in 



Fig |6.3[ b) for transport along three QDs, Xi = Xi and Vtp, 0.$ as in Eq. (6.17) with 
T = 150. The initial states Pj{to) = |l of the TLSs are taken to be a complete mixture. 
As a measure of purity we use TV (p"^)^ . 

On the other hand, as described earlier this section, if the electron starts in |1)^ 
it will get transfered to the pure state \3)a despite being in a real mixture during the 
transfer, as is shown in Fig |6.3| ^a). In this case, ip is an undetectable global phase. 

Now we briefiy discuss the inclusion of measurements once again. In the previous 
subsection the adiabatic theorem was assumed to hold and therefore we first have to 
make sure that it holds also with coupling to TLSs. This is best done by solving the 
Schrodinger equation without measurements as was done to get Fig |6.2[ The solution 
shows whether the transport is adiabatic (almost no fast oscillations) or not (much 
oscillations). If not, the driving fields Qj as well as the transfer time can be increased, 
until the solution shows adiabatic behavior. Then equations (6.9)-(6.13) derived in 
the previous section can be applied here as well, if one uses the appropriate eigenstate 
|'0o(O)i which now depends on the state of the TLSs. However, \ipQ(t)) can only be 
calculated numerically which makes an exact treatment very difficult. Help comes from 
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Figure 6.3: The purity of the reduced density matrix p^{t) during the transfer \1)a |3)a- 
The couphngs to environmental TLS are homogeneous, i.e. Xi = X2 = X3 = X; ^'^d their states 
are in a complete mixture. 

(a) : Transport of a position eigenstate |'0(^o)) = \^)a to IV'(ii)) = |3)a with couplings to envi- 
ronmental TLS X = 0.02 (solid), 0.05 (dashed), 0.15 (dotted). Because of the dependance of 
the adiabatic state |'!/'o(0) on the states of the TLSs, the electron gets entangled with the TLSs 
during transport. However, this entanglement disappears with the completion of CTAP, and the 
electron on the quantum dot rail finishes the transport in the desired pure state. 

(b) : Transport of half of a superposition state. \^lJ{t = 0)) = ^ {\0)a + |1)a) and x — 
0.005 (solid), 0.02 (dashed), O.l (dotted). The transported half of the superposition state picks 
up a dynamical phase which depends on the state of the TLSs. Hence the electron gets entan- 
gled with the TLSs and finishes CTAP in a statistical mixture of |0)a and |3)a, rather than in 
a superposition. 



the fact that the effects of dephasing do not depend much on the exact form of this 



state. Therefore one might prefer to use the unperturbed eigenstate from Eq. (6.3) as 



in section 6.1 to calculate transport fidelities. This leads essentially to two independent 
examinations of transport losses, one for measurements and one for TLSs. This separate 
treatment of the two decoherence effects is justified if they do not influence each other, 
which is generally the case in applications where decoherence is weak. 
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Chapter 7 

Discussion and Conclusion 



We first analyzed decohcrcncc effects on a cliarge qubit without transportation. The 
main result is that while storing information as a charge qubit on a quantum dot rail, 
some of it will leak because of decoherencc due to measurements performed by the 
QPCs. Information loss from decoherence arising from coupling to TLSs, on the con- 
trary, returns periodically. However, as the period is usually not known, the information 
is essentially lost as well. 

The differences between Markovian and non-Markovian noise becomes more striking 
in the electron transport by CTAP. The transport fidelity decrease due to the measure- 
ments is essentially linear in the transport time, therefore the time should be kept as 
short as the adiabatic theorem allows. But quite to the contrary, the reason for less 
fidelity in the presence of TLSs is because the spacing between energy levels decreases. 
This however can be offset by a longer transfer time, as long as the coupling to the 
TLS is not too strong compared to the driving field Qj, and the QDs involved are few 
enough. This is good news for non-Markovian decoherence during adiabatic processes. 

One negative effect of the TLSs is that the transported state will acquire a generally 
unknown dynamic phase during the transport. This limits applications were one wants 
to transport a charge qubit, which might be in a superposition of position eigenstates. 
But as the couplings to the TLSs are unknown, any diabatic transport will face similar 
problems. 

Summarizing we can say, if we can reduce Markovian and non-Markovian dephasing 
sufficiently to be able to store a charge qubit on a rail of QDs, then we can also trans- 
port it using CTAP. Furthermore, adiabatic transport protocols also have the natural 
advantage that they are robust against experimental parameter variations. 

It is often tried to get only a small population in the intermediate QDs, by using weak 
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pulses i^s and Op compared to the coupling Qn between intermediate QDs [H]. This is 
possibly due to the belief, that decoherence is suppressed if only two QDs are populated. 
This is in contrast to our results, because we find that for both, Markovian and non- 
Markovian decoherence, the transport fidelity mainly depends on the energy spacing. 
A large spacing, in turn, is achieved by large couplings between all QD, including the 
pump pulse and Stokes pulse, which results in some population of intermediate QDs. 
The above statement might stem from analyses of STIRAP, where the intermediate 
atomic states are exposed to spontaneous decay, whereas the initial and final states are 
not. This is not the case for CTAP and there is no advantage of not populating the 
intermediate states. 

We should note, that often the use of the spin degree of freedom as qubit is consid- 
ered to be more promising than the position degree of freedom, because of the generally 
larger decoherence time. Then the transport (as well as storage) of position superposi- 
tions is not necessary, therefore avoiding the biggest challenge for CTAP. As we saw in 
subsection |6.2[ the transport fidelity for position eigenstates remains very large despite 
appreciable coupling to TLSs. Therefore, this study provides further evidence, that the 
spin should be the preferred degree of freedom. 



Part II 

Quantum Brownian Motion 
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Chapter 8 

Introduction 



Well over a hundred years ago, Boltzmann derived the linear Boltzmann equation [33] 
to describe a tracer particle affected by molecules of a dilute gas. The special case of a 
heavy tracer particle, or Brownian particle, was studied in more detail by Einstein |44j 
and Smoluchowski jJSj , to describe the phenomenon of Brownian motion, first observed 
by botanist Brown |46j . 

In more recent times, physicists became interested in quantum generalizations, which 
lead to the well established fields of collisional decoherence [21 [371 HH] , as well as the 
Caldeira-Leggett model of quantum Brownian motion [49[ [50] , where the Brownian par- 
ticle experiences a random "force" due to the coupling to a bath of harmonic oscillators. 
Later, quantum versions of the linear Boltzmann equation (QLBE) |51 [ I52 | [53] were de- 
veloped, to describe the interplay between collisional decoherence, friction and diffusion. 
Not only is such a quantum description of collisional Brownian motion desirable, but it 
also sheds light on the non-classical process of decoherence, which is believed to be of 
importance in the quantum classical transition. Furthermore, experimental observations 
of the quantum nature [5l] and decoherence [55] of large molecules improve rapidly, and 
are reviewed in |56j . 

In the study of collisional Brownian motion (classical and quantum), one usually 
assumes that each collision is independent of earlier collisions, which requires that the 
collision time is small compared to the collision frequency. This requirement is typically 
satisfied in the high temperature (small collision time) and low density (low collision fre- 
quency) limit. Also some coarse grained time scale (large compared to collision time) is 
used, such that collisions appear to be instantaneous and the dynamics become Marko- 
vian. 

To fully appreciate the quantum discussion in this field, we state some well known 
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results of the classical theory [57J. On short time scales, the momentum of a particle 
initially at rest diffuses according to Ap^ oc t, whereas the position variance behaves like 
Ax^ oc t^. The reason is that an almost instantaneous collision changes the momentum, 
but not the position of a particle. Only on a much larger time scale, when the momentum 
distribution is close to thermal equilibrium, the position diffuses according to Ax^ oc t. 
Although for large classical particles one is more interested on time scales on which many 
collisions occur, in the study of decoherence of a quantum particle the short time scales 
are of importance, because no quantum behavior is observable after a large number of 
collisions. 

Interestingly, microscopic derivations of QLBEs (which are all Markovian, i.e. do not 
have memory) predict an additional quantum contribution to position diffusion (QPD), 
acting already on the short time scale of momentum diffusion. Although such a process 
is impossible in classical dynamics as it derives from finite position jumps, the very same 
process is needed for a QLBE to be of Lindblad form, necessary for completely positive 
Markovian dynamics. Therefore, QPD is currently believed to arise as a quantum effect 
accompanying collisional quantum friction |5H l58l [59] . In this work, one of the main 
aims is to show the contrary: We find that QPD is not a real physical process, but 
results from using a coarse grained time scale on which collisions appear instantaneous. 
Unproblematic in classical dynamics where collision times of hard-core particles are 
indeed short, the assumption of instantaneous collisions has to be used with care in 
quantum dynamics, where collision times depend on the widths of the colliding wave 
packets. 

Although position diffusion is generally needed to get a positive and Markovian 
master equation, one has to keep in mind that any quantum Brownian motion (QBM) 
Markovian master equation can only be valid in the high temperature and low density 
limit, and on a coarse grained time scale. It can hence be expected, that every derivation 
of a QBM master equation in Lindblad form will find position diffusion in some form, 
even if it is not an actual physical process. Indeed, in the limit of high temperature and 
low density, the predicted QPD vanishes in all proposed master equations. 

To clarify our point, it is convenient to draw a quick comparison to classical BM. If 
no coarse grained time scale is used at all, the collision of a gas particle and the Brownian 
particle will not be instantaneous, but instead, the momentum transfer will take some 
finite time, which depends on the steepness of the interaction potential. However, to 
a very good approximation one can use a hard core potential which is infinitely steep, 
resulting in instantaneous collisions. The effect are random momentum jumps and 
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momentum diffusion of the Brownian particle. In quantum mechanics on the other hand, 
even a hard core potential does not justify instantaneous collisions, as the colliding wave 
packets have a nonzero width. Therefore, there is a finite collision time depending on 
the width of the wave packets divided by their relative velocity, and every Markovian 
master equation can only be valid on a coarse grained time scale large compared to the 
collision time. 

Important contributions to QLBEs include E21 E3 EHl EI]; and are reviewed 
in [58] . Mostly, the analysis is based on using scattering theory to describe the effects of 
a single collision with a gas particle, which is assumed to be in an momentum eigenstate. 
Unfortunately, such collisions with momentum eigenstates impose a major unphysical 
feature regarding decoherence. For simplicity, we use one spatial dimension to point 
out our concern, but the three dimensional case is along the same lines. Assume that 
before the collision, the tracer particle is in a superposition of two momentum eigenstates 
{\p)+ \p'))/V2, whereas the colliding gas particle is in the state \pg). A collision will lead 
to 

{\P{p,Pg))'^ \Pg{p,Pg))+e"'' \p{p' ,Pg))<S> \pg{p' ,Pg)fj , 

where (j) depends on the interaction potential, and the momenta after the collision are 
dictated by momentum and energy conservation 

2mpg + (m — mg)p _ 2mgp — (m — rng)pg 

P ~ i ' Pg ~ \ • 

m + nig m + nig 

Because the gas particle's momentum after the collision depends on the tracer particle's 
momentum before the collision, it follows that the two particles become entangled during 
the collision process. In fact, {pg{p-,Pg)\ Pgip' ,Pg)) = for p 7^ p' , and tracing out the gas 
particle one then finds that all coherences between p and p' are lost in a single collision, 
no matter how small \p — p'\. Furthermore, it follows from linearity that for a general 
incoming state of the tracer particle, the collision leads to a density matrix which is 
perfectly diagonal in the momentum representation. Therefore, a well localized state 
will be spread out infinitely after a single collision. In a three dimensional collision, 
one finds that the tracers particles density matrix will be perfectly diagonal only in the 
momentum component parallel to the momentum transfer, which was also discovered in 
|61| independently from the present author. 

This seemingly unphysical situation can be understood, if the collision time is con- 
sidered. As the collision time of two colliding wave packets depends on their width, it 
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must be assumed to be infinity if the gas particle is in a non-localized momentum eigen- 
state. After an infinite time, it is acceptable that all momentum coherences disappear, 
even if the interaction is only with a single gas particle. 

One might wonder how |51l \52\ [531 EO] were able to derive a QLBE to describe the 
tracer particle at finite times, without encountering the problem of infinite momen- 
tum decoherence in a single collision. For this purpose, we briefly review Hornberger's 
contribution [53], because it seems to be the most complete one. Due to momentum 
conservation, a two particle collision reduces to a one particle problem, where the tran- 
sition operator T with momentum matrix elements related to the scattering amplitude 



(p/|T|pj) = (5(p^ — p^)/(p/, pj)/(7r?i) is of importance. In particular, an expression of 



the form 

X = ^^^(P/ + Ps|f |p, + ps){pi - Ps\f^\pf - Ps), (8.1) 

appears behind an integral over the momentum transfer (p/ — Pi). Here y — t- cxd is 
a box-normalization volume which is taken to infinity. The diagonals of the tracer 
particle's density operator after the collision correspond to ps = 0, whereas Ps 7^ 
describes coherences. For ps = this term contains an ill-defined square of the Dirac 
delta function and a physically motivated replacement rule 



{2TTh) 

V 



3 



Kp;irip.)|2 ^ H^-^^] 'Tj;r"i' • (8-2) 




is used. One should mention that for any representation of the delta function in terms 
of a series of functions, the square of the Dirac delta function goes to infinity. Hence, 
for Ps = 0, X contains an infinite term divided by an infinite volume V, justifying the 
use of a replacement rule to assign a finite value to the expression. On the contrary, for 



Ps 7^ the product of delta functions in Eq. (8.1 ) is well defined. In fact, using standard 
relations for Dirac delta functions, 

X = y;^S{pi-pf)6{ps,\^^) f{pf +Ps,Pi+Ps)f*{Pf -Ps,Pi-Ps) (8.3) 

is found, where q = pj — pj is the momentum transfer, qyp^ is q projected onto p^, 
and q = |q|. The product of delta functions encountered here is well defined, and upon 
integration will lead to terms of finite value. Therefore, X is now a term of finite value 
divided by an infinite volume V, which can only result in zero. This confirms our earlier 
argument, that momentum coherences vanish completely upon a single collision with a 
gas particle in a momentum eigenstate. 
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Diosi |i51J and Vacchini [52j did not encounter this decoherence problem because 
they assumed "quasi diagonahty" of the density operator in momentum basis and only 
studied Ps = 0. Although Hornberger did consider q 7^ 0, he substituted the square root 



of the replacement rule Eq. (8.2) for (pj ibps|T|pj itps) in Eq. (8.1). This assigns a non- 



zero value to X, effectively bringing back the otherwise lost coherences of momentum 
states. The use of the replacement rule for a term which would otherwise vanish is not 
easily justified. Whether decoherence in momentum bases, and the closely related QPD 
are correctly described by approaches of this kind, is surely not certain. 

Also physically, the use of non-localized momentum eigenstates for the gas particle 
seems problematic, because of the resulting infinite collision time. First, any Markovian 
QLBE relies on the assumption of short collision times. Second, successive collisions are 
not independent of each other if the collision time is long compared to the inverse collision 
frequency, leading to the necessity of studying multi particle collisions. Considering that 
diffusion processes usually depend on the considered time scales, an approach where a 
collision time diverges might not be appropriate to discuss QPD. 

To resolve this matter of momentum decoherence as well as of QPD, we study a 
single collision in terms of localized gas states, which avoids the problems of a squared 
delta function. Then a collision time can be precisely defined, and the low-density and 
high-temperature limit will be quantified by ^collision time^ x ^collision rate^ ^ 1. To 
examine a collision event time resolved, and to investigate under which conditions a 
complete collision of two wave packets occurs, we analytically solve the two particle 
Schrodinger equation. This is in contrast to scattering calculations, which postulate 
complete collisions of wave packets. Although solving Schrodinger's equation requires 
more effort, the result will lead us to conclude that quantum collisions can not contribute 
to QPD. 

As our main aim is to clarify fundamental issues of QBM, and not specific details 
involving the interaction potential between the particles, we use the simplest possible 
model: QBM in one dimension and a delta type interaction between Brownian particle 
and gas particle. Furthermore, the gas particles do not interact with each other and 
are assumed to be in a thermal state using Boltzmann statistics. These are both valid 
approximations in the high temperature and low density limit. 

Our approach is based on the theory of generalized measurements, and was intro- 
duced by Barnett and Cresser [60] in a heuristic way. Although the observation that 
a colliding gas particle carries away information about the position and momentum of 
the Brownian particle is certainly correct, they had to guess the measurement operators 
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and needed some unspecified parameters. Here we derive a measurement approach to 
QBM from first principles. In addition to rigorously justifying the approach, this also 
determines various parameters and corrects some flaws of the earlier treatment, such as 
the collision rate. 

This work is outlined as follows. We start by deriving the most important results of 
classical Brownian motion in chapter [9j and introducing some mathematical concepts 



used in this part of the thesis in chapter 10 This should prepare us to draw sensible 
conclusion when studying the quantum counterpart in the following chapters. In chapter 



11 the existing literature is reviewed. Our examination of QBM starts in chapter [121 with 



the study of a single collision of two particles with Gaussian wave packets. Here we find, 
that, if the widths Wg and W of the wave packets of the colliding particles are related 
according to nigWg = mW'^, where rUg and m are the masses of the respective particles, 
then the two particles are not entangled after the collision. This presents us with a 
preferred basis, which will significantly simplify the analysis of QBM. Furthermore, the 
centers of the wave packets in position as well as in momentum representation behave 
as in a classical collision. 



In chapter 13, we work out the measurement, which the gas particle performs on the 
Brownian one. For this we imagine a simultaneous position-momentum measurement on 
the gas particle, which serves as an indirect measurement of the position and momentum 
of the Brownian particle. We use the theory of generalized measurements IQ [8] to derive 
the appropriate effect and Kraus operators, describing the measurement in the Brownian 
particle's Hilbert space. We can then use the Kraus operators to describe the collision 
induced change of a general Brownian particle density operator. 



In chapter 14 we study the statistics of the collisions, assuming a gas in a thermal 
equilibrium. To this end, we decompose the thermal gas density operator in Gaussian 
wave packets. It turns out, that there is a considerable freedom in choosing the width of 
these wave packets, which we will do as to satisfy the conditions introduced in chapter[T2j 
Then we use a projection type operator to regard only the wave packets which collide 
with the Brownian particle in some time interval 5. Putting things together, we finally 



arrive at our master equation in chapter 15 We will explain why this is a valid approach 
to QBM, even though the measurements on the gas particle are not actually carried out. 
Furthermore, we will review all approximations used in the previous chapters and discuss 
under which circumstances they can be fulfilled simultaneously, which will naturally lead 
us to the high temperature and low density limit. 



In chapter 16, we use the master equation to derive equations of motion for the two 
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lowest order moments of position and momentum. We will find classical behavior for all 
but the square displacement in position, where we recover QPD as in earlier approaches 
to QBM. However, we are able to track down the origin of this unphysical feature and 
find its connection to a coarse grained time scale and instantaneous collisions. We 
will establish that position jumps do not really happen, but instead are the result of 
approximations. These approximations in turn are necessary to derive a memory free 
QBM master equation. 



Chapter 17 is devoted to the important study of collisional decoherence, which is 
carried out using the Wigner function of the Brownian particle. We will see that de- 
coherence is caused by phase averaging, as a collision will change the relative phase 
of superposition states. This phase change depends on the state of the colliding gas 
particle, and is therefore random. The often cited information transfer between the 
Brownian particle and the colliding gas particles, which is due to the measurements 
performed by the colliding gas particles, turns out to be small in the range of validity of 
our work. Furthermore, the decoherence rate depends on the spatial separation of the 
two wave functions making up the superposition state, and does not directly depend on 
the momentum separation. These results will lead us to the conjecture, that classical 
Fokker-Planck type equations might be applicable to the Wigner function of a quantum 
particle in a dilute gas. 
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Chapter 9 

Classical Brownian Motion 



We devote this chapter to the classical counterpart to our work. This discussion will 
later enable us to compare quantum Brownian motion with classical Brownian motion, 
and to debate whether the effects discovered in quantum Brownian motion are of true 
quantum nature. 

All of the calculations and results below can most likely be found in the literature. 
However, it turned out to be extremely difficult to find this sort of work. First, most 
work in standard textbooks is in three spatial dimensions. Second, usually the Brownian 
particle is assumed to be much heavier and much slower than the gas particles. 

To draw a comparison between quantum and classical dynamics which is as complete 
as possible, we also need classical results for our situation, i.e. one spatial dimension and 
non vanishing mass ratio between gas and Brownian particle. Although such work was 
most certainly done, possible around a century ago, it is apparently not readily available 
in current standard textbooks. Consequently, we will present a detailed, classical deriva- 
tion of the dynamics of a classical Brownian particle, based on statistical mechanical 
concepts introduced in [62] . 

9.1 Classical collision 

The study of a single collision in one dimension is a very simple one if a hard core 
interaction potential between the two colliding particles is assumed. Then a collision is 
instantaneous and its effect is an exchange of momentum of the two colliding particles. 
If Pg and p are the respective momenta of the gas and the Brownian particle before the 
collision, and pg and p the momenta after the collision, then momentum conservation 

Pg+P = Pg+P (9.1) 
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and energy conservation 



+f = (9.2) 



2mg 2m rua m 



lead to 



^ 2mgp - (m - mg)pg 
m + nig 
- ^ '^rnpg + (m - mg)p 
m + rUg 



Therefore the Brownian particle will gain the momentum q = p — p 

2mpg — 2mgp 
^ 171 + nig 

which depends on both, the Brownian particle's and the gas particle's momentum 



(9.5) 



9.2 Collision probability 

The probability of finding a gas particle with momentum between pg and pg + dpg and 
position between Xg and Xg + dxg is 

Pg{xg^Pg)^^g^Pg = UgHriPg) dxgdpg, (9.6) 
where Ug is the gas particle density, and 

MPg) = ^ , e-^g/^^"'^^^^) (9.7) 

InnigkB-l 

is the Maxwell-Boltzmann momentum distribution for a thermal ideal gas at temperature 
T in one dimension. 

During a time interval dt, a gas particle with a momentum pg will move a distance 
(vg — v) dt relative to the Brownian particle. Therefore, it will collide with the Brownian 
particle, if its position is within the interval (x — {vg — v)dt, x). The probability of the 
Brownian particle colliding with a gas particle with momentum between pg and pg + dpg 
is then 



Pp{Pg)dPgdt = / dXg Pg{Xg,Pg)dpg 

■J x—{vg—v)dt 

= ngiiT{Pg)\vg - v\dtdpg 



p] 



— = exp - 

^2'KnigkBT \ 2nigkBT 



P^_ P_ 
nin m 



dpgdt, (9.8) 



9.3. MASTER EQUATION 



67 



which, of course, depends on the momentum of the Brownian particle itself. 

A number of articles |5H \6U[ [631 El regarding QBM use the Maxwell-Boltzmann 
distribution for the probability distribution for a collision with a gas particle with mo- 
mentum pg. This happens to give the correct thermal stat^ but fails to give the correct 
friction and diffusion constants. Therefore we want to emphasize, that for classical as 



well as quantum Brownian motion, one strictly has to use Eq. (9.8) for the collisional 
probability distribution, if one aims for quantitative results. 



To get the probability for the Brownian particle to experience a momentum "kick" 
q, we substitute pg 
find 



m+rrig 
2m 



q+ -^p, which is according to Eq. (9.5), into Eq. (9.8), and 



Pp{q) dgdt 



(m + nig) 



y/2iTmgkBT im'^mg 



\q\ exp< 



[q{m + nig) + 2nigpy 
Sm?niqkBT 



dqdt. (9.9) 



That is, a momentum transfer with g — )• is impossible. The reason is, that zero mo- 
mentum transfer means that the colliding particles have the same velocity, therefore not 
being able to approach each other. 



Dotted: The exponential in Eq. (9.9) is 



shifted in the opposite direction of the 
Brownian particle momentum, because ac- 



cording to Eq. (9.4) it looses some of its 



momentum during a collision. 



Dashed: Because of \q\ in Eq. (9.9), more 



gas particles with momentum directed in 
the opposite direction of p collide with the 
gas particle. 

Solid: Both effects contribute to the mo- 



mentum transfer Eq. (9.9) being in average 



directed opposite to the Brownian particle 
momentum, therefore resulting in friction. 




m + m, 



9.3 Master equation 



We denote the phase space probability density for the Brownian particle by p{x,p,t). 
The aim of this section is do derive a differential equation, which is called a master 
equation, according to which p{x,p,t) evolves. If the particle evolves freely, then its 



^The reason is, that using Eq. (|9.7[) instead of Eq. (9.8 1 is equivalent to dividing the scattering cross 



section by the relative velocity of the colliding particles, which, of course, also results in the correct 
thermal state. 
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position moves according to its velocity, resulting in p{x,p,t + dt) = p{x — p/mdt,p,t). 
If, on the other hand, the particle experiences a momentum kick q, then we find p{x, p, t+ 
dt) = p{x,p — q,t). Strictly speaking, the position will also change. But because the 
probability of a momentum kick is of the order dt, we do not have to take into account 
the infinitesimal position shift. Adding both possibilities together, we find 

^- j dqPp{q)dt + j dq p{x,p - q,t)Pp-q{q)dt, 

(9.10) 



p{x,p,t + dt) = p[x — —dt,p,t 



m 



where [l — / dq Pp{q)dt\ is the probability that no collision occurred, and Pp^q{q)dt is 
the probability of a momentum kick q, conditioned on the Brownian particle momentum 
p — q at time t. 

Next we substitute 

pL.^dt,p,t) = p{x,p,t)-^dt^P^^^ (9.11) 
\ m / m ox 



into the right hand side of Eq. (9.10), to arrive at the master equation 
dp{x,p,t) p dp{x,p,t) 



dt m dx 



p{x,p,t) j dqPp{q) + j dqp{x,p- q,t)Pp-g{q). 

(9.12) 



The first term represents the free evolution, the second one corresponds to momentum 
jumps away from p, and the third to momentum jumps towards p. 

This master equation is correct for finite mass ratios nig/m, and not limited to a slow 
Brownian particle. The only approximations used are the hard core interaction potential, 
and a low density and high temperature gas such that we can assume independent 
collisions. In particular, we did not use a time coarse graining approximation. 

9.4 Expectation values 



We can use the master equation Eq. (9.12) to derive equations of motion for the ex- 



pectation values of the lowest order moments of position and momentum. Using the 
definition of an expectation value 

{fix,?)) = // dxdpf{x,p)p{x,p,t) (9.13) 



we find for the position 

— (x) = // dxdpxi- — h / dq[p{x,p - q,t)Pp-q{q) - p{x,p,t)Pp{q)]y 

(9.14) 



9.5. FOKKER-PLANCK EQUATION 
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In the second term we change the order of the p and q integrations, and substitute 
u = p — q. It is then seen that the second and third terms cancel each other. The first 
term is integrated by parts and we find 



ip) 



m 



(9.15) 



While integrating by parts, we assumed liiiix^±oo[xp{x,p, t)] = 0, which is certainly true 
for any normalized probability distribution. 
For the momentum we find 



— {p) = // dxdppi- — h / dq[p{x,p - q,t)Pp_q{q) - p{x,p,t)Pp{q)]\. 

(9.16) 



This time the first part vanishes because of limx^±oc [pix, p, t)] = 0, but the second term 
does not cancel the third upon substituting u = p — q, and we find 

^{p) = jj dxdp pix,p,t) j dqqPp{q) 

'JdqqPp{q)Y (9.17) 



In much the same way we derive 



^{x') = 2M, (9.18) 
dt^ ' m ' ^ ' 



(p2> = l^j dq{q^ + 2pq)Pp{q)Y (9.19) 
+ (x j dqqPp{q)Y (9.20) 



d_ 
d± 



dt m 

The integrations over q can not be carried out analytically in general, unless one assumes 
the limit that the Brownian particle is either much slower or much faster than a thermal 



gas particle. We will discuss both limits in chapter 16 

9.5 Fokker- Planck equation 

In order to derive a Fokker-Planck equation, we have to assume that the collisional 
momentum transfer q is small compared to the momentum uncertainty of /9(x,p, t)j^ 

^That is typically satisfied if the gas particles are very light compared to the Brownian particle. 
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Then we can approximate 



p{x,p- q,t)Pp^g{q) ^ p{x,p,t)Pp{q) - q-^[p{x,p,t)Pp{q)\ + — -^[p{x,p,t)Pp{q)], 

(9.21) 



which leads us to the Fokker-Planck equation 
dp{x,p,t) 



dt 



p dp{x,p, t) d 
m dx dp 



p{x,p,t) / dqqPpiq) 



+ 



f g2 
p{x,p,t) / dq—Pp{q) 

(9.22) 



The first term on the right hand side shifts the position according to the Brownian 
particle velocity. The second term is responsible for momentum damping, and the third 
one for momentum diffusion. 

The integrations over the momentum transfer q can not be carried out analytically 
without further approximations. Therefore, we now assume that the Brownian particle 
is slow compared to the thermal velocity of the gas particles, i.e. p{x,p,t) is only non- 



vanishing for p/m <^ ^JkBT jrag. We substitute Eq. (9.5) for q and then evaluate the 
integrals according to the appendix of this chapter. We then find for the friction constant 
7, and momentum diffusion constant mkBTj 

4ngy/2mgkBT p 



IP 



dqqPpiq) 



m 



/q^ y/'^nigkBT 
^QTrPpiQ) ~ ^ngkBT j= , 
2 y/TT 



(9.23) 
(9.24) 



respectively. Substitution into Eq. (9.22) leads us to the Fokker-Planck equation for a 
slow and heavy particle 

dp{x,p,t) p dp{x,p,t) d[pp{x,p,t)] rr9'^p{x,p,t) 

fu = a hjmkBT —2 , (9.25) 

01 m ox op Op'^ 

which is often referred to as Kramer's equation. It satisfies the equipartition theorem. 

9.6 What is a diffusion process? 



It is interesting to note that while we expanded Eq. (9.11 ) only to first order, we needed 



a second order expansion in Eq. (9.21 ) to derive a sensible Fokker-Planck equation. The 



reason is that the first expansion is exact because dt is an infinitesimal, whereas q is finite 
and therefore the second expansion is an approximation. In fact, if we had expanded 



Eq. (9.21) only to first order, then the Focker-Planck equation would not incorporate 
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momentum diffusion. This in turn would lead to a steady state phase space distribution 
p{x,p, t — )• oo) with exact momentum p = 0, and therefore violate the very condition for 
the justification of the expansion Eq. (9.21) (see beginning of this section) 

Although the friction term is due to random momentum jumps of finite size, which 
are in average directed opposite to the momentum of the particle, one would derive a sim- 
ilar term from a continuously applied force which depends on the particles momentum. 
But then the momentum change in an infinitesimal time would also be infinitesimal, 
and the expansion of p{x,p + dp,t) to second order in p would be an infinitesimal of 
second order and therefore vanish. For this reason, for a momentum diffusion process 
to be possible, a single collision has to change the momentum by a finite valuej^ 

Because an instantaneous collision does not change the position of the Brownian 
particle, the Fokker-Planck equation can not exhibit position diffusion. In fact, solving 
the equation for a Brownian particle with initial position x' and initial momentum p' , 
i.e. p{x,p,0) = S{x — x')5{p — p'), we find for short times t <^ 7^"*^, that the position 
variance increases according to Ax^ oc t^, while the momentum variance diffuses with 
oc t. 

The well known position diffusion Ax^ oc t of a Brownian particle is only valid on 
a time scale large compared to the momentum damping time 7~^, that is, once the 
momentum distribution is a thermal one. On this large time scale, the first term on 
the right hand side of the Fokker-Planck equation is responsible for position diffusion, 
because of the randomness of the particle momentum. 

During the thesis, we use the terminology generally used in literature about quantum 
linear Boltzmann equations (see [51] and [58]). That is, a diffusion process has to act on 
a short time scale of single collisions. Therefore, classical dynamics include momentum 
diffusion, but position diffusion is an impossible process as it would correspond to finite 
position jumps of the Brownian particle. A further reasoning for this terminology is that 
in quantum dynamics, the short time scale is of interest because after a large number 
of collisions, any quantum behavior is lost. 



■^If the momentum of the Brownian particle is much larger than mkBT, then the diffusion of the 
momentum distribution will be small compared to the damping, and important physical results about 
friction can indeed be obtained without considering the second order of the expansion in Eq. (9.21 1. 

*Of course, this value can be small and even the limit to zero can be taken after the expansion is 
performed, but that has to be accompanied with an infinite collision rate. 
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Appendix 

We will need several integrals of the form 



/ 



Pg 



pol{pg)exp\ 



ap 



dpg - 



dpg 



ap 



P^ 

mo 



2mgkBT 

P 



m 



pol{pg) exp 



-pI 



2mqkBT 



(9.26) 



where pol(pg) is a low order polynomial. The integral operators in the left square 
brackets act on the term to its right in the obvious way. We now split the term in the 
right square brackets in a symmetric term S and an asymmetric one A, and find 

-ap 



/ap rOG 
dpg - / dpg 
-CO Jap 

/ap poo 
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(9.27) 
(9.28) 



The integration of the asymmetric part can be carried out analytically, in particular 
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(9.32) 



4m^ 



Here, we defined a = mg/m. The symmetric part can not be integrated analytically. 
Therefore, we first expand the exponential and then integrate to find 
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(9.33) 
(9.34) 



Chapter 10 

Tools and Definitions 



We have already discussed general aspects of quantum mechanics in chapter [2| but 
there are more specific results associated with simultaneous position-momentum mea- 
surements, as well as quasi-probability distributions, that find frequent use in later 
chapters. Here, we are concerned with certain mathematical results, grouped into one 
self-contained chapter, from the point-of-view of their use as tools in later calculations, 
rather than providing a detailed discussion of their properties. 



10.1 Simultaneous position- momentum measurements 



We will see in the following chapters that in some sense, a gas particle which collides with 
the Brownian particle performs a simultaneous measurement of position and momentum 
of the latter particle. If the position and momentum uncertainties of the measuremenlj^ 
is at least as large as given by the Heisenberg uncertainty limit ApAx > h/2, then these 
are perfectly valid quantum measurements which can be described in the language of 
POVM which was introduced in section 2.4^ To be consistent with following chapters, 
we denote measurement outcomes with tildes. 



^That are the uncertainties related to the measurement, which are on top of the intrinsic uncertainties 
of the measured state. 

■^Although Heisenberg said that it is not possible to measure both, momentum and position with 
arbitrary precision, what he really had in mind was that the momentum and position of a state can 
only be simultaneously well defined up to some uncertainty. In his time, the simultaneous measurement 
of non-commuting observables was not yet formally developed, and therefore the difference of both 
statement was not clear. Today, some people feel that uncertainty limit related to the measurement 
itself should not be called Heisenberg limit. 



73 



74 



CHAPTER 10. TOOLS AND DEFINITIONS 



10.1.1 Quantum-limited measurements 

We start with the special case of a quantum-Hmited measurement, in which the Heisen- 
berg uncertainty principle is satisfied with an equal sign. With such a measurement we 
associate effect operators which are proportional to projectors onto Gaussian minimum 
uncertainty states. If the system to be measured is in a state p, then the probability 
that a measurement of position and momentum yields results in the respective intervals 
(x, X + dx) and {p,p + dp) is 



where 



and 



Po{x,p)dxdp = Tr[p7ro(x,p)] dx dp, (10-1) 



TTo{x,p) = ^ |x,p)vy(x,p| , (10.2) 



dx'^^=e*"'p/V("-"')'/2W'' \x'). (10.3) 

is a minimum uncertainty state with mean position x, mean momentum p, and posi- 
tion variance W . The uncertainties introduced by these measurements on top of the 
uncertainties Ax(/5) and Ap(/9) of the state itself are 

W 

Ax = (10.4) 

Ap = (10.5) 

Therefore, the parameter W specifies whether the measurement is more precise in posi- 
tion (small W) or in momentum (large W). 

10.1.2 Imperfect measurements 

What we will need in the following chapters are imperfect measurements, which means 
the uncertainties in position and momentum are greater than required by the Heisen- 
berg uncertainty principle. These can be obtained from the quantum-limited ones by 
convolving the ideal distribution Pq{x,p) with a weighting factor w{x,p) 

P{x,p) = jj dxdpw{x,p)PQ{x + x,p + p). (10.6) 



By using Eq. (10.1) it follows that the effect operators for such a "smeared-out" mea- 
surement are 

7r(x,p) = ^— ^ jj dx dpi(;(x,p) |x + x,p + pXx + x,p + p| . (lO-^) 



10.2. PHASE SPACE DISTRIBUTION FUNCTIONS 



75 



In particular, we will be interested in a Gaussian convolution factor 



w{x,p) 



1 



exp 



1 f X 



+ 



(10.8) 



where n specifies the uncertainty due to imperfectness of the measurement. Now we 
find 



Ax 



1 



W\ n + 



h _ 1 



(10.9) 
(10.10) 



For n — )• the convolution factor becomes a Dirac delta function and we regain the 
quantum-limited case. 

We can relate the effect operators according to 



tt{x,p) = D{x,p)ab\x,p) 



(10.11) 



where a = ■7r(0,0), and D{x,p) = e*^^^ xp)/h -g unitary displacement operator. We 
will also need the square root of the effect operators 



\/7r(x,p) = D{x,p)V^D\x,p). 



(10.12) 



10.2 Phase space distribution functions 

In classical dynamics, the phase space probability density p{x,p) is of great importance. 
It gives the probability density of finding the particle at position x with momentum 
p. Because the Heisenberg uncertainty principle forbids that a particle has definite 
position and momentum at the same time, such a probability density can not exist 
in quantum mechanics. Nevertheless, there are several possibilities to assign a phase 
space distribution function to any given quantum mechanical density operator, but each 
of them will lack at least one property of a phase space probability density. For this 
reason, they are also called quasiprohahility distributions. 

10.2.1 Wigner function 

The Wigner function Wp{x,p) is defined by |,8. .65] 

Wp{x,p) = \ [ dx' e-'^'P'"' p{x + x',x-x'), (10.13) 
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where p{x,x') = {x\ p\x') are the position matrix elements of the density operator. Its 
marginals give the position and momentum densities 

j dxWp{x,p) = p{p,p), (10.14) 
J dpWp{x,p) = p{x,x), (10.15) 

where p{p,p') = {p\ p \p') are the momentum matrix elements of the density operator. In 
this sense, it is the closest possible to a phase space probability density, but it is lacking 
the property of positivity. It is well known that only the Gaussian wave functions and 
mixtures thereof have positive Wigner functions. 



(a) (b) (c) 




Figure 10.1: (a): The Wigner function of the incoherent mixture l/2{\x,piXx,pi\ + \x,p2Xx,P2\) 
reminds of a classical phase space density. 

(b) : The Wigner function of the coherent superposition l/2{\x,pi) + \x,p2)){{x,pi\ + {x,p2\) 
oscillates from negative to positive values in the phase space region between the two "classical" 
humps. 

(c) : Position and momentum marginals of the same superposition state as in (b). Although the 
momentum distribution is very similar to an incoherent mixture, the position distribution shows 
an interference pattern which is absent in the mixed state. 



The Wigner function can be generalized to operators other than density opera- 
tors. The Wigner function is real for all Hermitean operators p, and it is normalized 
jj dx dpWp{x,p) = 1 for operators p with trace one. Another remarkable property of 
the Wigner functions is the overlap formula 

Tr(/>ip2) = 27rh jj dxdpWp,ix,p)Wp,{x,p), (10.16) 

which asserts that similar density operators also have similar Wigner functions and vice 
versa, and is not true for other popular phase space probability distributions. 
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Negative and/or oscillating regions in the Wigner function can be seen as signatures 



of nonclassical behavior (see figure 10.1). Therefore, decoherence is associated with the 



vanishing of these negative regions. 

A further usefull property of the Wigner function of a free particl^is that it evolves 
according to the classical equation for the phase space probability distribution 

dWp^t){x,p) ^ p dWp!^t){x,p) ^^^^^^ 

dt m dx ^ 

which means that the Wigner function VFp(t)(x,p) "flows" according to its velocity v = 



p/m. An example is shown in figure 10.2 



t = to 



t = U 



t = t, 




Figure 10.2: The Wigner function of the initial coherent superposition l/2{\x,pi) + 
\x,p2)){{x,pi\ +{x,p2\) during free evolution at three different times. Each point in the dis- 
tribution moves according to the velocity v ~ p/m. 



10.2.2 Husimi Q function 



One may obtain a new phase space distribution function by convolving the Wigner 
function with a Gaussian phase space distribution having the marginals as a minimum 
uncertainty wave packet which is located at the origin 



Qp{x,p) 



1 



dx'dp'Wp{x' ,p') exp 



{x-x'y {p-p'yw 



(10.18) 



As it turns out, the Q function is precisely the phase space distribution obtained from 



a quantum-limited position- momentum measurement Eq. (10.1) 

Qp{x,p) = Po{x,p) 

= {x,p\p\x,p). 



(10.19) 



'This property is also true for a particle in an Harmonic potential. 
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The convolution accounts for a smearing of the marginal distributions, establishing the 
positivity which was missing in the Wigner function. 

Similar to the Wigner function, there is a one to one correspondence between state 
operators and Q functions. Therefore, knowing the precise measurement statistics 
Po{x,p) is equivalent to knowing the state p. This property will turn out to be use- 
ful at times, because if we need to show the equivalence of two operators, it is sufficient 
to compare their diagonals with respect to Gaussian wave packets. However, contrary 
to the Wigner function, very similar Q functions can represent very different states. For 
example, it can not distinguish well between superpositions and mixtures, which makes 
the Q function not suited to the study of decoherence. 

10.2.3 Glauber-Sudarshan P function 

The P function Pp{x,p) is a phase space distribution, which diagonalizes a density 
operator in terms of minimum uncertainty wave packets 

dx dp Pp{x,p) \x,pXx,p\ . (10.20) 

Although this formula is valid for general states p, the above formula seems like a 
statistical mixture of Gaussian wave functions. For states which are not such statistical 
mixtures, the P function will become negative, and might include delta functions and 
derivatives thereof. In fact, the Wigner function, which itself can become negative in 
small regions, is obtained from the P function by a smearing process 



Wp{x,p) = 1 1 da;'dp'Pp(x',p') exp 



The thermal state of a free particle is an important example where the P function 



is positive, as we will see in section 14.1 For this reason, we can assume that a thermal 
gas particle is in a Gaussian state with a probability given by Pp{x,p). 

10.3 Useful formulas 

Representations of position and momentum operators 

It will turn out to be useful to write position and momentum operators in a representa- 
tion diagonal in the Gaussian pure states, resembling the Glauber P representation of 
density operators. 

f f dx dp 

X = // -^-^x\x,p\x,p\ (10.22) 
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f + 



for — )• cxD 
for — )• oo 



P 

2 

fip) 
{xJiP)} 



dx dp 
27:h 



p \x,pXx,p\ 



dx dp n 

T 

27rh 



X \x,p){x,p\ 



dx dp 2 



2TTh 

dx dp 
~2Tph 
dx dp 

2'Kh 
dx dp 

2'Kh 



p \x,py^x,p\ 

xp |x,pX^)P| 
f{p) \x,p\x,p\ 
xf{p) \x^p\x,p\ 



(10.23) 
(10.24) 
(10.25) 
(10.26) 
(10.27) 
(10.28) 



The correctness of these formulas is easily shown by taking their expectation values for 
states described by Gaussian wave functions, i.e. 



(x,p|x|x,|?) 



dx'dp' I 



2TTh 

dx'dp' I 



x'(x,p| x',p')(x',p'| x,p) 

{x-x'f W^{p-p'f 



2'Kh 



-X exp 



2PF2 



2/i2 



X. 



(10.29) 



Because Gaussian states form an overcomplete set of states, an operator is uniquely de- 



fined by its diagonal values (see subsection 10.2.2 ), and therefore the operator Eq. ( 10.22 ) 
is indeed the position operator. The other operators are shown along the same lines. 

Commutators of position and momentum operators 

We will also frequently need the commutator of arbitrary analytic functions of x and p 

[f{x),g{p)] = -f2^-^&Hs:)g^''\p), 

^-^ n! 



n=l 



(10.30) 



where /^"'^ denotes the n-th derivative of /. This can be derived from 

min(fe,Z) 



^ nUk — n)l(l — n)\ 

n=l 

which in turn is shown by induction. 



(10.31) 
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Chapter 11 

Literature Review 



In this chapter, we review some of the existing hterature on QBM. We can roughly 
divide the topic into three categories, each of which is reviewed in one section. The 
first assumes, that the Brownian particle in infinitely heavy. This simplification leads 
to the study of collisional decoherence, where any change of the Brownian particle posi- 
tion is neglected. The Caldeira-Leggett model in comparison, incorporates friction and 
diffusion of the Brownian particle, but achieves this using an environment of harmonic 
oscillators rather than gas particles. Furthermore, the interaction potential between 
Brownian and gas particle is linearized. Finally, there are quantum versions of the lin- 
ear Boltzmann equation, which were developed to study the full interplay of friction, 
diffusion and decoherence of a Brownian particle in a gaseous environment. This ap- 
proach can be considered the genuine quantum counterpart of classical Brownian motion, 
as studied, e.g. by Einstein |44| and Smoluchowski |45j . Because our own work belongs 
to the latter category, we will review publication of this type in more detail. We close 
this chapter with a brief section, to draw some critical conclusion which are of relevance 
to our work. 

There are yet other approaches to QBM. Some do not specify the Brownian environ- 
ment of the particle, but rather postulates some properties of QBM. These properties, 
like spatial invariance or momentum damping, are related to classical Brownian mo- 
tion. The master equation is then constructed such that the dynamics of the Brownian 
particle fulfills these properties. The articles of Gallis [66] and Holevo |67j use such a 
postulation ansatz. 

Our work questions some of the popular results of quantum linear Boltzmann equa- 
tions. Therefore, we will review some of the papers in more detail, to point out where 
difficulties in each approach might occur. However, the literature is too large to account 
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for every contribution in such detail. We therefore restrict a comprehensive review to 
only a view of these, and mention only some of the others rather briefly, referring to the 
well written review article of Vacchini and Hornberger |58] . 

In order to avoid difficulties arising from the variety of notations used in these papers, 
we changed some of the notations of the original articles. Variables and values with index 
g correspond to gas particles whereas no index corresponds to the Brownian particle. 

11.1 Collisional decoherence 

The paper of Joos and Zeh [2\ initiated the quantitative study of environmental de- 
coherence due to scattering of photons and molecules. They discuss this phenomenon 
on various macroscopic systems, like the position of a macroscopic particle and the chi- 
rality of macro molecules. Due to our interest in QBM, we restrict this review on the 
first of these examples, and refer the reader to the very well written original article if 
interested in the non-observability of coherences in our macroscopic world. We should 
mention that the authors were mainly interested in a qualitative understanding of the 
decoherence process, and therefore made use of several qualitative formulas. 

Their discussion about collisional localization of a macroscopic particle is restricted 
to recoil- free collisions, which effectively means that the Brownian particle is infinitely 
heavy. The scattering of a gas particle by the Brownian one is described by 

|r)|V^,) ^ |r)|V<,(r)) = |r)Sr (11.1) 

where Sr is the gas particle scattering matrix for a Brownian particle at position r, and 
\ipg)is the state of the gas particle before the scattering event. Therefore, the Brownian 
particle's density matrix after the scattering event is 

p(r,r') = poir,r'){iJg\Sl,Sr\^Pg), (11.2) 

where po{r, r') is the Brownian particle density matrix prior to a collision. If the scatter- 
ing interaction is invariant under translation in space, 5'r(kg, k^) = 5(kg, ]i'g)e^^^^^^^9^'' 
can be used. They continue with 

ii^glSlSrli^g) = J d3fc,d3A;;d3A:;'S(kg,k;)5*(k3,k;')e''^«(r'-')e*'^«'-<r'c(k;)c*(kp, 

(11.3) 

where c(kg) is the momentum representation of \4^g). Now they approximate the gas 
particle state by a momentum eigenstate c(kg) = L~^/'^6^(kg — Pg), where is the 
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momentum of the incident particle and L^ is the box normahzation volume. Next the 
scattering matrix is expressed in terms of the scattering amplitude 



f{kg,k')S{kg - k'). 



:il.4) 



Inserting this relation into Eq. (11.3) leads to a problem shared by all approaches which 
use momentum eigenstates of a gas particle to study a collision: The occurrence of an 
ill defined squared Dirac delta function. Joos and Zeh therefore simply substitut^ 
[<53(kg - k^)] /L3 ^ 6^{\ig - k!g), as well as [6{kg - k'g)f /L 6{kg - k'g). 

Expanding the exponential in Eq. ( 11.3[ ) to second order and performing the integral 
over the angles of kg eventually yields 



P(r,r') 



Po(r,r') 1 



87r2L2 



-CF, 



polr,r)exp' 



-cr, 



where 



eff 



dcosG|/(cose)|2 [(2-cosG)^ 



1 



(11.5) 



:ii.6) 



The effect of N subsequent collisions is given by multiplication of the exponent by A^. 
The number of collisions in a time interval 6 can be given in terms of the flux Ug/vg 



N 



L'^UgVgS, 



[11.7) 



although the appearance of a surface related to the box normalization volume instead 
of the total scattering cross section seems questionable. Taking the limit 5 — )• 0, they 
finally find 

dp{r, r') 



-A(r-r')V(r,r'), (11. 



dt 



with the localization rate 



A = ^-JL^H_^_ (11.9) 



It is important to note that because of the expansion of the exponential in Eq. ( 11.3 ), 



the result Eq. (11.8) is only valid for small distances (r — r') -p^ ^ 0. For large distances 
the localization rate in Eq. (11.8) increases to infinity, which is clearly not physical as 
the localization rate should be bound by the collision rate. 



^They do not mention this procedure, but it is clearly carried out. 
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They continue to calculate A for various special cases by averaging over the mo- 
mentum distribution pg, including not only scattering of gas particles, but also thermal 
photons. Coherence length after t = Is are given. Next Joos and Zeh include the internal 
system Hamiltonian H = /2m. As this part is quite long and does not qualitatively 
change any results, we skip its discussion an refer to the original paper [2]. 



A problem in Eq. (11.8) is that it leads to localization without bound, which due to 
Heisenberg's uncertainty principle corresponds to an increase of mean energy towards 
infinity. This problem is discussed in [68], and occurs because any movement of an 
infinitely heavy Brownian particle is neglected from the beginning. The inclusion of the 
unitary free particle dynamics later in their paper can not avoid this "heating" of the 
Brownian particle. 

Joos and Zeh's results were refined by the article of Gallis and Fleming |47j . 

where they discuss the problem of the transition from the quantum world to the classical 
world. In one part of their work, environmentally induced decoherence is discussed. In 
particular, Gallis and Fleming assume an environment which scatters off the system of 
interest, much as in QBM. 

The explicit calculation is very similar to the one of Joos and Zeh, and not repeated 



here. The main difference is that they do not expand the exponential in Eq. (11.3), and 
therefore their result 
dp{T,r') 

— = -F(r-r)/)(r,r ) 

F{r) = /dp,MP.)^/^^(l-e<P-^^)^)l/(p„p;)P, (11.10) 

where ^xiPg) is the Maxwell-Boltzmann distribution, should be valid for general spatial 
separations of superposition states. Indeed, for large r, the exponential oscillates quickly 
and does not contribute to the decoherence function F{r), guaranteeing that the latter 
is bound for large separations. For small separations r — )■ 0, they obtain the familiar 



result Eq. ([1L8|). 

Hornberger and Sipe 08] finally managed to avoid the ill defined square of the 
Dirac delta function, by applying scattering theory to normalized states. In doing so. 



they corrected Gallis and Fleming's result Eq. (11.10) by a numerical factor. 

Performing a careful scattering calculation, and assuming the limit of an heavy gas 
particle, they find that the scattering operator S transforms the two particle state |r)|i/)g) 
according to 
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= (11.11) 

Here, \tpg)is a general, but normalized gas particle state, Sg is the one particle scattering 

operator acting on the gas particle, and is its momentum operator. Therefore, if the 
two particle state before the scattering event is ^ = po ® iV'sXV'sl with 

p = y"dVdVp(r,r') |rXr'| , (11.12) 

then the scattering of a gas particle state \tpg) results in 

p{r,r') = « <Vo(r,rO. (11.13) 



Now they proceed by introducing the transition operator Sg = 1 + iTg. The scalar 
product is then 



) = l+{lPg\A\i^g), (11.14) 



where the operator A can be expressed in terms of the transition operator Tg. To 
calculate the change of the Brownian particles density matrix due to one collision, they 
insert two complete sets of momentum eigenstates to find 

ii'Mli^g) = jd%dyg{^g\p'g){l>'g\A\Pg){l>g\i^g). (11.15) 

In this equation, the advantage of using momentum eigenstates instead of |'0g)is appar- 
ent, as only diagonals of A would be needed. But A involves terms like fgTg, and if 
a complete set of momentum eigenstates is inserted between these transition operators 
(because the transition operator involves an energy conserving Dirac delta function), 
one would again encounter the square of a Dirac delta function. 

To avoid this problem, they chose Gaussian minimum uncertainty wave packets 
\rg, pg) for the gas particle states. The consequence is a very long calculation, which we 
drop here, but rather present the result 

{^g\A\^g) = / A^^id), (11.16) 

where b is the momentum variance of the state \rg, pg). Furthermore they defined 
A^^id) = rq(rg-R)y"dne'('i-«")-('-"-')/^|/(9n,q)|2 

[rq(r,-r) + rq(r,-rO] J dn|/(gn,q)|2 
+ — [rq(r, - r) + T^ivg - r')] Re[/(q, q)], (11.17) 
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where n is a unit vector, R = (r + r')/2, 

exp{-[i?2_(R.q/g)2] /«2} 



rq(R) 



(11.18) 



and a is the position variance of the state \rg,pg). 

Now they turn to a colhsion with a thermal Boltzmann gas. Most important for the 
apphcation of the previous part of their paper is the convex decomposition of a thermal 
gas particle state into minimum uncertainty wave packets \ipg) = \rg,pg), with mean 
position Tg and mean momentum Vg. This is reviewed in subsection 14.1, and the result 
is 



Pa 



V 



(11.19) 



where V is some normalization volume. The momentum distribution ^^(pg) is similar to 
the Maxwell-Boltzmann distribution /UT(pg) = {2-KmgkBTy^/'^e"^y^'^'^^'''''^\ but with 
the temperature T substituted by a slightly smaller value, as some of the particles energy 
is related to the momentum uncertainty of the states |xg,pg). In the following these 
states are chosen such that this momentum uncertainty is small, and the approximation 
fifiPg) ~ /^r(Pg) can be used. 



B 



A 



Spg/ma 



Figure 11.1: The volume between the dashed lines is R{pg). During a time interval (0,(5), a 
collision with the Brownian particle (blue) can occur, if a gas particle (grey) with momentum 
Pg is in R{pg). 



Now a collision can be pictured in an almost classical way. Assuming that the 
Brownian particle is located within a small area, as are the gas particle states \ipg) = 



11.1. COLLISION AL DECOHERENCE 



87 



\rg,pg), then we can think of a volume R{Pg) depicted in Fig. 11.1 A gas particle will 



collide with the Brownian one during a time interval 6, exactly if it is in this volume. 



Although particles like the particle B in Fig. 11.1 will miss the Brownian particle, 
here it is considered to be a scattering particle. That is acceptable, as the authors 
show that the scattering of a gas particle which "misses" the Brownian particle does not 
contribute much to the decoherence of the Brownian particle's state. 

Therefore, the change of the density matrix during a time intervall 6 is 



Ap(r,r') = Np{r,r') [ d%fiT{Pg) [ ^ {rg,pg\ A\rg,pg), (11.20) 



where A'^ is the number of gas particles in the volume V . After substitution of Eq. ( 11.16 ) 



and some calculation, they arrive at Gallis and Fleming's result Eq. (11.10), with a small 
correction in that the right hand side is multiplied by 2ti. 

Finally, they repeat the calculation of Gallis and Fleming |47j with a modified and 
physically motivated replacement rule for the squared Dirac delta function 

where cr(pg) is the total scattering cross section for an incoming gas particle momenta 
Pg. The recovery of their earlier results shows that the much simpler approach, us- 
ing delocalized momentum eigen states can also be successful. It further gives some 
confidence in the application of the replacement rule, at least for an infinitely heavy 
Brownian particle. 

The work of Hornberber and Sipe is, up to now, the only one found in the literature to 
avoid squares of Dirac delta functions. Their results should therefore correctly describe 
collisional decoherence. However, their calculations are only valid for an infinitely heavy 
Brownian particle. As such an particle does not move in position space, issues around 
friction and position diffusion can not be resolved. 

A different replacement rule is motivated by Adler [69] ■ He uses a particular repre- 
sentation of the Dirac delta function 

5{E) = ^ f d^e'^*/^ (11.22) 



from which he follows 



m = ^, (11.23) 
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where A is an elapsed time interval large compared to the collision time. He continues 
along the lines of Hornberger and Sipe's replacement rule calculation, and obtains their 
results. 

It is however interesting to note, that, by using a different, but equally valid rep- 
resentation 5{E) = 2^ J dt exp{iAEt/h) of the delta function, Adler's argumentation 
leads to a different replacement rule 5{0) = where A is an arbitrary real number. 

Dodd and Halliwell [6l] used non-relativistic many body quantum field theory, to 
essentially derive the result found by Hornberger and Sipe. They also need to rely on 
the replacement rule. At the end of their article, extension to a finite Brownian particle 
mass are considered, but these are even in need of a replacement rule for the derivation 
of a Dirac delta function. 

11.2 Caldeira-Leggett model 

The Caldeira-Leggett model of QBM, which is readily available in textbooks [5j , consists 
of a Brownian particle in a potential V{x), coupled to a bath of harmonic oscillators via 
a linear interaction. The Hamiltonian of this system is 

H = Hs + He + Hi + H, (11.24) 
Hs = ^ + V{x), (11.25) 

i 

Hi = -f J^KiXi, (11.27) 

•i 

He = i'E^^ (11-28) 

i 

where uJi is the frequencies of the i-th harmonic oscillator, rrii is its mass, and Kj is 
its coupling constant to the Brownian particle. The interaction term leads to a renor- 
malization of the Brownian particle potential, and the counter-term is included to 
counteract this potential. 

In the derivation of the master equation, the spectral density 

2m,- a;,- 



J(a;) = y^-^5{u:-Ui) (11.29) 



is of importance. The spectral density can be a smooth function, if the limit of a 
continuous distribution of oscillator frequencies is used. In particular, an Ohmic spectral 
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density 



J(uj) = —^oj (11.30) 
vr 



is often used, where 7 turns out to be the damping rate. In the limit of high temperature, 
it is possible to derive a simple master equation 

^ = (11.31) 

which was first achieved in the seminal paper of Caldeira and Leggett |49j . Although 
they used the complicated Feynman- Vernon theory of the influence functional, there are 
now simpler derivations of these results available [5]. 

Because this master equation is not of Lindblad form, it can not be completely 
positive. However, it can be cast into Lindblad form by adding a position diffusion term 



8mkBT 



which is small for high temperatures. It is interesting to see, how approximations used 
to derive Markovian dynamics can lead to unphysical results, as the lacking positiv- 



ity of Eq. (11.31) shows. To establish positivity, a position diffusion term has to be 
added, even if there is no physical reason to do so. However, as long as any unphysical 
property is small under the conditions of validity of the master equation, the master 
equation can perfectly be used to calculate physical quantities, which can be tested 
against experiments. 

The exact master equation for the Caldeira-Leggett model has also been derived |50j . 
and is known as the Hu-Paz-Zhang master equation. As an exact equation valid at any 
temperatures, it is non-Markovian and does not suffer any problems such as the violation 
of positivity. 

11.3 Quantum linear Boltzmann equation 

The contribution of Diosi [51j was the first one to address the interplay of friction, 
diffusion and decoherence in collisional QBM. He also started by considering a single 
collision in the interaction picture, and uses an initially uncorrelated state g = p® Pg, 
where pg is a thermal state with a Maxwell-Boltzmann momentum distribution pxiPg) 



from Eq. (9.7). He introduces the transition operator by = 1 + iT, where S is the 



unitary scattering operator. He finds for the change of the two particle density matrix 
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due to one collision 



(11.33) 
(11.34) 



For the evolution of the Brownian particle, he approximates dp/dt by 

At At ^ ^' 
where At is a time long compared to the scattering time. 

Now he substitutes the standard form of the transition operator 

T = ^y"d3pd3p,dX/(p^pr)'5(^PS-v)ip-p;;+p.'P>'P.i' 

(11.35) 

where quantities marked by stars are centre of mass quantities (see [51] for precise 
definition). After performing the partial trace and changing the integration variables, 



he finds for the second term in Eq. (|11.33|) 
27rn„ f m + m 



1*2 



|/(P9'P9. 



d^pdV d^P^ dX* /i(p; + ^(p + p^') 

2 



m„ 



|P - Pg + Pg)\P\p\v>'\v' -Vg + Vg\ 

(11.36) 



where Ug is the density of the gas particles. In the derivation of this formula, he assumed 
that the density operator of the Brownian particle is almost diagonal in momentum 
representation. 

In Diosi's approach, there also appears the square of the Dirac delta function. He 
removes this ill-defined term by applying 5{E = 0) = At/27r. By rewriting the previous 
equation in operator form he finds 



TiJfQT^ 



Ug / m + nig 



m 



*2 



m 



y dVg - e^t) I/(PP' pDI' ^PSPr^^pVr 



p; + ^(p + p';))e-^(p«*-p^)* 



m 



(11.37) 



The third term on the right hand side of Eq. (11.33) can be calculated in a similar 
way, whereas the first term can be neglected in rarified gases. Performing an integration 
over the Dirac delta function, he finds the master equation for the Brownian particle to 
be 

3 



dp 
dt 



rin 



m + niQ 



m 



dE*dn*dn'* k*^ Mo*,E*) 

dQ' 



(11.38) 
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where da/dil = |/p is the differential cross-section. This equation is Diosi's central 
result. From the condition dp/dt = 0, he obtains the correct thermal equilibrium state 

p(t^oo) oc exp(-^) (11.39) 

at temperature T, and therefore, Diosi was the first to avoid the unbound heating of 
infinite mass approaches. 

Diosi also derives the familiar interaction picture master equation for QBM, 

^ = -^^[i,{p,p}]-^[i,[i,p]]-^[p,[p,p]], (11-40) 
by using the limit of a heavy Brownian particle, with a small coherent extension. 



Both, Eq. (11.38) and Eq. (11.40) show position diffusion which led Diosi to the con- 
clusion that quantum friction always has to be accompanied by finite, random position 
jumps. We will show in our approach that this is in fact not true. 

Despite rather nice results, there are some shortcomings in Diosi's calculations. First, 
he assumes that the probability of a collision of a gas particle with momentum pg with 
the Brownian particle is independent of its momentum. This can certainly not be 
quantitatively correct, as e.g. two particles with the same velocity can never collide. 
Gallis and Fleming accounted for this difficulty by introducing a particle flux, which 
is proportional to the particle's velocity. For a Brownian particle with finite mass this 
issue becomes less trivial. 

Second, the approximation of the Brownian particle density operator being almost 
diagonal in the momentum basis is not trivial, as in the calculation, it is rather assumed 
that the density operator is in fact diagonal. Although this assumptions seemed appro- 
priate at the time (it was also used by Vacchini [52] ) , we [TO] as well as Diosi [61] himself, 
showed later that a collision with a gas particle in a momentum eigenstate removes all 
coherences parallel to the momentum transfer (see also chapter [s]) . To account for this 
effect, one has to consider off-diagonals of the Brownian particle density operator in 
momentum basis. 

Vacchini |52| uses the formalism of non-relativistic multi particle quantum field 
theory. Similar to Diosi, he also assumes that the Brownian particle's density matrix is 
almost diagonal in momentum basis to derive a master equation. By taking the limit of 



light gas particles, he obtains the familiar form Eq. (11.40). However, his coefficients 7, 
Dpp, and Dxx are quite different to the ones found by Diosi |51j . 

In the derivation of the master equation, there appears no square of an energy 
conserving Dirac delta function, because EinUiai — Ejinai = is avoided by adding te 
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to the energy change. There are no reasons given for this procedure, and its vahdity is 
difficult to judge. 

The letter of Hornberger |53] seems to be the most complete approach to collisional 
QBM to date, in that he is the first to explicitly consider off-diagonals in the momentum 
representation of the density operator. To do so, he introduces an extension to the 
replacement rule, as to be applicable to off-diagonals. 



As Diosi, he starts from Eq. (11.33). The incoherent part can be expressed in terms 
of the kernel 

(p|Tr,(f[|p"Xp">P,]ft)|p'>. (11.41) 



Using the standard representation of the transition matrix, he finds, that the momentum 
transfer q has to be the same for both cohering momenta, i.e. p — p" = p' — p'"- Inserting 
a momentum diagonal representation of the gas particle density operator, Hornberger 
arrives, for p = p', at the familiar squared Dirac delta function. He uses the replacement 
rule 

where is a box normalization volume, Tq is the single particle transition operator, 
/(PfjPi) is the scattering amplitude and cr{pg) is the total scattering cross section. 
Next, Hornberger turns to the the off-diagonals p 7^ p', which were avoided by previous 
studies. He finds terms like 

(P/ + Ps\To\pi + Ps){Pi - PsITqIp/ - Ps), (11.43) 

where = rng{p' — p)/[2(m + mg)] corresponds to the momentum separation of 
the cohering momenta. He then substitutes twice the square root of the replacement 



rule (11.42), which yields the square root of a product of two energy conserving 6- 



-.2 „2 



functions, with arguments ^^2^' =t (pj — Pi) ■ Ps- He argues therefore, that Ps should 
be replaced by its projections Ps,±q onto the plane vertical to the momentum transfer 
q = pj — Pi in all expressiona^l and is left with a single, proper Dirac delta function 
6 



Pf-pt 



2 

To account for the collisional statistics, he introduces a two particle rate operator 
R = y dpdpgngw(p,pg)cj(rel(p,pg)) IpXpI (g) |pg)(pg| , (11.44) 



^To give a one-dimensional analog, instead of the rule J dx f(x)5{x ~ y) ~ f{y), one uses 
J dx f{x)5{x — y) ~ J dx f{y) = f{y) J dx, and cancels the now infinite integral, with an infinite nor- 
malization length. 
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to determine the collision rate with the gas. Here, Ug is the particle density of the gas, 
v{p,Pg) is the relative velocity of gas particle and Brownian particle, and rel(p, Pg) = 



m+m„ 



Pa 



p is their relative momentum. He then returns to Eq. (11.33) and 



(11.41), and substitutes T with Ty R. Eventually, he finds the master equation 

^2 



dp 

dt 



1 



p 

2m' 



+ /dq 



1 



dk , , 



(11.45) 



where -£'q,k = e''^''^^^F{'k, p; q) and 

^(m + rUg) 
F(k,p;q) = ^ — 



/ 



q q 

rel(k_Lq, p±q) - -,rel(k_Lq,p_Lq) + - 



k_Lq + 



nig + m q 



m 



O + ^P||q 

2 m " 



The second integration is over the plane perpendicular to the momentum transfer q, 
and a vector with index ||q denotes the part of the vector parallel to q. We wish to 
note, although it is not pointed out in the article, that this QLBE also displays QPD. 

We showed in the introduction (chapter [s]), that the use of the modified replacement 
rule for momentum coherences is problematic, because the term (11.43) is perfectly 
well defined without any replacement rule. However, if Hornberger had not used the 
replacement rule, he would have found that a collision with a momentum eigenstate 
destroys all momentum coherences, which is also explained in the introduction. These 
observations indicate, that one might have to use normalized gas particle states, to 
describe decoherence aspects of QBM. 

There is a further questionable element in the use of the modified replacement rule. 
To outline this objection, we write for the density matrix before a collision 



p(p",p" 



2(m + mg)ps/mg) 



(11.46) 



A strict application of Hornberger's argumentation, that the momentum separation p^ 
should be replaced by Ps_Lq in arguments of all functions, would require it to be used 
also in the argument of the density matrix. However, replacing the off-diagonals of the 
density operator by entries which are on the diagonal in the direction q, clearly does 
not make sense, and was not done in the article. It seems unsatisfactory, that one has 
to choose when to use a projected vector, such as not to obtain unphysical results. 



On the other hand, applying the delta functions Eq. (8.3) without a replacement 



rule, one finds that a collision with momentum transfer q leads to 

p(p,p') oc 6{ps-(i)p{p" -ci,p" -ci + 2{m + mg)ps/mg) 



(11.47) 
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which should be vahd for 7^ 0. That is, a mathematically strict treatment using 
the delta functions shows, that only coherences with momentum separation orthogonal 
to the momentum transfer survive a collision, which is exactly what we derived in the 
introduction by a different method. 

Despite this criticism, we would like to emphasize, that Hornberger was the first to 
address the very complicated problem posed by the off-diagonals of the density opera- 
tor. As such, his valuable article is the starting point for a number of further investiga- 
tions [7lJ, and heavily influenced the work presented in this thesis. A somewhat more 
detailed derivation of his master equation, including an alternative, although not much 
more stringent motivation of his replacement rule, is found in [72j. 

A somewhat alternative approach to QBM was offered by Barnett and Cresser |60j . 
Motivated by the difficulty of collisional approaches, they formulated a measurement ap- 
proach to QBM, recognizing that a colliding gas particle performs a measurement on the 
Brownian particle. They could then use the machinery of generalized measurements to 
derive their Lindblad operators, and did not encounter the dubious squared Dirac delta 
function. Most reasonable seemed simultaneous position-momentum measurements, and 
as a gas particle can not be considered as a perfect measurement apparatus, imperfect 



measurements of the form described in subsection [10.1.2 where postulated. 

These measurements are specified by two parameters. First, W determines whether 
the measurement is more precise in position or momentum, and second, n specifies the 
imperfectness of the measurements. In fact, n = corresponds to a quantum-limited 
phase space measurement. 

The Kraus operators for the collisional transformation of the Brownian particle den- 
sity operator are then given by 



and include the back action of the measurement, which is the classical momentum 
transfer due to a collision. Furthermore, they included the statistics of the collision 
due to the thermal momentum distribution ^xiPg) of the gas particle. 

Assuming that the measurements (or collisions) are instantaneous, the master equa- 
tion 



dp I 
dt ~ h 



p2 
2m 



Rp + R /// dpgdxdpAp^{x,p)pAl {x,p) (11.49) 



is derived, where R denotes the total collision rate. They applied the limits — )• 00, 
mg —7- 0, and n — t- 00, such that 7 = 2mgR/{m + mg) and R/n = are kept finite, to find 
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a master equation of the well known form Eq. (11.40). Their diffusion coefficients are 



= 7H.T+-^A./) + — , (11.50) 



i^xx = (11.51) 



Rh^ 

where the variances Ao-x^ and A^-p^ correspond to the uncertainties associated with the 
simultaneous measurement of position and momentum, over and above the particle's 
intrinsic variances. 

It is interesting to note, that this approach gives some physical insight to the mys- 
terious position diffusion, as it can be traced back to the state reduction postulate 
of measurements in quantum mechanics. The better the momentum resolution of the 
measurements, the larger the position diffusion. 

There are some inconsistencies in their derivation. As in Diosi's article, the collision 
statistics do not account for the relative velocity of the two colliding particles. Better 
collision statistics would be obtained by replacing uriPg) with \pg/mg — p/m\^T{Pg)- 
Next, the collision rate R is not given in terms of physical properties of the gas. In fact. 



by using correct collision statistics, we will show in chapter 14 that the total collision 
rate is R = ng\/2kBT j ^T^rrig. Furthermore, there are several more parameters, which 
can not be obtained from first principles. 

More important is that the Brownian particle, if described by the master equation 



Eq. (11.49) (or the limiting one for rrig — )• 0), does not approach the correct thermal 
state. The authors claim that the difference compared to the expected thermal state 
is due to quantum corrections. However, by considering the special case of a Brownian 
particle, which is of the same type as the gas particle, i.e. rag = m, it is apparent that 
this can not be true, as then the Brownian momentum distribution must approach the 
very same momentum distribution of the surrounding gas. 

Nevertheless, many of their ideas, especially in regard to the measurement interpre- 
tation of a collision, are of importance throughout our work. 

Breuer and Vacchini [7T] performed a quantum trajectory analysis of Vacchini's [52] 
QLBE. The quantum trajectory method is a popular technique of simulating general 
master equations, and is well described in the book of Breuer and Petruccione [5]. 

They define the scaled Brownian particle momentum 

u = p/{mvmp), (11.52) 
where fmp = \/'^kBT/mg is the most probable velocity of a thermal gas particle, and 
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m is the mass of the Brownian particle. Therefore, u = 1 means that the Brownian 
particle has the same velocity as an average gas particle. Although the results about 
dissipation seem very nice, the physical parameters chosen in the decoherence chapter 
are not of much relevance for QBM. For the most part, the masses of the gas particle 
and the Brownian particles are the same. When they finally take the limit of a light 



Brownian particle, they also assume ti ^ 1, which upon definition Eq. (11.52) means 
that the Brownian particle is fast compared to the gas particles. For a typical gas, even 
at cold temperatures, one would have to prepare a Brownian particle in a superposition 
of velocities itv, where v is of the order of several hundred meter per second. 

In particular, superposition states of the form |uo)+ |— uq) are considered. Although 
one would expect, that a single collision with a particle of the same mass destroys 
all coherences, they find (their figure 10) that coherences vanish only once the scaled 
momentum approaches one. We further wish to point out, that the 'decoherence per 
collision' to lowest order in uq is 2ug/3, with the only temperature dependence being due 
to the scaled momentum. This is obtained by dividing their decoherence rate from their 
Eq. (67) by their collision rate from their Eq. (31). For a given Brownian momentum 



superposition state |po)+ |~Po)i that means that from Eq. (11.52), the 'decoherence per 
collision' decreases as the temperature of the gas increases, which certainly can not be 
true. 

We believe, that the reason for these counter intuitive results could be found in 
their definition of coherence in their Eq. (62), which does not seem to take into account 
dephasing effects. As we will see in chapter [iTj dephasing is the main contribution to 
collisional decoherence. 



11.4 Difficulties with current theories 

There are two main observations from this literature review: First, it seems well ac- 
cepted, that quantum dissipation goes along with position diffusion, and hence finite 
position jumps. This opinion originated even before detailed studies on QLBE emerged 
(see e.g. [73J), and is obtained throughout later articles by different authors, including 
the recent review article [58j by Vacchini and Hornberger. 

The other observation is about the ambiguity of current derivations. Despite the fact 
that QPD is clearly related to momentum decoherence, so far there is no satisfactory 
treatment of the change of momentum coherences due to a collision. Furthermore, 
all derivations of QLBEs assume instantaneous collisions, which are at odds with the 
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diverging collision time associated with colliding gas particles in delocalized momentum 
eigenstates. This is especially problematic, because diffusion processes often depend on 
the considered time scales. 

It is interesting to note, that, the master equation of the Caldeira-Leggett model in 
the Markov limit predicts slightly negative probabilities for certain initial states |74j . 
Yet, it is clearly the case that these negative probabilities are unphysical, i.e. that they 
simply arise as a consequence of approximations used to derive the master equation. 
Similar, the master equation for collisional decoherence predicts slow, but unbound 
heating of the Brownian particle. Also this effect is clearly not of physical nature, but a 
result of mathematical approximations. Nevertheless, to our knowledge, the possibility 
that the same could be true for QPD never occurred in the literature, despite such a 
position diffusion process demands non-continuous evolution of the density probability 
(in the sense described in detail in Chapter [9]), and therefore also for the wave function. 
In this sense, the discovery in the following chapters, that QPD can not arise from 
random collisions, is not really that surprising. 
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Chapter 12 

Single Collision 



The obvious starting point of the study of a particle undergoing random colhsions with 
gas particles, is the discussion of a single collision. Therefore, this chapter is devoted 
to the examination of a one dimensional collision of two point particles. For reasons 
discussed in chapter [8j we assume that both particles are initially described by Gaussian 
wave functions 

„—ixp/ (2h) 

{x'\x,p)w = , (12.1) 

The relative phases of these states are chosen such that = D{x,p) |0, 0) is satisfied, 
where 

i)(x,p) = e*(P^-^P)/^ (12.2) 

is the Glauber displacement operator. The index W indicates the position variance, and 
will be dropped at times for shorter notation. 

We assume a hardcore interaction potential such that the two particles cannot tunnel 
through each other, i.e. V{xg,x) = a5{x — Xg), a — >• oo. This potential is a good 
approximation to real interaction potentials if the wave length of the particles is large 
compared to their interaction range. In regards to the topic of this part of the thesis, we 
call the two colliding particles the Brownian particle and the gas particle, respectively. 
In the following, the index g will be used for the gas particle and no index corresponds 
to the Brownian particle. The two particle Hamiltonian is then given by 

H = + ^ + a6(x-Xq), a^oo, (12.3) 
2mg 2m 

where x and p are the position and momentum operators for the particles. 

We will find that if the width of the two Gaussian wave packets relate in a certain 
way, then the collision does not produce any entanglement. In this case, the mean 
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positions and momenta of both particles behave hke in a classical coUision, and the only 
non-classical feature will turn out to be the usual spreading of the wave packets. These 
features will enable us to think almost classically of QBM, which in turn will guide us 
in setting up the QBM master equation. 

There are two approaches to tackle the described collision. The first and easier one 
is along the lines of Hornberger and Sipe [H], but generalized to finite mass ratios of 
the colliding particles. One decomposes the Gaussian wave functions into momentum 
eigenstates \pg) and \p), and applies scattering theory to their collision. For particles 
moving only in one dimension, the only possible scattering operator is of the form 



which is easily seen from conservation of energy and momentum. The parameter a = 
m/nig denotes the mass ratio of the colliding particles. For the hard core interaction, 
one finds that the phase ip{p,pg) = vr is independent of the momenta. For completeness 
and for a first taste of what follows, we start with such a scattering calculation in the 
first section of this chapter. 

In scattering theory, a complete collision of the two Gaussian wave packets is pos- 
tulated. But there are situations in which the wave packets collide only partially, in 
particular when the velocity uncertainty of the colliding particles is comparable to their 
relative mean velocity. Therefore we will employ the more complicated approach of 
solving the time dependent Schrodinger equation. This will enable us to precisely link 
the situation in which complete collisions occur with the high temperature and low 
denisty limit. Furthermore and contrary to the scattering theory calculation, this ap- 
proach enables us to look at the two particle state not only before and after the collision 
takes place, but also during the collision event. In position representation this will serve 
as proof that no position jumps take place in collisional Brownian motion, whereas in 
momentum representation we will find the expected momentum jumps. 

12.1 Scattering theory calculation 
12.1.1 Some remarks about scattering theory 

Scattering theory is a very powerful tool for the following type of question. Assume 
at time t = we have two sufficiently far separated particles whose initial state can be 
approximated by a product state (asymptotic-in state) g{t = 0) w gin = pi®P2- Because 
of their mutual separation, one can neglect their initial interaction which justifies the 
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approximation of an initial product state, but if they travel towards each other they will 
interact at some time. What is their state sufficiently long after their interaction? 

Quantum scattering theory provides an answer to this question without looking at 
the actual interaction time. For this one introduces the two particle scattering operator 
S which maps asymptotic-in states to asymptotic-out state^ 



Qout 



SQinS^ (12.5) 



At no time is Qout the actual two particle state, but if the time t = t is such that the 
two particles do not interact any more, the actual two particle state can be obtained 
approximately by using the free particle evolution operator for each particle 

q{t) « Ul{T)0U2{T)QoutUl{T)(S)Ul{r) 

« Ul{T)®U2iT)SQiO)S^UliT)^Ul{T). (12.6) 

The approximation is good if the interaction of the particles at times t = 0,t can be 
neglected. 

The scattering operator preserves the energy of the two particle state, and as such it 
commutes with the two particle free evolution operator [Ui{t)<^U2{t), S] = 0. Therefore, 
we get the same result if we first apply the free particle evolution operators and then 



the scattering operator as is shown in Fig. 12.1, The calculation in this section is along 



Fig. pJ] (b). 

The popularity of scattering theory is due to its simplicity compared to calculations 
using the Schrodinger equation. On the other hand, scattering theory calculations have 
several drawbacks in that they do not tell us 

• what happens during a collision? 

• what is the collision time tc (sometimes called interaction time)? 

• if the collision is complete {tc < oo)? 

• if a collision occurs at all? 

In typical applications of scattering theory, one prepares an incoming beam of particles 
which scatters off a targelj^ and one is interested in the momentum distribution of the 

^Of course, the scattering operator depends on the interaction potentiaL But once the scattering 
operator is known, either from experiments or from theory, conclusions can be drawn without considering 
the microscopic details of the interaction potential. 

^Often two beams of particles are prepared, which scatter off each other. 
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Figure 12.1: Different representations of the same physical process. Each time the same initial 
two particle state {t = 0) evolves into the same final state {t — t). Blue and red lines represent 
free evolution of the respective particle and arrows represent the action of the scattering operator, 
(a) The initial states evolve freely up to the interaction time, and afterwards continue their free 
evolution towards the final state, (b) The scattering operator acts on the initial state and the 
transformed state undergoes free evolution, (c) The initial state evolves freely until t = r, and 
then the scattering operator is applied to yield the final state. 

particles after they collided with the target. The beam is prepared long before it hits the 
target, and the momentum distribution is measured long after the target was hit, often 
by using a screen. In such a setup one is not interested in any of the above questions, 
which explains the huge success of scattering theory. 

The situation is different in QBM, where we aim to derive the time resolved dynamics 
of a Brownian particle. In principle, once the approximation of instantaneous collisions is 
applied, it should be possible to derive a master equation for the Brownian particle only 
using scattering theory. However, such approaches lead to some questionable results in 
the past. Therefore, we choose to instead apply Schrodinger's equation to gain a deeper 
insight in the scattering event, eventually leading to a better understanding of QBM as 
well as the conditions under which we can expect Markovian behavior. 

Nevertheless, we present a scattering theory calculation in this section. On the one 
hand this serves as confirmation of the later Schrodinger type calculation, on the other 
hand it indicates where problems might occur in such an approach. We will assume an 
initial product state of gas and Brownian particle, both being described by a Gaussian 
wave packet. 

Although the incoming gas particle is assumed to be reasonably localized in momen- 
tum, we will see that the outgoing Brownian particle state depends on the precise nature 
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of the incoming gas particle state. In fact, from linearity of the scattering operator as 
well as the partial trace, we fincj^ 



1 



(12.7) 



Therefore we can approximate a gas state \ipg) hy IV'g) if their overlap is close to 
one, without changing results for the Brownian particle. But because the overlap 
Wg{xg,Pg\ Xg,pg)w^ of two Gausslau wave packets [Eq. (12.1)] approaches zero if Wg 
goes to infinity (for any fixed value of Wg), it can not be guaranteed to find correct 
results by approximating the incoming gas particle state \xg,pg)wg with a momentum 
eigenstates (which is liuiw^^oo \xg,Pg)w^ up to normalization). 

If one is only interested in the momentum distribution after a collision, it might be 
sufficient to consider incoming gas particles in momentum eigenstates. But as we are 
interested in the entire Brownian particle state, including the coherences of different 
momenta as well as the position, this can no longer be assured and one has to consider 
the actual incoming gas particle statej^ 

12.1.2 The calculation 

We use an asymptotic-in state in the following momentum representation 



\x,p)w ^ \Xg,Pg)Wg 



9_^i{xp+XgPg)/2h 



Applying the transformation Eq. ( |12.4 ) and substituting 

_ (1 - a)p' + 2p'g _{a- l)p'g + 2ap' 



1 + a 



l + a 



(12. 



:i2.9) 



we find 



1^' 



out/ 



^yww. 



9^i{xp+XgPg)/2n UudUg \u)\Ug) 



x[{1 — q)u + 2ug\ + Xg[{a — l)ug + 2au] 



X exp 



(1 + a)h 



■^This is correct for any decent distance measure of density matrices. 

^Of course, one can write J dpg tp{p) |p5)for the incoming gas particle state and perform the calculation 



in the momentum basis, as is done in the next subsection. 
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X exp 



W^[{l-a)u + 2ug - (l-a)]p2 W^[{a-l)ug + 2au - {l-a)pgf 



2{1 + a)^h^ 



2(l + a)2;i2 



(12.10) 



To see whether this is a product state, we now separate u and Ug in the last exponential. 
To this end we find a term 



exp 



A{l-a)uug {W^-aW^) 
2(l + a)2^2 



(12.11) 



which shows that u and Ug can in general not be separated, and therefore, Qout is not a 



product state. But it is also apparent from Eq. ( 12.11 ) that the entanglement of the two 



outgoing particles is lifted, if the width of the colliding wave packets relate according to 
their masses 



maWi 



(12.12) 



This relation will hugely simplify the calculation as well as the physical interpretation, 



and therefore the (over) complete basis |x,p)vF with W satisfying Eq. (12.12) is the 
basis of choice for studying a single collision with a gas particle state \xg.,Pg)wg- If the 
initial state of the Brownian particle is not of the form then is can be written 

as a superposition of such states and one can still use the results from the simpler 
calculation with For such an initial state, of course, the two particle state 

after a collision will show entanglement. An important example is a so called cat-state 



\'X'a^Vo)'w + \'^h^Vh]'W^ which is presented in section 12.7 
Next we introduce 



V 



2axg + (1 — a)x 
1 + a ' 
_ 2pg + (1 - a)p 



2x 



:i 



Q jXn 



1 + a 

and after replacing u by p' we find 



P9 



1 + a 
2ap — (1 — Oi)pg 



1 + a 



(12.13) 



1^' 



out/ 



9_^i(xp+XgPg)/2h 

x^dpMp;e-(^p'+^«p«)/^^(p'-^^)'^'+K-p-.)^'^«1/2r^^ |p'>|p;^> 

-\X,P)W ®\Xg,pg)Wg (12.14) 



This surprisingly simple result shows that all the collision does, is to change the position 
and momenta of the particles according to Eq. (12.13). We want to point out that this 
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does not mean that the position necessarily changes during the cohision process, because 
as discussed in the previous subsection, the state l^'out) is not the actual state at any 
time. To find the physical state at time t = t sufficiently long after the collision, we 
need to apply the free particle evolution operators as in 

\^{t)) = U{T)\x,p)w(^UgiT)\Xg,Pg)w,- (l2.15) 

Taking into account the free particle evolution operators, which shift the position of the 
wave packets according to their momentum, we will find in the following sections that 
the position distribution follows a continuous flow. 



12.2 Eigenfunctions 



The following sections will be concerned with solving the two particle Schrodinger equa- 



tion Eq. (12.3). The first task will be to find the eigenstates of the Hamiltonian, which 



will be done in position representation. The time independent Schrodinger equation 
reads 

Two orthogonal sets of solutions are found with the ansatz 



^E,jx'g,X') 

'^Etotix'g,x') 



jH^'+»^'3)sm{k{x'-x')), 



Jk{x'+ax' ) , 



kGR, k£R+ (12.17) 
keR, keR+. (12.18) 



"s-* sm{k\x'g - x'\), 

Note that in the second set the relative coordinate is replaced by its absolute value. 



Substitution in Eq. (12.16) gives 



a 
Etot 
E 

E 



m 



E + Ea = E + E 



{l + a 



(1 + a 



k^h^ 

I 

2m g 

. 

2m 



(12.19) 
(12.20) 
(12.21) 

(12.22) 



We will refer to the eigenfunctions Eq. (12.17) and Eq. (12.18) as "^fr and ^'|r, where 



the superscripts stand for antisymmetric and symmetric. With 



[dx'g dx'J 1^'^ 



hk(l + a)^^^' 

kk 



(12.23) 
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these states are eigenstates of the total momentum operator ptot = Pg + P and have 
therefore definite total momentum 



Ptot = hkil + a). (12.24) 

es c 

{ap — pg)/mg with eigenvalues 



Similar, are also eigenstates of the square of the relative velocity operator Vrei 



"rel 



1 +l\\2~k\ (12.25) 



12.3 Initial state 

In this section, we construct the appropriate initial state for the wave packet scattering. 
The problem is that the most obvious candidate for this state, 

\^{t = 0)) = \xg,pg)0\x,p), (12.26) 

does not satisfy constraints implied by the interaction between the two particles, so a 
more careful derivation is required, and is detailed below. 



The state Eq. (12.26) can be written in position representation 



„ — i{XgPg+Xp)/{2h) 

Note that this state is only a good approximation to the real two particle state if the 
Gaussian wave functions have little overlap. To make this statement precise, we define 
the overlap as 

dx'Mt = ^,x'x' = x'f = , ^ . expf--^!^) (12.28) 



from which we find the condition for small overlap 



\xg-x\ > ^Wj^W^. (12.29) 

To simplify the calculations, we will use the centre-of-mass reference frame, and 
assume the gas particle approaching the Brownian one from the left side, i.e. 

= -p > (12.30) 
axg = -X < 0. (12.31) 
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The task at hand is to expand Eq. ( 12.27) in the eigenfunctions and ^fr- Here we 
encounter a problem as the two particle state Eq. (12.27) does not satisfy the boundary 
condition ^'(0,x^ = x'). Hence it is not possible to expand this state into a linear 
combination of and ^I'|r, and we have choose a more elaborate approach to find an 
initial state which is close to Eq. (12.27) and has the property ^'j(0,Xg = x') = 0. 
This will prove quite tedious, but the very nice result of an almost classical collision will 
make this effort worthwhile. 



To do so, we first expand Eq. (12.27) into a set of functions which is related to the 
set energy eigenfunctions by substituting the sine function with an exponential function. 
In a second step we will remove the overlap of the two Gaussian wave packets and by 
doing so we will recover the energy eigenfunctions. 

We start by comparing two functional expansions, one being the two-coordinate 
Fourier transform, the other one being closely related to our eigenfunctions by replacing 
the sine function with an exponential function: 



*(0,x' x') 



Here we introduced x' 
is easily found to be 



1 
1 

2^ 



x' + ax'g and x' 



dkg dk e*'=«^9e^'=^'^'(0,/fcg,fc) 



dk dk e 



ikx' ikx' 



^{0,k,k). 



(12.32) 
(12.33) 



- x' . The Fourier transform of Eq. ( 12.27) 



^(0,^^,^) 



^g'^ ^i(x,p„+xp)/2h^-i(x,k„+xk)^-W^(k,-^-f)V2^ 
IT 



whereas Eq. (12.33) can be written as 

1 

2^ 



■^{0,x'g,x' 

Now we substitute 



dk dk e^Ck+»kW,e^Ck-kW ^(0, k). 




With the determinant of the Jacobi matrix 
J = 




l+a 



(12.34) 



(12.35) 



(12.36) 




det(J) 



l + a 



we find 

^(0,4,x') 



1 

2^ 



dkgdke'^^'^'^e'^''' 



1 



l + a 



VI/ I 0, 



ak kg -\-k 



l+a l+a 



(12.37) 



(12.38) 



108 



CHAPTER 12. SINGLE COLLISION 



Comparison with Eq. ( 12.32[ ) then results in 

^(0, k, k) = {l + a)5 (0, k + ak, -k + k) (12.39) 
where we can substitute Eq. ( |12.34 ) to finally find 



m{0,k,k) = (i + a)y^^^^e^(^«?'«+^P)/2/ig-ifc{x,-x)g-ifc(axg+^^ 

X expl-h\w' + W^) + k(^^-?^] -kk{aW^-W')\. 



(12.40) 



In the centre of mass reference frame (Eq. ( 12.30 ) and Eq. ( 12.31 )) this can be simplified 
to 



^{0,k,k) = (i + a)y^^^V^^e*^P(i+")/2"V^^{i+")/a 

X exp |-^P(a2vF'2 + w'^) + k[k + ^) {W'^ - aW^] 



X expj-^l^fc + lj (Wg^ + W')'j. (12.41) 

Let us have a closer look at this equation. We see in the first line, that the phase is 
independent of k which tells us that the expectation value of the position of the centre 
of mass is zero, as should be expected in the reference frame of the centre of mass. In 
the second line (in combination with the third line) we find that the expectation value 
of the total momentum is also zero, as it should be in the reference frame of the centre 
of mass. Furthermore, the second term in the second line describes an entanglement of 
the values k and k (there is no entanglement between gas and Brownian particle, but 
entanglement of the total momentum and the relative velocity), which is lifted if and 
only if 

^2 ^ ^^r2_ (12.42) 
This is the situation we are interested in, as it will turn out to be the condition for 



Brownian and gas particles not being entangled after the collision. Using Eq. (12.42) 
we can now write 

%{0,k,k) = (l + a)^Le'^P(^+°^/2ay^a;(l+a)/a 

y ^^^/a 

X e^j>{-l~k\l + a)W'-l{k+^y^-±^wA. (12.43) 
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Fi gure 12.2: Position representation |(a;g| a;gPg)p (dashed) and |(a;|a;p)p (solid). The overlap 
is not physical and has to be removed before the initial state can be written as a superposition 
of energy eigenstates. 



With Eq. (12.33) and Eq. (12.43) we have constructed an initial state Eq. (12.27) 



which is of the form seen in Fig. 12.2 As pointed out earlier, there is some probability 



of finding the gas and Brownian particle at the same position which is not physical for 



the hard core interaction. Although, if relation (12.29) is satisfied this overlap is small, 



it is non-zero and hence it is not possible to decompose this state into energy eigenstates 



Eq. (12.17)-(12.18). Therefore we have to construct a state \^i{t = 0)) which is close to 
\xgPg)i^ \xp) for \xg — x\ ^ Wg + W , aud has the property 

^-(0, x'g, x) = for all x' = x'g. (12.44) 

As a first step we use the antisymmetric superposition 




X 



Figure 12.3: Position representation of l^'sup) (gas: dashed. Brownian: solid). There is no 
overlap and therefore this represents a physical state. The black parts represent the gas particle 
coming from the left, whereas the gray part represents the gas particle coming from the right. 



|^S«p)= C{\XgPg)l^ \XP)- \-Xg,-Pg)l^ \-X,-p)) 

which automatically fulfills Eq. ( |12.44 ). The normalization constant 



C 





1 + a / 


2 - 2exp 





+ 



(12.45) 



(12.46) 
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approaches -v/l/2 if initially Brownian and gas particle are sufficiently separated. With 



Eq. (12.43) it is an easy task to find 



^supiO,~k,k) = 0(1 + a) 



W 



^ixp{l+a)/2ah 



_^-ikx{l+a)/a 



hJ 2a 



(12.47) 



Substituting in Eq. (12.33), we find that the cos(A;x') cancels and therefore "^supiO, k, k) 



expands |^'s„p)into energy eigenstates 

i r°° _ r°° ~ --, ~ — 
^sup{0,x' x') = - dk dke'^'' sm{kx)'^supiO,k,k). 



(12.48) 



The possibility to decompose ^'^^^(O, x', ) into energy eigenstates is a direct conse- 



quence of Eq. (12.44). The state \^sup) can be pictured as in Fig. 12.3 (this is not 
completely correct since in this state Brownian and gas particle are entangled). To 



achieve Eq. (12.44) we paid the price that the gas particle does not approach the Brow- 



nian particle from one definite direction. In order to return to the original picture where 
the gas particle approaches the Brownian one from the left, we want to get rid of the 
parts in Fig. |12.3| indicated in grey. This is simply done by adding the symmetric set of 



eigenfunctions Eq. (12.18) to the antisymmetric one: 

dk e'^^' (sin(fci') - sm{k\i'\)] ^,„p(0, k, k) (12.49) 



'^i{{),x' x) 



i 



dk 



oo 
oo 



dA; 



dk Ul-Ax'x')-m-Ax'x')) ^sup{Q,~k~k)- (12.50) 



This state has all the required properties and is the initial state we will use in the 
following. 



There are of course other ways to construct initial states which fulfill Eq. (12.44) and 



are approximately the tensor product of two Gaussian states. But the choice here has 
the appealing property that \^ sup{i = -co)) is an exact product state as can be seen 
in the next section. This should be required since the two incoming particles should 
not be entangled a,t t = — oo. This property is of course lost for all finite times, but if 
|xg — x| ^> Wg + W then the two particle state is still approximately a product state at 
t = 0. 
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12.4 Time evolution 



Since '^"f!^ are eigenfunctions of the Hamiltonian with energy Eq. (12.20) the time 
evolution is an easy matter: 



^ i ^ q 1 ) 



^/27r 



dA; / dk e 



'iEtotik,k)t/h 







^-^ix'g, x') - ^-^{x'g,x')^ ^sup{0, k, k), 

(12.51) 

where Etot{k, k) is given in Eq. (12.20). To perform the integrals for x' > x'g, we use the 
exponential instead of the sine (remember that the cosine cancels upon multiplication 
of the exponential with ^'s„p(0, k, k)) 




k- I-oc I-oo dfe e-'^^o^Ck^^t/h^ikx' ^ikx' ^^^p(o, fc, k) for x'g < x' 



'^i{t,x'x') 



and find after a straightforward integration 



for x'g > x' 



(12.52) 



a 



X exp 



exp 



. XgPg + XP 



iht ^ 



exp 



exp 



) {w^ + 

( pW* _ 



xht 
m 



J 



xht 

m 



(12.53) 



for x' < 0. Let us have a closer look at the two terms in curly brackets. The first term 
dominates if (2; + ^) > 0, which is the case for t = because we chose x > at the 
beginning of the calculation. The second term dominates for t — )• 00 as p < 0. That the 
smaller term vanishes for t — )• ±00 requires a more detailed comparison which can be 
found in the appendix of this chapter. It turns out that the smaller term only vanishes 
if 

That is, the absolute value of the mean momentum of each particle has to be greater 
than its momentum uncertainty. This is actually quite intuitive, as otherwise the wave 
packets would spread faster than they move towards each other, and parts of the wave 
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packets would not collide at all, which we will refer to as "incomplete collisions" . If this 
condition holds we are also able to define a collision time (see appendix) 



1 + a 



rrin 



(12.55) 



as the time it takes from the first term to dominate to the second term to dominate. 



Let us assume that condition (12.54) is fulfilled. The first term in the curly brackets 



can then be neglected sufficiently long after the collision and the two particle state can 
be brought into the form 



Pa 



,e-iP9^V'^exp 



1 



iht 



rUgW^ 



■nig j 



e 2ft V 



iht 
Wm 



iht 



ni 



2iW^ + 



(12.56) 



where C 



as well as x' 



x' and x' = ax'g + x' were used. 



This is a surprising simple result since it is a product of a function of x'g and a 
function of x' and therefore the position representation of a product state. Even more, 
these are the Gaussian states 



\^{t))= Ug{t) \-Xg, -Pg)® U (t) \-X, 



(12.57) 



where Ug{t) and U{t) are the unitary evolution operator for a free gas and Brownian 
particle, respectively. Of course, this equation also applies if the gas particle approaches 
the Brownian one from the right. Therefore, in the centre-of-mass reference frame, all 
the collision does is mirroring the initial values Xg, pg, x, and p which is exactly what 
happens in a classical collision. The only extra requirement in a quantum collision 



is that the widths of the colliding wave packets fulfill Eq. (12.42). This result agrees 



with Schmiiser and Janzing [75j who are looking at entanglement creation during the 
collision of two hard core particles in one dimension. They find that the collision does 
not produce any entanglement if the width of the Gaussian wave packets relate to their 
masses like mW"^ = nigWg. 



Now we briefly return to the remark at the end of the last section. If Eq. (12.54) 



-oo. 



holds, then the second term in the curly brackets in Eq. (12.53) vanishes for t - 
We then find 

\^it))=Ug{t)\Xg,Pg)0U{t)\x,p), (12.58) 
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which is a Gaussian product state as promised. 



12.5 General reference frame 

So far we used the centre of mass reference frame. In this section we turn to a coUision 
of the two particles with an initial state 



\qf{t = 0))= \xg,pg)(g) \x,p) 



(12.59) 



with general Xg, pg, x, and p. It would seem that we could just change the reference 
frame as in classical mechanics, but that would not necessarily give us the phase of the 
final state. Therefore we will perform the change of reference frame in a more formal 
way by using the Glauber displacement operators Dg{Axg, Apg) and D{Ax, Ap), such 



that the state Eq. (12.26) transforms to the state Eq. (12.59) 



= 0))=Dg{AXg, Apg)D{AX, Ap) \Xg,Pg)\x,^. 



(12.60) 



To avoid confusion, in this section we use a wide tilde to denote quantities in the centre- 
of-mass reference frame. 

Since we have a number of 'x'-symbols, it is worth summarizing them at this point 
to avoid confusion. The 'p'-symbols are defined in an analogous way. 



gas particle 



Brownian particle 



position operator 
initial position of wave packet 
initial position of wave packet 
in centre of mass reference frame 

X — X 

variable in position representation 

X I G!.X g 

X'g — X' 

position of wave package induced 
by collision 



X 

X 

X 

Ax 

x' 



X 
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The following relations will be useful: 



Xg ~\~ 


— Xg 


X + Ax 


= X 


Pg + Apg 


= Pg 


p + Ap 


= P 




= X 


-Pg 


= P 


AXg 


= Ax 


^Pg 


= a Ap 



Ax„ 



Xg 
X 

Ax 

Pg 
P 

Apg = a Ap 



l+a 



a 



l+a 
x+aXg 
l+a 

l+a 

Oip-Pg 
l+a 

Pg+P 
l+o 



(12.61) 



a 



The first four relations on the left hand side follow directly from Eq. (12.60), the next 
two are the properties of the centre-of-mass reference frame, and the last two are because 
D = Dg{Axg, Apg)D(Ax, Ap) represents a change of reference frame. Actually, it is 



because of the last four relations that there is no additional phase in Eq. (12.60). For 
given initial values Xg, pg, x, and p, we can use the right hand side of Eq. (12.61) to 
find the appropriate values Axg, Apg, Ax, and Ap as well as Xg, pg, x, and p to use 



in Eq. (12.60). 



Let us return to the time evolution 

\^{t))= e-'^^/''\^{0))= DD~^e'^*/''D\^{0)y (12.62) 

where we also used a wide tilde to denote the two particle state in the centre-of-mass 
reference frame. We therefore need to calculate the operator 



exp 



.PgAXg—ApgXg . f) Ax — Ap X 



X exp 



X exp 



a2 p2 



it Pg , P 



h \ 2mg ' 2m 



■ ^Pg ^g ~ Pg '^^g _j_ • ^P x — p Ax 



= ei,Ty\B + [A,B] + -[A,[A,B]] + 
In the last line we used the Baker-Campbell-Hausdorff formula 



(12.63) 
(12.64) 



exp(x) exp(y ) = exp ( X + y + ^ [X, y] + 1 [X, [X, y]] - ^ [y, [x, y]] + • • 



(12.65) 
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With the last two relations of the left hand side of Eq. ( 12.61 ) the following commutators 
are found 



ih 



[A,B] 



PgApg pAp , 

"I ~ ^ ~ AXg)ad [X — Xg) 



m„ 



m 



PgApg _^ pAp 



ih 
~t 



[A,[A,B]] 



nig 



m 



Ap'g 

H 



(12.66) 
(12.67) 



where 6' denotes the derivative of the delta function in its argument. Substituting the 



commutators into Eq. (12.64), we arrive at 

-it ({pg + Apgf {p + Apf 



exp 



exp 



h 



2m„ 



+ 



2m 



+ aS{xg — x) 



-it I Apl 
h \ 2mg 2m 



exp 



-it ( Pg Apg ^ pAp 



h V m„ 



m 







exp 





(12.68) 



where (a Ap — Apg)5'{x — Xg) = was used in the second line. Since we know from the 



previous subsection that g-*m/fi ^(0)^= k'(t)^with \^{t)J given in Eq. p2.57L this 
operator can now be substituted in Eq. ( 12.62| ) to give 



\^{t)) 



exp 



-it f Apl ^ Ap^^ 



h \ 2mg ~^ 2m 



X Dg{AXg,Apg)Dg 



0]D{Ax,Ap)D 



ma 



Dg ( Axg + ^-^^,Apg ] D [ Ax + ^— ^ , 



mn 



tAp 
m 



tAp 



m 



(12.69) 



Next we need to determine the effect of the displacement operator on a single particle 
wave function. For that we apply the operator to the position eigenstates 



t){x,p)\x') = e^^P/2V^^P/V^P/^ 

_ ^ixp/2h^ix'p/h^^i _^ 



(12.70) 



Therefore the transformation of the position representation of a wave function can easily 
be shown to be 



(12.71) 
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By substituting the position representation Eq. (12.56) into Eq. (12.69) and by using 



Eq. (12.71), 1^(0))= \xg,pg)(g) evolves to 



\'^{t))=Ug{t)\Xg,Pg)^U{t)\x,p) 



(12.72) 



where x and p are given in terms of the initial values x and p (using Eq. (12.61)) as 
follows 



2x — (1 — a)xg 



Pa 



l + a 




2ap - (1 - 


a)pg 


l + a 




2axg + (1 - 


- a)x 


l + a 




'iPa + (1 - 


a)p 



l + a 



(12.73) 
(12.74) 
(12.75) 
(12.76) 



As the evolution operator U (t) shifts the position of a wave packet according to 
x(t) = X + pt/m (in addition to spreading the wave packet), these equations show that 
the centres of the wave packets in position as well as momentum space behave precisely 



the same as in a classical collision (see chapter 9.1). This remarkable result will allow 



us to think almost classically when setting up the QBM master equation in chapters 14 



and 15 The only difference will be that the interaction of the colliding particles is not 



instantaneous. 



12.6 Momentum and position jumps 



After we have examined a single collision, we are in the position to discuss the appearance 



of position and momentum jumps (recall the definition of a jump from section 9.6 ) during 
the collision process. The two particle wave packet during the collision in position 



representation is given in Eq. (12.53) (without loss of generality we use the center- 



of-mass reference frame in this section). In the first term the Brownian momentum is 
localized around p and in the second term around —p. During the collision the first term 
decreases continuously while the second term increases. As the momentum distribution 
is effectively zero between p and —p, this process accounts for momentum "jumps" 
(although the jump is not instantaneous), and if these occur at random we are lead 



to momentum diffusion. In fact, already Eq. (12.25) shows that the absolute relative 



velocity of gas and Brownian particle is a constant of motion (this is a consequence of 
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the delta type interaction potential). Therefore, during the collision the eigenvalues for 
the relative velocity can swap the sign, but not change continuously. 

The situation is different in position space, as during the collision (x + — ~ 0) both 



terms of Eq. ( 12.53 ) are localized around the same position. To be more precise, we have 



a look at the position probability distribution p{t, x') of the Brownian particle which we 



obtain from Eq. (|12.53|) 

p{t,x') = 



dx'g \^{t,X'g,X')\ 



1 



vr + 



+ Er 



ax' + X + 



Er 



pt 



ax 



pt 

m 



a(W^ + 



exp 



(x' + X 



pt ^ 



J2t2 



a 



1^2 + 



exp 



Ty2 + 



_ 2*V 

m ) 



(12.77) 



Here we defined Er{x) = -j^ J- qq e *^ which is closely related to the error function 

u hi rok : ...:_x_J _„/ , U r..„., J Wa 



and appears because Eq. (12.53) is restricted to x'g < x' . We further used pg ^ ^ to 



avoid the interference of the two terms in Eq. ( 12.53 ), but this is only done to simplify the 



expression and not necessary to avoid position jumps. A plot of the position probability 

(a), and reminds of a 



distribution over time during a collision can be seen in Fig. 12.4 



classical collision. The probability distribution "flows" continuously in position space 
and no jumps can be associated with such an behavior. As such, random collisions can 
not result in position diffusion according to (Ax'^) oc t on the short time scale of few 
collisions. 

This has to be compared with the momentum probability distribution p{t,p') for 



which we again use Eq. (12.53) 



pit,p') = J dx'g \^{t,x'g,p')f 
^{t,x'g,p') = -j=j dx'e-'P'-''^^{t,x'g,x'). 



(12.78) 



The lower bound of the integration is because Eq. (12.53) is only valid for x' > x' and 



^'(t,x',x') is zero for x' < x' . The integrations are rather lengthly and involves non- 



analytic functions. We therefore only show the resulting plot in Fig. 12. 4| (b) obtained 
numerically. We see that the momentum distribution does not "flow" from the initial 
value towards the final value, but rather continuously decreases at the initial value and 
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Figure 12.4: Position (a) and momentum (b) probability distributions for the tracer particle 
during a collision in the reference frame of centre of mass. While the position probability 
distribution "flows" , the momentum distribution "jumps" from the initial value to the final 
one. This disproves QPD. The solid lines in (a) are the corresponding classical trajectories. 
The collision time indicated in (b) is taken from Eq. (12.55). Parameters in the centre of mass 
reference frame are: x = 10, p — —2, m — 1, a — 0.3, h—1, and W — 4. 



increases at the final value. It is this sort of behavior, which, if occurring at random, 
results in a diffusion process according to (Ap^) cx t. 



12.6.1 Discussion of position diffusion 

In light of the conclusions of the previous subsection, one might wonder why all pro- 
posed QLBE exhibit QPD. In particular the necessity of QPD to achieve positivity of 
any Markovian QLBE seems to challenge our results at first glance. This apparent con- 
tradiction can however be resolved by taking into account that any Markovian master 
equation is only vahd on a coarse grained time scale 5, which has to be large compared 
to the collision time tc, such that a collision can be considered as an instantaneous event. 
During a time interval 6, a collision does not only change the momentum distribution, 
but also significantly influences the position distribution as is indicated by the arrow 
in Fig. 12. 5[ That is, looking only at a coarse grained time grid, a collision results in 
apparent position jumps, therefore leading to fictitious position diffusion. 

This is by no means a quantum feature, as also classical dynamics lead to position 
diffusion if they are derived using a coarse grained time scale. In particular Smoluchowski 
|44j used the damping time as coarse grained time scale to derive Brownian motion, one 
of the most fundamental diffusion processes in nature. However, in classical dynamics 
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Figure 12.5: During a time interval S, there can either be a collision (red) or no collision (green). 
On a coarse grained time grid a collision results in a discrete change of the position distribution 
(indicated by the arrow) compared to a no collision trajectory. As this discrete change occurs 
at random, it gives raise to position diffusion, ft is also clear from the picture, that the same is 
true for classical trajectories. 



it is very much accepted that position diffusion results from a coarse graining time 
approximation and is not apparent on the short time scale of momentum diffusion. 
In fact, if the interaction between gas and tracer particle is of the hard core type, 
it is possible to derive classical Markovian dynamics without using a coarse grained 
time approximation [57] (see also section [9]). The resulting Fokker-Planck equation for 
the phase space probability distribution (also called Kramers equation) then exhibits 
momentum diffusion, but not position diffusion. 

The wave nature of quantum particles forbids instantaneous collisions even for hard 
core interaction potentials, and therefore the introduction of a time coarse graining 
approximation is necessary to derive Markovian quantum dynamics. For the same reason 
as for classical calculations, the coarse grained time scale then leads to fictitious position 
diffusion on the short time scales. But again, it has to be realized that this contribution 
to position diffusion is not a physical process, even if any possible QLBE would indicate 
so. 



120 



CHAPTER 12. SINGLE COLLISION 



12.7 Scattering off a cat-state 



To get a first qualitative feeling about the process of collisional decoherence, we consider 
a gas particle in the state \'>pgiO)) = \xg,pg) to be scattered off the Brownian particle, 
which is initially in a superposition state \tp{0)) = \xa,Pa)+ \xb,Pb) (not normalized), 
known as a cat state in the literature. Because the gas particle state after the collision 
depends on the initial position and momentum of the Brownian particle, we expect that 
the collision creates entanglement if the Brownian particle is initially in a cat state. 



Using Eq. (12.72) - (12.76) we find for the two particle state after the collision 



= Ug{t)\Xg{Xg,Xa),PgiPg,Pa))'^U{t)\x{Xg,Xa),p{pg,Pa)) 

+ Ug{t)\Xg{Xg,Xb),Pg{Pg,Pb))^U{t)\x{Xg,Xb),piPg,Pb)), (12.79) 



where x{xg,Xa), ■ ■ ■ are given in Eq. (12.73) - (12.76). Building the two particle density 



operator and tracing out the gas particle we find the Brownian particles density operator 

p{t) = U{t)\x{Xg,Xa),p{Pg,Pa)Xx{Xg,Xa),p{pg,Pa)\U^ 
+ U{t) \x{Xg,Xb),p{Pg,Pb)\x{Xg,Xb),p{Pg,Pb)\ U'^ 
+ U{t) \x{Xg,Xa),p{Pg,Pa)Xx{Xg,Xb),p{pg,Pb)\ 

x{Xg{Xg,Xb),Pg{Pg,Pb)\ Xg{Xg,Xa),Pg{Pg,Pa)) 

+ h.c, (12.80) 
where h.c. denotes the Hermitian conjugate. We see that due to the collision, the 



diagonals get shifted according to Eq. (12.73) - (12.76). That also means, that the 



two humps in the position probability distribution get closer to each other, as is shown 



in figure 12.6 One might think that after a large number of collision, both humps 
will eventually be on top of each other, which would contradict uniformity in space. 
However, this apparent problem is resolved by taking into account gas particles which 
are initially located between the two humps, therefore only colliding with one of them. 
When constructing up our master equation, we will therefore have to be careful with 
the collision statistics, which will be determined by a rate operator. 

In addition to a change of the distributions of momentum and position, a collision 
also results in a multiplication of the coherences by 

{Xg{Xg,Xb),Pg{Pg,Pb)\Xg{Xg,Xa),Pg{Pg,Pa)) = 



exp 



2a{xbPa - XaPb) + (1 - a)Pg{Xa - Xfe) - a(l - a)xg{pa - Pb) 
' {l + aYh 
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Figure 12.6: Position distribution of an initial superposition of two Gaussian wave functions, 
both having the same mean momentum, but different mean positions. The lower hump experi- 
ences the collision first and then moves towards the upper hump, until this one also experiences 
the collision. The white lines correspond to the mean positions of the colliding gas particle. 



X exp 



a 



+ 



(12.81) 



The collision has two effects on the off-diagonals. In the first line of the right hand side, 
there is a relative phase of the two wave packets, which depends on the state of the 
gas particle. In the second line, decoherence occurs in both, position and momentum, 
which reflects the post collisional entanglement of the two particles. Later, when we 
discuss a measurement interpretation of Brownian motion, the second line corresponds 
to decoherence due to a measurement, which the gas particle performs on the Brownian 
particle. 



It is important to realize that, if the state of the gas particle is taken randomly 
from a thermal distribution, also the first part will contribute to decoherence due to 



phase averaging and, in fact, in chapter 17 it will turn out to be by far the dominant 
contribution. 
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Appendix 

To compare the two contributions of Eq. (12.53), we take their L^-norm defined by 



dx'g dx' \f{x'g,x')\'^. 



(12.82) 



Rewriting the absolute value of each term of Eq. (12.53) we find 
2C\/S 



vr Ty2 + 



'2 



dx' dx' exp < 



x'>x' 



l+g 



n 



(1 + a) J- 



dx / dx' exp < 



2C^ 



\ m 



dx e 



(12.83) 



where the plus corresponds to the first term in the second line of Eq. (12.53) and 
the minus to the second one. We used the transformation of the integration variables 
according to 



l+a 
x' —ax' 
l+a 



det J 



l + a 



(12.84) 



where J is the Jacobi matrix. A term is small if the upper integration boundary is 



significantly negative. As x > 0, the second term in Eq. (12.53) is small at t = exactly 



if Eq. (12.29) is satisfied. For t — )• cxd we find that the first term vanishes exactly if 

(12.85) 



p > 



a h 



l + aW 



If inequality (12.85) does not hold, then the relative velocity of the wave packets is not 
large compared to their spreading, and therefore, there is a non-vanishing probability of 



having no collision at all. For later purposes it is useful to put Eq. (12.85) into a form 



which is invariant in respect to Galilean transformations. Relation (12.85) was derived 
in the center-of-mass reference frame in which (1 + a)\p\ = \ap — Pg\ is valid. Therefore 
we find 



l^^P-Pgl > \/l + 



a- 



(12.86) 



which is valid for all p and pg (not only in the centre-of-mass reference frame) . 
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We can also deduce the collision time tc from Eq. (12.83) as the time it takes for the 



upper integration boundary to change from ±1 to =Fl- Using also relation (12.85) we 
find 



l + a \p\ 



(12.87) 



or in a general reference frame 

tc 



m 



\ap-pg\ 



2JW| + 



V - v„ 



{12.1 



This is a very intuitive and expected result in that it is the time non-interacting wave 
packets need to cross each other. 
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Chapter 13 



State transformation from a 
measurement interpretation 



In this chapter we will derive a Kraus operator representation of the collisional trans- 
formation of the Brownian particle state. In doing so, we will work out measurements 
which the gas particle performs on the Brownian one. For that, we take a well prepared 
gas particle state \xg,pg) at t = 0, and let the gas particle collide with the Brownian one 
which is assumed to be in an unknown state. Then we ask ourselves what information we 
can one obtain about the initial Brownian particle state, by performing a measurement 
on the gas particle after the collision. The answer depends, of course, on the type of 
measurement performed on the gas particle, and we have to find some kind of optimal 
measurement. 

The outline is much the same as usual in indirect quantum measurements (see sec- 
tion 2.4): One wants to get information about the system's state at time t = 0. For that 



one first prepares some known state for the meter (here the gas particle). Then, the 
meter is made to interact with the system (the Brownian particle), and finally, at some 
time tf after the interaction one performs a measurement on the meter. If the initial 
state of the meter and the interaction are chosen appropriately, then the result of the 
meter measurement gives the sought information about the system. 

Several standard textbooks discuss measurements of this outline (see e.g. [5J and ^). 
However, the examples are usually limited to two simplifications: First, the interaction 
of meter and system is taken to be instantaneous, meaning that it is switched on and 
off quickly enough such that the free evolution of meter and system can be neglected 
during the interaction time. This can certainly not be true in our case, as actually the 
free evolution of gas and Brownian particle "switches" on the interaction by bringing 
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the particles together. Second, the measurement on the meter is usually of projection 
type (or a smeared out version), whereas the measurement we will apply on the gas 
particle is not of this type, but a quantum-limited position- momentum measurement. 

We are not aware of any examples of the most general type (none of the above 
simplifications) of indirect measurements. Because of this we feel that the study of 
collisional measurements is a very interesting one by itself. But most importantly, we 
will employ results of this chapter later when setting up a master equation for the 
Brownian particle. 

13.1 Effect operators 

Although the state of the Brownian particle prior to the collision is not known, it is 
sufficient to consider Brownian particle states of the form \x,p) to uniquely determine 
the effect operators, because is a (over-)complete set of states. We recall from 

the preceding chapter that if mW"^ = nigWg, the scattering process can be described by 

scattering 

\Xg,Pg)^\x,p) Ug{tf)\Xg,Pg)0U{tf)\x,p) (13.1) 

where Ug{t) and U{t) are the evolution operator for the free gas particle and the free 
Brownian particle, respectively, and tf is any time after the collision took place. The mo- 



mentums and positions after the collision relate with the initial values as in Eq. ( 12.73 ) 



(12.76). Since at time tf the two particles are in a product state, a measurement on the 
gas particle can not change the state of the Brownian one. 

We want to perform a measurement on the gas particle, which is optimal (see efficient 



measurements in section 



2.4) in obtaining information about the Brownian particle. The 



right hand side of Eq. (13.1) suggests to use Ug{tf) \xg,pg) as measurement basis on the 



gas particle. The corresponding effect operators are 

^gi^a^Pg) = i^Ug{tf)\xg,pg){ig,pg\Ul{tf), (l3.2) 

where measurement outcomes are denoted with a tilde. Keep in mind that the mea- 
surement on the gas particle is done after the interaction took place, and therefore the 



gas particle effect operators Eq. (13.2) are applied at time tf. As the process should 
give information about the Brownian particle's state before the collision, the Brownian 
particle's effect operators are applied at time t = 0. 



After we have chosen the effect operators Eq. (13.2) acting on the gas particle's 



Hilbert space, we can determine the corresponding effect operators ■n{x,p) for the Brow- 
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nian particle by comparing the respective probabihties (see Eq. (3.16) of [S]) 

{x,p\ 'k{x,p) \x,p)dx dp 

I U^itf) ^ {X,p\ U^tf)\ \TTg{Xg,Pg)0l\ \U{tf) \X,P)0 Ugitf) I 



dXg dpg 



1 



2'Kh 

dXg dpg 



2-nh 



Xg,Pg\Xg,Pgl dX g dp g 

exp 



[2x - {l-a)xg - {l+a)xgf [2ap - {l-a)pg - {l+a)pgfW'^ 



2{l+afW'i 



2{l+afh'^ 



(13.3) 



where Eq. (12.73) - (12.76) were used in the last equality. This equation can be solved 
for all {x,p) with the ansatz 



n{x,p) 
w{x,p) 



2TTh 



dx' dp w{x' ,p ) \x + x',p +p'y^x + x\p +p'\ 



2a 



■Kh{l-aY 



exp 



2a 



(1 -a)2 



X 



+ 



(13.4) 



with W"^ = aW^. Substituting Eq. (13.4) into the left hand side of Eq. (13.3), we find 
after performing two Gaussian integrals, that 

dx dp 4a 



{x,p\ tt{x,p) \x,p)dxdp 



27Th (l+a)2 



exp 



2(x — X 
PF|(l+a)2 



2a^W^{p-pf 



h'^n+a)^ 



This has to be compared to the right hand side of Eq. (13.3) to find 

1 — a 1 + a _ 



-Xo + 



Xn 



p 

dx dp 



2 

1 — a 1 + a 
-Pa + -^P9 



2a 
{l + af 
4a 



2a 
dxgdpg. 



.(13.5) 

(13.6) 
(13.7) 
(13.8) 



That is, from the measured values {xg,pg) we can use Eq. (13.6) and (13.7) to infer the 
approximate position x and momentum p of the Brownian particle. The measurement 



on the Brownian particle described by the effect operators Eq. (13.4) is a smeared out 
version of a quantum-limited simultaneous position-momentum measurement. Whether 
the measurement is more precise in position or momentum depends on the choice of Wg 
of the initial gas state. The smearing out results in additional uncertainty on top of 
the Heisenberg uncertainty limit, and can be reduced by choosing a gas particle with a 
similar mass as the Brownian particle. 



In the limit of a light gas particle we see from Eq. (13.5) that the position precision 
of the indirect measurement is half the position variance Wg of the initial meter state 
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\xg,pg)wg- On the other hand, the momentum precision is reduced by a factor a com- 
pared to half the momentum variance h/Wg. The reason is that in this hmit the gas 
particle measures twice the velocity of the Brownian particl^ and the precision is its 
own velocity variance. As momentum is velocity times mass, the momentum precision 
is reduced by the ratio mg/m. 

For a = 1 we find w{x,p) = 6{x,p) and the effect operators take the form of a 
quantum-limited position-momentum measurement. Of course, if a = 1 then a collision 
just swaps the states of the two colliding particles, and therefore a measurement on 
the gas particle after the collision is equivalent to a measurement on the Brownian one 
before the collision. Hence, a quantum-limited indirect measurement is no surprise if 
both particles have the same mass. Furthermore, because both particles swap their 
state, the Brownian particle carries no information about its pre collisional state and a 
subsequent measurement can not reveal any more information, therefore assuring that 
the Heisenberg uncertainty principle is not violated. 



Note that the effect operators Eq. (13.4) were already used (but not derived from 
first principals) in [60] to phenomenologically describe collisional measurements. But 
the authors used free parameters to specify the precision of momentum and position 
measurements, and were not able to give a relation to the mass and the initial state of 
the gas particle. 

It is often stated in the literature (e.g. [58]), that a colliding particle measures the 
position of a tracer particle, but not the momentum. In this section we showed clearly 
that this statement is not generally true. Instead, we find a beautiful symmetry between 
measurements of momentum and position, despite an interaction potential which solely 
depends on position variables. 



13.2 Kraus operators 

In this section we are concerned about the change of the Brownian particle's state vector 
due to a measurement. Without loss of generality we assume that the Kraus operators 
acting on the gas particle's Hilbert space are the square root of the effect operators 
(using more general Kraus operators give the same result for the Kraus operators acting 
on the Brownian partcle) 

Ag{Xg,Pg) = -^l=U g {t f) \x g , P g^X g , P g\ {t f) . (l3.9) 

^This is because in the limit a — >■ the gas particle velocity after the collision is Vg = —Vg + 2v. A 
measurement of Vg with some uncertainty Avg then results in half of this uncertainty for the velocity v. 
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Therefore, a measurement with result {xg,pg) changes the two particle state Eq. (13.1) 
to (omitting normalization) 
1 



{Xg,Pg\ Xg,Pg)Ug{tf) \X g , f) g)U {t f) \x , p) 



(13.10) 



The Kraus operators acting on the Brownian particle depend on the measurement 
result (x,p) as well as on the state of the gas particle which was used for the mea- 
surement, i.e. {xg,pg). Measurement theory [8j states that the Kraus operators can be 
written as 

A{xg,pg;x,p) = V{xg,pg;x,p)^/7r{x,p), (13.11) 

with some unitary operator V{xg,pg;x,p) which is yet to be specified. Together with 
Eq. ( |13.10 ) we see that this operator has to satisfy 



uHtf)V{xg,pg;x,p)^TT{x,p) \x,p) 



1 



{Xg,Pg\ Xg,Pg)\x,p). 



(13.12) 



The phase of the scalar product depends on the state of the Brownian particle and 
because the Kraus operators will act on general states (i.e. on superposition states) we 
have to keep track of it. We find for this phase 

2a{px — px) + a(l — a)xg{p — p) + {1 — a)pg{x — x) 



arg(Xg,Pg| Xg,Pg) 



(13.13) 



where the relations Eq. (12.73), ( |12.74 ), (13.6), and (13.7) were used. We further need 



\/t!:{x,p) \x,p) 



cexp I 



,a{px — px) 
h{l + a) 



1 — a 2a 1 

X + z X, 



1 + a 1 + a'l + a 1 + a 



P + ^^—p), (13.14) 



where c is a real number. The derivation of Eq. (13.14) is a bit longer and can be found 



in the appendix of this chapter. Next we substitute Eq. (12.75) and (12.76) into the 
!;ht hand side of Eq. 

U''{tf)V{xg,pg;x,p) = exp(^- 



right hand side of Eq. (13.12), and Eq. (13.14) into the left hand side, to find 



. aXgP — XPg 

'V' h{l + a) 



D 



( Mxg-x)2{pg-ap) \ 
\ 1+a 1+a ) 



where Eq. (13.13) was used, as well as the Glauber displacement operator D{x,p) for 
he Brownian pj 

A{xg,pg;x,p) 



the Brownian particle. Substituting into Eq. (13.11) we finally get the Kraus operators 

x) 2{pg 



expii^lf^-^ ) u{tf)D 



h{l + a) 



(2a{a 



ap] 



V 1 + Oi 



1 + a 



\/7r(x,p). 
(13.16) 



As required, the Kraus operators do depend on the measurement outcome {x,p) 
and the state of the measurement apparatus, i.e. {xg,pg), but not on the state of the 
Brownian particle before the measurement. 
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The displacement in position seems to be odd, but that can be resolved by considering 
the finite time tf the measurement (or the collision of the wave packets) takes, as well 
as the free Brownian particle evolution operator U{tf). In fact, if one assumes the 
collision takes place instantaneously at the time when the centers of the wave functions 
of the gas and the Brownian particle are equal, i.e. tf = and x = Xg, then the 



displacement operator Eq. (13.15) just shifts the state Eq. (13.14) back to the original 
position X . Therefore the non-zero position value in the displacement operator is needed 
to counteract the position movement due to the effect operator ^y ■fr{x,p). 

One might wonder why the effect operator decomposes with only one Kraus oper- 
ator rather than the more general decomposition tt = Yli^l^i- The reason for this is 
twofold. First, we used a gas state which is initially pure, thus not introducing classical 
uncertainties. Second, we used an optimal measurement on the gas particle in the sense 
that we measured the basis in which correlations between the two particles are most 
pronounced. 

Now we are in the position to write down the state of the Brownian particle at time 
tf after a collision with a gas particle, which initially is in the state \xg,pg). Because in a 
real collision we do not actually perform any measurement on the gas particle, we need 
to use the general formula for a non-readout measurement (not measuring the meter 
accounts for not reading out the measurement result) 

,M = //c.^dpi,.„p,.,p,„0)it,.,,p,;*,p). ,13.7) 



A nice and important property of this transformation is that it is valid for general states 
of the Brownian particle, and not just for Gaussian ones considered so far. 

Later on we will need the Kraus operators without the free evolution part. Further- 



more, as Eq. (13.17) shows, we do not have to take care of phase factors which only 



result in a global phase. Therefore we define 



- f 2a 2 \ / 

B{xg,pg;x,p) = ^-p— (xg - x), j-p^(p<; - ap) 1 V7r(5,p) 



1 -f{l-a)x + 2axg {l-a)p + 2pg 



^\ T — ^ ^ / VaD{x,p). (13.18) 



where = denotes equal up to a phase, and 



a = 27r^7i-(0,0) 

dx' dp' w{x' ,p) \x' ,p')(x' ,p\ (13.19) 
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is a trace one operator with vanishing first moments of position and momentum. Then, 



Eq. (13.17) reads 



p{tf) = U{tf) 



JJ dx dp B{xg,pg- x,p)p{0)&{xg,pg; x,p) 



Uitj) 



(13.20) 



Appendix 



Here we will derive the action of the square root of the effect operators on the Gaussian 
states. This can be done either in position representation, resulting in a very messy 
calculation, or, by a much neater algebraic calculation, which we choose here. To do so, 
we introduce some harmonic oscillator algebra along the lines of Barnett and Cresser |6Uj . 
They defined annihilation, creation, and number operators via 



it 



h = a^a 



X . Wp 
X . Wp 



V2W 



+ 



V2h 

W^f 1 



2^2 2h? 
and the eigenstates |n), n = 0, 1, 2, . . . of n. 



(13.21) 
(13.22) 
(13.23) 



n \n) 



n \n). 



(13.24) 



They then found that (see Eq. (10.11) for the relation of a to the effect operator) 



l-e-^)e"^" 



n 



1 + n ^ ( 1 

n=0 ^ 



+ n 



\n){n\ 



(13.25) 
(13.26) 



with A = ln(l + 1/n). From the first line we deduce 

Va = Vl - e-^e-^"/2. 



(13.27) 



Also useful is that the Gaussian pure states are eigenstates of the creation operator 

fx , .Wp\ 

(13.28) 



a\x,p) 



w 



\V2W V2h 

a{x,p) \x,p)w, 
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which is easily shown, e.g. in position representation. From tliis, we can derive the 
fohowing formulas 



{n\ x,p) 



1 / \ a{x,p) 

{n — l\a\x,p) = i^—{n — l\x,p) 

n \ n 



-\h/2 



x,p 



n=0 

exp 



|a(x,p)p 



-An/2 



n){n\ x,p) 



(13.29) 



We can now use these formulas to derive the action of the square root of the effect 
operators on the Gaussian states, which will turn out to be necessary to determine the 
Kraus operators corresponding to each effect operator. To this end we use ae"'^"/^ = 
g-A/2g-An/2^ to show that e^"^"/^ |a;,p)is an eigenstate of a with eigenvalue e~'^/^a(x,p), 
and therefore e~^"/^ \x,p)(x le^'^^'^x, e^^^'^p). The proportionality constant is obtained 



from Eq. (13.29) and with Eq. (13.27) we arrive at 

|a(x,p)p 



a \x,p) 



' exp 



1 



(13.30) 



By using Eq. (10.12) it is now straightforward to find the action of the square root of 
the effect operator on a Gaussian state 



^/7r{x,p) \x,p) 



V27rh 



X exp 



exp 



—I 

2h 



px — px — e ^^'^\p{x — x) — x{p — p)] 
/^W{p — p)^'^ 



W 



-r^l'^ix - x) + er^l'^ij) - p) + p^. (13.31) 
To use this result for the measurement a gas particle performs on the Brownian particle. 



we compare Eq. (13.4) in chapter 13 with Eq. (10.8) to relate 

, 2 



-A 



n 



n + 1 



1 — a 
1 + a 



(13.32) 



where a = mg/m is the ratio of the masses of the two particles. Then we find 



\/7r(x,p) \x,p) 



2Ja 



-a 



exp 



ia(xp — xp) 
[ (1 — a)n 



; — x\^ / W{p — p) 



(l + a)2 

2ax + (1 — a)x 2ap + (1 — a)p 



1 + a 



1 + a 



(13.33) 



Chapter 14 

Rate Operator 



Since we know the effect of a collision if the gas particle is in a Gaussian state \xg^pg\ we 
now have to ask what is the rate of such collisions. As in classical collisions, the rate will 
depend on the state of the Brownian particle. For example, the Brownian particle will 
never collide with a gas particle moving with the same velocity. Because of this state 
dependence, the quantum mechanical collision rate will have to be operator valued. 

In the first section we will discuss the probability of finding a gas particle in a state 
\xg,pg). In the next section we show, for a given gas particle state, under what conditions 
a collision with the Brownian particle will occur. 



14.1 Gas density operator 

Before we derive the probability of finding a gas particle in a state \xg,pg), we start 
with some important remarks about the density operator of a gas particle of an ideal 
Boltzmann gas in a thermal state. In particular we wish to point out that there are 
many different convex decompositions 



Pg = j daP{a)\^g{a)li,g{a)\ (14.1) 

of such a thermal density operator. Here a is a finite set of parameters, \ipg{a)) is a set 
of normalized gas particle states (not necessarily orthogonal) spanning the Hilbert space 
of a gas particle, and P{a) > is the probability density of the gas particle to be in the 
state \'ipg{a)). This freedom has some important consequences, one of which is that the 
notion of a collision time tc, which is roughly the time required for a wave packet of the 
tracer particle to cross the width of the wave packet of a gas particle (a more precise 



definition is given in section 12.4), is not to be viewed as a physical parameter whose 



value is intrinsic to the gas, but rather a parameter whose value can be chosen to suit 
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the needs of the analysis. But already here it is clear that momentum eigenstates which 



are often used in Eq. (14.1) are not appropriate to define a collision time. 

The thermal state of an ideal gas particle decomposed in momentum eigenstates is 
usually given in terms of the Maxwell-Boltzmann distribution urip) 

2TTh 

Pa = 




;i4.2) 



\p9lp9 



where T is the temperature of the environment and L is a normalization length. Horn- 
berger and Sipe [38] showed how this density operator can also be decomposed in pro- 
jectors to minimum uncertainty wave packets \xg,pg). For this purpose we need use the 
equality 



exp 



-pI \ _ / T 



2mgkBT J V 2-KmgkBTT J \2mgkBT 



[dpgeJ^^]e.J -^P^-P^' ) (14.3) 
J ^ \2makBT V 2m,kBT J 



which is true if T + T = T as is easily shown by performing the integration. Then we 
split the second term and include the identity operator 1^ = J^^dxg \xgXxg\ in position 
representation to find 



jTTh? f dXg f f -^p p ^ ^ f -^p p )2 ^ 



(14.4) 



'2TTm„kBT 



where 



is the Maxwell-Boltzmann distribution with modified temperature T. Finally we have 
to show that ^ / ^^^'^77. exp ( , ') Ixn) are the desired Gaussian states by going to 
the position representation 



n^nksT J V ^rUgkBT 



^gr 

2TTh^ Y [ {-{V-Pgf^ 



, „, dp exp [ /"^ \{x'\pIp\ 
mgkBT J J V -imgkBT J 



1 



1 A' [^,,:,J-j£-t£^,M^ 



{2TT)'^h?mgkBT J J V ^TagkBT 
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dp exp 



P-Pg ^x'g-Xg'^^ 



{2'KYn?mgkBT J J \4mgkBT " 2h 

Pgi ^'g-Xg) _ mgkBT{x'g-Xgf 



2 



2mgkBTY ( Pgi^'g-Xg) mgkBT{x'g - Xg) 

/ {X'g-Xgf \ 

exp (14.6) 



9 



where Wn = , ^ _i is defined. Up to a global phase, that is exactly the position 

y/2mgkBT 

representation of \xg,pg), and we arrive at the desired decomposition of the thermal 
density operator into Gaussian states 

P9 = J -jf-J (^Pg f^f(P9)\^9^P9)i^9^Pg\ ■ (14.7) 

We see that part of the thermal energy is associated with the wave packets themselves, 
while the rest is in the motion of the centres of the wave packets via fij,{pg). As our 
collisional calculation is only valid if the mean momentum of the gas particle is large 



compared to its momentum uncertainty, we will use Eq. (14.7) with 

T>T, f^T (14.8) 

in the following. 

14.2 Definition of the rate operator 

The probability for the Brownian particle to collide with a gas particle state \xg,pg), 
during a time interval 5, should be determined by the expectation value of a rate operator 
Rs{xg,pg) times the time interval 6. Because we know the state of the thermal gas 
particles, and we are looking for an operator acting on the Brownian particle, it is 
advantageous to rephrase the question: // at time t we are given a gas particle in the 
state \xg,pg), what is the probability that it collides with the Brownian particle in the 
interval (0, 6) ? 

In classical physics, for a given gas particle with position Xg and momentum pg, there 
will be a collision with the Brownian particle exactly if its position and momentum are 
within the phase space region 

Pg/mg-p/m J 



Ss{Xg,Pg) = < {X,p) 



136 



CHAPTER 14. RATE OPERATOR 



That is, the particles either coUide or do not colhde. There are no partial collisions 
because classical hard core collisions are instantaneous. This phase space region is 



shown in figure 14.1 (a) 



(a) 



m 



(b) 



S5{Xg, Pg) \^ 




m T / 



Wm 



X 



Figure 14.1: (a): If a classical particle is in the phase space region Ss{xg,Pg), then it will 
collide with a classical gas particle with position Xg and velocity Pg during a time interval 6. (b): 
Ts{xg,pg) acts like a projection operator well within Ss{xg,Pg), but not near the edges of this 
phase space region which is indicated by the blurry area. If S is sufBciently large, the edges can 
be neglected in comparison to the interior. 



In a quantum collision we have to be more careful. Because the collision of wave 
packets is not instantaneous, we have to account for the possibility that at time 5 a 
collision is not complete. However, if we choose 6 sufficiently large 



(5 > 



2^(1 + a)Wg 



(14.10) 



where tc is the collision time Eq. (12.88), then the probability of a complete collision 
will outweight the probability of an incomplete one, and we can neglect the latter. Note 
that relation (14.10) can not be fulfilled for v close to Vg, but for sufficiently large 5 



(and/or for high temperatures) relation (14.10) can be achieved for most gas particle 
velocities. We can then use guidance from the above classical statements also for a 
quantum collision, at least in an approximate manner. 

The aim now is to restrict the collision to Brownian particle states which are located 
within the phase space region Ss{xg,pg). To this end we define the quantum mechanical 
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"phase space projection operator" 



1 



2TTh 



dx dp \x,p){x,p\ 



(14.11) 



ISs{Xg,Pg) 

Strictly speaking, phase space projections can not exist in quantum mechanics because 
of the Heisenberg uncertainty principle. But if the Brownian particle is in a Gaussian 
state \x,p), with values of x and p well within Ss{xg,pg), then 



^s(.Xg,pg) \x,p)^ \x,p) 



(14.12) 



and there will certainly be a collision. On the other hand, if {x,p) is well outside of 
Ss{xg,pg), then 

ts{xg,pg)\x,p)^0 (14.13) 



and there will be no collision. If 6 obeys relation (14.10), then the weight of the values 
X and p which are near the boundary of Ss{xg,pg) becomes small, and 

^2/ 



^si^g^Pg) ~ '^s{xg,Pg) 



(14.14) 



is indeed an approximate projector operator. The action of this operator can be visual- 
ized in Fig. |14.1| (b). 

To get the rate operator we also have to consider the probability of having a gas 
particle in the state \xg,pg). If the gas is an ideal one and in a thermal state, we saw in 
the previous section that this probability is 



n 



gf^fiPg) 



2Trm„kBT 



exp 



-pI 



2mgkBT 



(14.15) 



^2 



where Ug is the particle density of the gas and T = T 
divide the operator by the time interval 6 to get a rate, and arrive at 

,'^s{xg,pg) 



2mgkBwr Furthermore, we 



R5{Xg,Pg) 



UgH^iPg)- 

n„ 



6 



: exp 



6\ 2TTmqkBT 



-pI 



2m„kBT 



SsiXg,Pg) 



dx dp 

2TTh 



x^p){x,p\ . (14.16) 



Note that the rate operator depends on the time interval 5, which is not uniquely fixed 



by physical parameters, but can be chosen freely as long as it satisfies Eq. (14.10). 
However, the trace norm Tv[Rs{xg,pg)] of the rate operator does not depend on 5. Also 
we will see that on a coarse grained time scale the resulting master equation does not 
depend on this parameter. 
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14.2.1 Integrated rate operators 

We also define a rate operator corresponding to any collision with a Gaussian gas state 
with mean momentum pg 



R{Pg 



n, 



dXg Rs{Xg,Pg) 



(14.17) 



dxdp 
ax„— — — \x,p}{x,p\ 



rriQ m 



n„ 



exp 



-Pi 



2TTmgkBT X^nigkBT^ 
exp 



n„ 



y^lTTnigkBT \2mgkBT 



dx dp 

2TTh 



Pa P 



P^ 



p_ 
m 



\x,p){x,p\ 



(14.18) 



where Eq. (14.8) as well as Eq. (10.27) was used. The integration over Xg is easily 



performed because the integrand does not depend on this value (see Eq. (16.82) for 
more detail). Furthermore we define the total collision rate 



R 



dpg R{Pg 



n, 



g y 2kBT rr dx dp 
.JwnTg JJ 2Trh 



ngy/2kBT 



TTm„ 



exp 




2kBTm? 



where we used again Eq. (14.8) and Eq. (10.27). The integration over is also straight- 
forward if the cases ^ — ^ ^ ^ ^'^^ considered separately. Note that neither R{pg) nor 
R depend on 5. In fact, these collision rates are precisely the ones found in the classical 
calculations in section [9l 

Also useful will be the total collision rate in the limit of slow Brownian velocities 



\v\ <C \Va 



nin 



R 



ngyj2kBT 



1 + 



as well as 



U\t)R{Pg)U{t) 

U\t)RU{t) 



nig p 
2/cBr^ 



R{Pg) 
R. 



(14.21) 



(14.22) 
(14.23) 



Chapter 15 

Master Equation 



Now, as we have developed formulas for the change of the state due to a collision 



Eq. (13.17), as well as for the rate of a collision Eq. (14.16), we are in the position to 



derive a master equation describing the Brownian particle experiencing random collisions 
with gas particles of an ideal gas in a thermal state. 

15.1 Discrete master equation 

Let us first consider the change of the Brownian particle state during some finite time 
6. We assume that 6 is small compared to the collision frequency 



dTT[Rp{0)] <€. 1, 



:i5.1) 



such that we can neglect the possibility of two collisions in (0, 6). If there is no collision 
in this time period, then the density operator undergoes the unitary transformation 

p{6) = U{6)pm\6), (15.2) 

whereas if there is a collision with a gas particle in the state \xg,pg) the transformation 



is according to Eq. (13.20|) 

p{5) = U{6) 



jj dx dp Bixg,pg; x,p)p{0)B\xg,pg; x,p) 



U\6). (15.3) 



To find the actual density operator at time t = 6, we have to account for collisions 
with all possible gas particle states with the correct probabilities. Quantum trajectory 
theory [76] suggests to multiply each Kraus operator by the square root of the probability 
operator 



Pi6) = U{6) 



l-SRp{0)yi-SR + 




dxg dpg dx dp 
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I3{Xg,Pg]X,p)j6Rs{Xg,Pg)p{0)j6Rs{Xg,Pg)B\xg,Pg-,X,p) 



uHd), (15.4) 



where 1 — 5R gives the probabihty of having no colhsion. The probability operator 
5Rs{xg,pg) ensures that only the part of p(0) which is within the "range" of \xg,pg) 



undergoes the transition Eq. (15.3), weighted by the thermal distribution fxxiPg)- The 



'range" is defined by Eq. (14.9). The integration in Eq. (15.4) is over all possible gas 



particle states {xg,pg), as well as over all possible "measurement outcomes" {x,p). 

It seems advantageous to use the interaction picture Pint{t) = {t)p{t)U{t) for the 
following 



Pint{S) = yl-6Rpint{0)yl-dR + 6 JJJJ dxgdpgdxdp 

X I3{Xg,Pg]X,p)\j Rs{Xg,Pg)pint{^)\j Rs{Xg,Pg)B\Xg,Pg\X,p). (l5.5) 



Because 6 is small compared to the inverse collision rate, we can expand 
l-SR/2 to find 

P^nt{S)-pUO) ^ -hRp^,^,(0) + pM0)R]+ [[[[dxgdpgdxdp 



X I3{Xg,Pg]X,p)\j Rs{Xg,Pg)pint{^)\j R5{Xg,Pg)I3\xg,Pg-,X,p) 

= Csipintm (15.6) 

which is the discrete master equation in the interaction picture. Let us have a quick 
look at this equation. The operators B and B\ integrated over x and p, transform 
— )• \x,p). The rate operator ensures that only states undergo this transformation 
which are in the "range" of a gas particle within the state during the time interval 

(0, 5). Furthermore, this operator includes the probability of finding a gas particle in the 
respective state. Finally, the integration over Xg and pg averages over all possible gas 
states. The anti commutator arises from the non-collision part, and ensures a constant 
trace of the density operator. 

15.2 Continuous master equation 



To derive a continuous master equation from Eq. ( 15.6 ), one would wish to take the limit 
5—7-0. Here we encounter a problem. From the discussion in section 14.2 we know 
that we can not take this limit in the rate operator Rs{xg,pg), because no complete 
collisions can occur in an infinitesimal time interval, and in the derivation of the rate 
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operator we required the ability to neglect incomplete collision. In fact, as can be seen 



from Eq. (14.11), taking this limit in the rate operator would account for performing 
quantum-limited position- momentum measurements on the Brownian particle, and these 
are by no means performed by a colliding gas particle. Therefore 6 strictly has to satisfy 
relation ( |l4.10[ )p] 

If we now perform the continuity transition 

Pint{dt) - Pint{t) Pinti^) - Pint{0) n 7\ 

dt ^ 6 ' ^^^-^^ 

we arrive at a master equation, which accounts for the free evolution in an exact man- 
ner (because we take the continuity transition in the interaction picture), and for the 
collisions in first order in 5. This transition essentially takes all the possible collisions 
which might occur in (0, 5) and performs the corresponding state transformations (with 
probability reduced by dt/6) at time t = rather than spread out during (0,5). It is 



this step where position diffusion is introduced. (See also figure 12.5 ) If we look 
at the real evolution of a Brownian particle (initial state \x,p)) colliding with a gas par- 
ticle (initial state \xg,pg)), we find that the position of the Brownian wave packet evolves 
continuously from x to x + 5p/m during the time 6, with x and p given by Eq. (12.75) 



and (12.76), as is clearly seen in Fig 12.4 However, in the continuous master equation, 
which is derived from the discrete master equation, there is a jump from x to x occurring 
at time t = 0, and therefore position jumps are introduced. In fact, it is well known 
that any attempt to derive a Markovian master equation (and with our set up we can 
only derive this type of master equation) describing QBM without introducing position 
jumps must be ill-fated |5T], as it can not be brought into the required Lindblad form. 



If we were to use the continuity transition (15.7), the position jumps would not be 
uniform. In fact we would find for the position expectation value (for simplification we 
take |f I ^ ^k^T jrug., more about the derivation of expectation values in chapter 16) 

4a\/¥ern„ 




(15.8) 

where the latter part is due to the non-uniform position jumps. This part is propor- 



tional to the number of collisions and therefore small and hence Eq. (15.8) would be an 
unorthodox, but acceptable result of a master equation. However, it is straightforward 
to avoid such unphysical contributions linear in b by noting that for any smooth function 



^One could say that the rate operator in the master equation is supposed to measure the coUision 
rate without disturbing the state much. This can only be achieved approximately if the rate operator 
performs a very imprecise measurement, which in turn is achieved by a large 5. 
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f{S)~f{0) 

s 



f{5/2+dt)-f{5/2) 



h 0{6 ) is valid. Therefore we will perform the continuity 

transition to second order in 5 



Pint (f + dt) 



Pint {2] 



dt 



PintjS) - PintjO) 
6 

Cs[pint{0)]- 



:i5.9) 



This transition essentially takes all the possible collisions which might occm' in (0, 6) and 
performs the corresponding state transformations (with probability reduced by dt/6) at 
time t = 5/2. As this time lies symmetrically in the above time interval, the (non- 
physical) position jumps will be uniform and not show up in the position expectation 
value. But they will show up in second order in 6 in the position variance, which is 
unavoidable if one is seeking for a Markovian description of QBM. 

After we have done the continuity transition, we can return Schrodinger picture by 
using 



Pint{^+dt) - Pintii) 



dt 



[/t(|; 



C/t(dt)p(f+dt) U{dt) 



Ml) , 1 

dt h 



dt 



H,p 



HI) 



(15.10) 



In the second step we expanded the free evolution operator U{dt) = 1 — IHdt. Finally 



we substitute the discrete master equation Eq. (15.6) in the right hand side of Eq. (15.9), 



noting that Pint(O) = p(0), and multiply with the free particle evolution operator from 
the left and right to find 

Ml) _ 



dt 



-^[H,p{l)]-U{l)U[Rpit) + Pm] 




dxg dpg dx dp 



B{Xg,Pg]X,p)jRs{Xg,Pg)p{0)jRs{Xg,Pg)B\Xg,Pg-,X,p) > U^' {^) . (l5.1l) 



As last step we substitute p{0) ^ {^) p{^)U'< {^) into the incoherent part of Eq. ( |l5.1l[ ) , 
which is a justified as long as two collisions during 6 are negligible. Finally we use ho- 
mogeneity in time to replace (1) by t to arrive at our final master equation (also using 



Eq. (14.23)): 



dp{t) 



dt 



[H,p{t)]--[Rp{t) + p{t)R] + 




dxgdpgdxdp 



U{l)B{Xg,Pg-X,p)^fl^i^^^)U\l)p{t)U{^ 

(15.12) 
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This equation is one of the main results of our work. 

Note that the rate operator stih depends on 6. However, as long as 5 is chosen 
appropriately and all other limits and approximations are satisfied (see next section), 



Eq. (15.12) should be a good description of the actual physical situation. The influence 
of the precise value of 6 on the evolution of the density operator is of second order in 
6 and negligible on a coarse grained time scale. A demonstration of this point can be 
found in the evolution of the expectation values of position and momentum as well as 



their variances in chapter 16 



If we had used (15.7) instead of the more correct (15.9), we would have found the 
same master equation except without free particle evolution operators. In this case the 
influence of the precise value of S would be of first order. As discussed earlier, that is not 



strictly forbidden, but using this slightly more complex master equation (15.12) results 
in a neater form of equations of motion for the lowest order moments. 

15.3 Review of approximations 

In the preceding chapters we used several approximations and limits. At this point 
we want to summarize them, and to discuss under what physical conditions they can 
simultaneously be fulfilled. 



The collisional calculation in chapter 12 is valid exactly if 



\xg -x\ > \/l + aWg (15.13) 
\vg-v\ > VTTa— |— . (15.14) 

mgWg 

The first relation states that there is negligible initial overlap of the two wave functions, 
and the second one states that their mean relative velocity has to be large compared to 
the uncertainty of their relative velocity (small overlap of the initial wave functions in 
momentum space). 

In the definition of the rate operator we required 

5 > VlTa , , . (15.15) 
\Vg-v\ 

to ensure that most collisions in 5 are complete. 

Finally, the derivation of the master equation required the low-density limit 



^/wmg 



5Tr[i?p(t)] w ^ 9^ ^ ^ (;L5_;Lg) 
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to ensure that we have to consider at most one coUision during S. 

Of course the above relations cannot be fulfihed for ah Xg and pg occurring in a 
thermal gas. For our description of QBM to be vahd, they should be fulfilled for most 
gas momenta pg, and for a given gas momentum they should be fulfilled for most gas 
positions Xg which can reach the Brownian particle's wave function in the time interval 
5. 

Using that a gas particle can reach the Brownian one only if |xg — a;| < 6\vg 



v\, 



relation (15.15|) follows directly from relation (|15.13|). For a typical collision we can use 

\Vg -V\- 



!^s21 and relation (15.14) becomes 

mg ' * ' 

\/l + ah 

^TTlgkBT' 



Wg » 



(15.17) 



which gives a lower bound to the parameter Wg. Substituting this result into rela- 



tion (15.15) gives a lower bound for the coarse grained time scale in terms of the tem- 



perature 



6 > 



[l + a)h 



(15.18) 



which can be interpreted as the high-temperature limit. Substituting into Eq. (15.15), 
we further find 



Wg « 



TT 



1 



2(1 + a) ug 



(15.19) 



Putting relations Eq. ( |15.16[ ) and Eq. ( |15.18[ ), as well as Eq. ( |15.17[ ) and Eq. ( |15.19[ ) 
together, we find see that the parameters 5 and Wg are bound from above and below by 
the respective relations (disregarding constants) 

n 



/m„ 



h 



■s/rUgkBT 



Wg « 



ng^/ksT 



n„ 



(15.20) 
(15.21) 



Both of these relations can be simultaneously satisfied in the high-temperature and 
low- density limit 



Ugh 



< 1. 



(15.22) 



^Jm^^k^ 

This is the crucial relation for the validity of our QBM approach. As such, it does 
not contain any freely chosen parameters like Wg or 5, but only depends on physical 
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parameters like the temperature and density of the gas. As expected, high temperature 
can be traded off with low density, but as we approach lower temperatures, the coarse 



graining time scale Eq. (15.18) on which our description is valid becomes larger. 

We close this section with a comparison to classical Brownian motion (see chapter [9]). 
There, the limit of a small momentum transfer per collision is often used together with 



the limit of frequent collisions. It is apparent from Eq. (15.22), that this limiting pro- 
cedure is not possible in a strict treatment of QBM. The underlying reason is that the 
wave functions of the individual gas particles would overlap in this limit. Then, the ideal 
gas can not be described by Boltzmann statistics any more, and has to be considered as 
a quantum gas using the appropriate quantum statistics (Fermi-Dirac for fermions, and 
Bose-Einstein for Bosons). 



CHAPTER 15. MASTER EQUATION 



Chapter 16 

Expectation Values 



In this chapter, we wih use Eq. (15.12) to calculate equations of motions for the lowest 



order moments of position and momentum. The general approach to find an equations 
of motion for an expectation value of an arbitrary time-independant observable O is as 



follows. Recalling the formula of an expectation value of an observable from section 2.3 
we find 



(O) = 

_ d(0) 
dt 

Now we can substitute the master equation 



Tr 



Tr 



Op{t) 
dp{t) 



o 



(16.1) 



dp(f) 
dt 



Cp{t) into the right hand side. As 



the master equation is linear in p{t) and local in time, it is always possible to use the 



invariance of the trace under cyclic permutations to bring Eq. (16.1) in the from 



d(0) 
dt 



Tr 



£to) p(t) 



where £^ is called the adjoint Liouville operator. 

In the following we will discuss the observables x, 
challenge will be to express, e.g. jC^x^, in terms of these observables to find a closed set 
of equations for the lowest order moments. 



(16.2) 



p, p^, and {x^p}. The 



16.1 Mean position 

In this section we will explain each step in detail, whereas in the following sections we 
will limit ourself to outline the calculation and to refer to this section for the details. 
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Substituting Eq. (15.12) into Eq. (16.1) and using the invariance of the trace under 
cycUc permutations, we find for the expectation value of the position operator 



d{x) 



dt 



Tr 



2hm 



X p 



Tr 



1 



xR + Rx)p 



+ Tr 



dxn dpg dx dp 



X U{l)^/M!il~[^)B\Xg,Pg-,X,p)U^ 

(16.3) 



Therefore we have 



C^x 



2hm 



[p 



a2 



x] - -{xR + Rx) + 




dxg dpg dx dp 



X U{l)^jRs{Xg,Pg)B\Xg,Pg-,X,p)U^{l)m{l)B{Xg,Pg;X,p)^Rs{Xg,Pg^^^ 

(16.4) 

The first term results from the unitary evolution of the Brownian particle. Evaluating 
the commutator, this term is the expected Brownian velocity operator p/m. On the 
second line we start in the centre at the position operator and work outwards. Using 
Eq. ( 10.30 ), we find [U{t),x] = -^U{t) from which we deduce 

p5 



?7t(|)x[/(|) = f + 



2m 



:i6.5) 



Next, we substitute Eq. (13.18) for B{xg,pg;x,p). Furthermore, we apply 



D"^ {x,p)xD{x,p) 



dx X D\x,p) \xy^x\ D{xp) 
dx' x' \x' — x^^x' — x\ 
dx' {x' + x) l^^'X^^'l 

X + X 



:i6.6) 



{x,p)pD{x,p) = p + p 



(16.7) 



to write 



Jj dxdpB''{xg,pg;x,p)U''{^)xU{^)B{xg,pg;x,p) 



dx dp 

2TTh 



b{x,p)vd 



1 — a 2axq 6 

X H X H H I p - 

1+a 1+a 2m' 



1-a , 2pg 

p H — 

1+a 1+a 



2axn 



+ 



1 — a 



-X 



1 + a 1 + a 2m \1 + a 1 + a 



2P, 



+ 



1 — a 



-p 



'db\x,p) 
(16.8) 
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In the second step, Eq. (16.76), Eq. (16.77), and Eq. ( 16.78| ) from appendix A of this 
chapter were used, as weh as Tr{a) = 1 and Tr:{xa) = Tic{pa) = 0. Substituting 



Eq. (16.8) into Eq. (16.4) and using Eq. (14.16) we find 
p 1 



m 



;{XR+RX) + jj dXgdpg^^^^ 

2axg ^ 1^^ + — f + 
1 + a 1 + a 2m \1 + a 1 + a 



^ts{Xg,Pg)U\l). 

(16.9) 



With Eq. ( |10.22[ ) we write 

^5{Xg,Pg)x 



II 



dx dp 

2'Kh 



xts{Xg,Pg) \X,P\X,P\ 



dx dp 



tng m 



X \x,pXx,p\ 



(16.10) 



where we also used relations (14.12) and (14.13). The same formula applies for xf5(xg,Pg) 
and \ /Ts{xg,Pg)xyTs{xg,Pg) and therefore we can use 



(16.11) 



XTs{Xg,Pg) Jrs{Xg,Pg)xJTs{Xg,Pg) ^ rs{Xg,Pg)x 



The same property, which is a consequence of Ts{xg,pg) being an approximate projection 
operator, also applies to the momentum operator and we can pull all R and Ts{xg,pg) 
across position and momentum operators 



p_ 

m 



xR+ II dxgdpg — — 



X U{l) 



2axa 1 — a ^ 6 f 2pa 1 — a 

H X H — ^ H p 

1 + a 1 + a 2m \1 + a 1 + a 

Next we use Eq. ( |14.23| ) and Eq. ( |16.5D for 



fs{Xg,Pg)U\l^). (16.12) 



U^{^)xRU{^ = U^a)xU{l)U^{l)RU{l] 



X + —p R, 

2m 



as well as R = ff dxgdpg ^T{Pg)Ts[xg,pg)/5, to find 



2a[xg — x) 6 f pg — ap 



1 + a 



+ 



m\ 1 + a 



(16.13) 



i^s{Xg,Pg)U^{l). 



^+Ua)lldXgdpg //d-^PM..) 



2TTh 



2a{xg — x) 5 f Pg — ap 



1 + a 



+ 



m\ 1 + a 



X |x,pX^,p|C/t(f). 



(16.14) 
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In the second step, we used the definition Eq. ( 14.11 ) of Ts{xg,pg), as well as Eq. ( 16.10 ). 
Upon integration over Xg, which is carried out in Eq. (16.82) and Eq. (16.83) in the 



appendix B of this chapter, we find that all terms in the integrand cancel. Therefore, 
we are left with the expected 



X 



m 



from which we deduce 



d(x) 
dt 



m 



(16.15) 



(16.16) 



That is, the mean position of the Brownian particle changes in time according to its 
mean velocity. 



At this point we wish to mention that if we had used the continuity transition (15.6) 



instead of (15.8), then, upon the integration over x„, the integrand in Eq. (16.14) would 



only vanish in lowest order in 5. The remaining part is of first order, and we would be 
left with a non-physical ~ + '^•^(^) ^'-'^ some function / involving Gauss and 
Error functions. See discussion in section [15.21 for more detail. 



16.2 Mean momentum 



Applying the adjoint Liouville operator Eq. (16.4) on p results in 
1 



p = - - {pR + Rp) + / / / / dxg dpg dx dp 



X U{^)\j R8{Xg,Pg)I3\xg,Pg]X,p)pI3{Xg,Pg]X,p)\j R5{Xg,Pg)U\^] 



Using a calculation analogous to Eq. (16.8), we find 

jj dxdpB\xg,pg;x,p)pB{xg,pg;x,p) 



2P, 



+ 



1 — a 



1 + a 1 + a 



(16.17) 



(16.18) 



which we substitute into Eq. (16.19) 



= _Rp+ JldXgdpgU{l)RsiXg,Pg)(^-:^ + ]-^ 



P 



^t(l) 



/ 



dp,R{P,)^^. 



(16.19) 



In the second step we used the definitions ( 14.17) and ( 14.19) as well as U {t)R{pg)pU\t) = 
R{pg)p (see Eq. (14.22)). Contrary to the previous section, here the contributions from 



p and Pg do not cancel each other. This was to be expected as the gaseous environment 
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has a damping effect on the momentum. Note that at this point there is no dependance 
on 6. 



Now we can use Eq. (|14.18|) to arrive at the final result 

-(Mp)) 



dt 



(16.20) 



where 



MP) 



2m„ 



rin 



1+a y/2TTmgkBT 
Am 



dpg exp 
2pkBT 



2mgkBT 



m 



rrir, 



P_ 
m 



Pg_ 
rua 



—ap^ \ 
2mkBT J 




(16.21) 



involves the error function. For the second line, results from the appendix of chap- 
ter |9] were used. This is the known result for classical Brownian motion, as derived in 
chapter M 



For small velocities p/m <^ y^kBT/nig we can expand fxiP) and write 



d{p) 



Tin 



y/2mgkBT 



dt 1 + a 

Therefore we find for linear friction 

dt 

f 



^{p)_ ^ 2mg {p-^) 
m 3kBT m? 



m„ 



(16.22) 



-/- 



m 



Ang^J2^mgkBT 
(l + a)V^ 



(16.23) 



Another interesting limit is the high velocity limit p/m S> ^Jlzj^TJmg^ as it occurs 
e.g. in a cloud chamber. Then we can neglect the first term in the square brackets of 



Eq. (|16.21|), and extend the integration in the second term to infinity. We find 

kB T{pl\p\) \ 



dt 



2ngmg 
1 + a 



(P \P\) 



+ 



rrin 



(16.24) 



"9 J 

which is the expected quadratic friction and can be explained as follows. In the high ve- 
locity limit, the velocity of the gas particles can be neglected and the momentum transfer 
per collision increases linear with the Brownian velocity. In addition, the collision rate 
increases linear with the velocity, and we end up with quadratic friction. 
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16.3 Mean squared position 



Applying the Liouville operator Eq. (16.4) on results in 

1 



2hm 



2 -21 
X 



x^R + Rx^] + 




dxg dpg dx dp 



X U{l)^Rs{Xg,Pg)B\xg,Pg;X,p)U\l)x^U{l)B{Xg,Pg;X,p)^Rs{^^ 

(16.25) 



We proceed in the same manner as in section (16.1) and find 



// 



dx dp B\xg,pg; x,p)U^j)x^U{^)B{xg,pg;x,p) 
dx dp 



2-Kh 



■D{x,p)y a 



1 — Q_ 2axg ^ 6 / ^ 1 — Q_ ^ 2pg 



X + X + 



+ 



I+q" ' 1+a ' 2m\^^ 1 + 



a 



2a 



(l+a)2 



(—)' 

\2mW J 



l — a)x + 2axg 6 ({\ — a)p + 2p, 



1+a 



+ 



2m 



1+a 



2 



(16.26) 



In the second step, we used formulas from appendix A of this chapter, as well as Tr((T) 

4f^, Tt[^X^x] = 1^2 1^; 



1, Tr(x(T) = Ti{pa) = Tr[a(xp + px)] = 0, Tr[ax2] = VF^ii^, Til^fdx^x] = VF^i^, 



^[^P'] = and Tr[^/apV^p] = 



With Eq. (16.26) and i[p'^,x'^] = 2^{x,p} we find for the adjoint Liouville operator 
applied on the squared position operator 



i^-Rx'+ J"" ^. R 



m 



' 2mW J 



[1 — a)x + 2axg 6 / (1 — ct)p + 2pg 
1 + a 2m \ 1 + a 



+ JJ dXgdpgU{^) R5{Xg,Pg) 
Vt(|) 



{x,p} ^ 2a ^ 
m (1 + a)2 



5h 

2mW 



4a 



+ (1 - a) x + 



5p 



(l + a)2 



[7(2) // dXgdpg Rs{Xg,Pg) 
2' 



Xn + 



2m / V ^ 2mr. 



X + 



6p 
2m 



c/t(f) 

(16.27) 



where we also used Rx^ = U{^)R{x + ^p) C/t(f). Next, we use the definition of 



the rate operator Eq. (14.16) and R = JJ dxg dpg Rs{xg,pg). Then, with Eq. (10.22) to 
Eq. (|10.26|) we find relations like t5{xg,pg)x^ = JJs,(.^^,p^) w (^^ " ^V2) \x,pXx,p\. 



All this substituted into Eq. (16.27), we proceed analog to section 16.1 and use the 
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appendix B of this chapter to perform the integration over Xg 



C)x? 



+ 



\aR 
(l + a)2 



2mW J 



2^2 



+ 



3(1 + q)^ 



p_ 

m 



(16.28) 



Contrary to section 16.1, not all contributions cancel out, but we are left with three 
terms, each representing position diffusion. The first term can be interpreted as the 
result of position measurements performed by the gas particles, because a measurement 
shifts the initial position towards the measured position. The second one is due the 
momentum shift of the momentum measurements, which over the coarse graining time 
6 also results in a position shift. Both these terms are multiplied by a because the 
measurements of light gas particles are weak ones. The third term is due to the (classical) 
collisional momentum exchange, which again, over a time interval 5 gives rise to a 
position shift. 



We have seen in section 12.6 that position diffusion can not be a physical process 
induced by random collisions. Indeed, the second and third terms are clearly due to 
the coarse graining time approximation, while the first term is due to the assumption 
of instantaneous measurements, which in reality take the collision time tc- Therefore, 
all position diffusion terms can be traced back to approximations we used to derive our 
master equation, and are fictitious. 

To show that these terms are small in the range of validity of the approximations 

used, we first discuss which one is the leading one. Here, we use the slow Brownian 

particle limit (noting that we need only qualitative estimates and that for fast Brownian 

particles the result is even more pronounced) to perform the integration 

3 / \ - 

Pg _ ng ( 2kBT\ 2 



/ 



dpgfiriPg 



rUg ) 



(16.29) 



and to use Eq. (14.21) for the collision rate R. Up to a factor Aang^ykBT/mg, we find 
W"^ for the first term, [6h/{2'mW)]'^ for the second term, and d'^ksT / {Qm) for the third 
term. From Eq. ( 15.15[ ) we deduce that the third is large compared to the first, whereas 
according to Eq. (15.17) the third term is also large compared to the first one. 

Therefore, in the range of validity of our master equation, we can neglect the first 
two terms and arrive at 

d(x2> 



dt 



m 



a 



3(1 + a) 



dPgf^riPg) 



p_ 
m 



P9_ 

nin 



(16.30) 
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There is a class of states with {{x,p}) = (e.g. the state \x = 0,p = 0)). If the Brownian 



particle's state is one of these, then the position diffusion term in Eq. (16.30) is the 



leading one. However, there are no such states for which {{x,p}){t) = for all times. 



As we will see in section 16.6, for times longer than 6 the first term will always become 
large compared to the second. Now we remember that our description of QBM is only 
valid on a time scale large compared to 6, and therefore we can indeed neglect the 
position diffusion term. The result is 

dt m 

We want to point out that the limit 5 — )• minimizes the artifically introduced 
position diffusion, although not to zero but to a value connected to the width W of the 
wave packets used to define the approximate phase space projection operator in section 
114.21 Here we also wish to note that the same limit maximizes errors in the decoherence 
rates, because the rate operator then acts as a quantum-limited position-momentum 
measurement on the Brownian particle, and this is not actually performed. To describe 



all aspects of QBM, one should therefore choose some finite 6 according to section 15.3 



16.4 Mean squared momentum 



Applying the adjoint Liouville operator Eq. (16.4) on results in 
1 



-{p'R + Rp') + 



dxg dpg dx dp 



R5{Xg,Pg)B\Xg,Pg]X,p)p^B{Xg,Pg]X,p)j Rs{Xg,Pg)U''{^) (l6.32) 



where [p'^,p^] = and U{t)p^W {t) = p^ were used. Now we apply Eq. (13.18) and (16.7) 



for a calculation analog to Eq. (16.8) to find 
jj dxdpB\xg,pg;x,p)p^B{xg,pg;x,p) 



2 + 2a2 



Tr(<Tp 



s2\ 



2a. 



-Tr 



apvap] + 



(l + a)2 ^ 1 + a 

Using Tr(crp^) = Tr ( \/^p\/^p] = in the limit W ^ oc^ 



[{l-a)p + 2pg]^ 
(l + a)2 



(16.33) 



we substitute into Eq. ( 16.32 ) 



C^p^ 



-Rp^ + 



dXgdpgUi^UXg^Pg) "^^^ ^ ^.l^" '""^ 



^To be precise, one should include these terms as W is finite. Then, we would find an extra momentum 
diffusion term, corresponding to the momentum measurements performed by the gas particles. At the 
end, this contribution would turn out to be negligible compared to the leading diffusion term. Because 



this is all analog to section 16.3 we remove this term already here. 



16.4. MEAN SQUARED MOMENTUM 



155 



dp. ^^^p 



(16.34) 



In the second step we used the definitions ( 14.17) and ( 14.19 ) as weh as U {t)R[pg)p^U\t) 



R{pg)p^ (see Eq. (14.23)). Now we can use Eq. (14.18) to arrive at the final result 



where 
/it(p) 



4mn 



8n„ 



dt 



dpg exp 



-{hriP)) 



-pI 



2mgkBT 



m m„ 



^_ P9_ 
m m„ 



(16.35) 



iP + Pg) 



exp 



[l + a)'^mgy^2TrrngkBT\ \2mkBT 



-ap 



[2a{l-a)mgkBTp^ - 2{mgkBTf 



+ 



ap 



dpg exp 



-pI 



2mgkBT 



[aV + (1 - 2a)plp] 



(16.36) 



For low velocities |f | ^ ^Jkj^Tjrrig of the Brownian particle we can expand hriP) to 
find 



d(p 



-2\ 



8n„ 



dt 



{l+aymg^y2-KmgkBT 



2{mgkBT) — a{2 — a)mgkBT(p' 



s2\ 



a 
12 



(4+a)(p^> 
(16.37) 



Note that this equation gives the correct value for in the steady state limit only 
if a ^ 1. The reason is, that, if the mass of the Brownian particle is comparable to 
the mass of the gas particles, then the Brownian particle can not be considered slow in 
the steady state, and therefore the approximation leading to the above equation is not 
satisfied anymore. 

We therefore confirm the steady state p oc exp^ ^m ^^r ^ foi^ general mass ratios by 
substituting the assumed steady state in Eq. (16.35), and show that d(p^) /dt = 0. 
Disregarding all constants, we have to calculate 



Tr[p/i2 



oc 



oc 



dp < exp 



-{l + a)p^ 
2mkBT 



[2a(l — a)'mgkBTp'^ — 2[mgkBTy 



+ 



ap 



dpg exp 



2 2 ' 

-Pg - ap 
2m„kBT 



[aV + (1 - 2a)p^p] 



(16.38) 
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While the first term is easily integrated to give 

'ia{mgkBT)'^\/27rmkBT 
(1 + a)3/2 ' 



(16.39) 



the second term is not that straightforward as the integration over pg can not be carried 
out analytically. Therefore we will change the order of integration. The integration area 



splits up in to parts as seen in figure 16.1 



Pg = C^P 




Figure 16.1: Integration area of the second term in Eq. (16.38) 



For the first integration area, we write 



dp / <iPgfiP,Pg) 
Jo 



poo poc 

/ dpg / dpfip,Pg) (16.40) 

JO Jvn/a 



'0 Jpg/c 

which is valid for arbitrary functions f{p,Pg)- For the second integration area we use 

/O /-ap f-0 f-0 

dp / dpgf{p,Pg) = - dp dpgf{p,Pg) 

-oo Jo J —oo Jap 



rPg/a 

dPg / dpf{p,Pg) 



oo J —oo 



/O ;"00 
dpg / dpf{p,Pg) (16.41) 

-OO Jn„/a 



~oo Jpg/a 

where the last equality is valid for functions f{p,Pg) which are odd in p. Putting the 
two relations together, 

//•ap p poo 

dp / dpgf{p,Pg) = dpg dpf{p,Pg) (16.42) 

Jo J Jpg/a 
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can be applied to the second term in Eq. ( 16.38 ). The p-integration is then easily carried 
out by substitution of p'^ = u and integration by parts, and we find 



2 J dpg exp 1^ 



/-(l + «)Pg 



2am„kBT 



[{l-a)plMkBT + a^{MkBT) 



21 4:a{mgkBT)'^\/27rmkBT 



(1 + a)3/2 



(16.43) 



where the p^-integration is now over the entire real axis. This cancels exactly with 



Eq. (16.39) and we find 



d{p') 



dt 







(16.44) 



as required for a thermal state. Note that the derivation of this result is valid for all a 
and therefore also holds if the Brownian particle is smaller than the gas particles. 

16.5 Mean position times momentum 



Applying the Liouville operator Eq. (16.4) on {x,p} results in 

£'^{x,p} = [j5^,{£,p}] - ^ (^{x,p}R + R{x,p}^ + jjjj dxgdpgdxdp 

X U{|)^J MXg,Pg)BHXg,Pg;X,p)U^ {|){x,p}U{|)B{Xg,Pg;X,p)^J Rs{Xg,P 

(16.45) 

Similar to the sections before we find 

JJ dxdpB''{xg,pg;x,p)U''{^){x,p}U{^)B{xg,pg;x,p) 

^ (1 _(^Q,-)2 |(1 ~ a)x + 2axg + ^ [(1 - Oi)p + 2pg] , (1 - a)p + 2pg^ 

(16.46) 



where we used Tr\^y/&py/&xj = and similar relations already used in section 
well as the limit W — )• oo. Together with [p^, = —Aihp^ we get 



16.3 



as 



2{f) 



m 



R{x,p} + 



(l + a)2 



dXgdpgU{^)R5{Xg,Pg) 



S^il-a)x + 2aXg + ^[{l-a)p + 2pg],{l-a)p + 2pg^U^{^). 

(16.47) 
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We continue with the help of 
R{x,p} 



2 / dpgfiriPg) II ^^^^\x,pXx,p\ 



P 



Pa 



m m„ 



xp 



(16.48) 



and similar relations already used in section 16.3, as well as Eq. (10.28) to get the final 
result 



d({^,j3}) 



2 if) 



m 



{{xjTip)}) 



(16.49) 



with fxip) from Eq. (16.21). In the slow Brownian particle limit, we can again expand 
/t(p) and find 

d{{x,p}) 2(p2) Ang^2mgkBT 



dt m ^{vfi + mg) 

16.6 Solutions for expectation values 



-({x,j3}). 



(16.50) 



In the limit of a slow Brownian particle (|f | ^ ^Jk^T jm^ , the equations of motion 
derived in the previous subsections are exactly solvable. We will also use the limit of a 
heavy Brownian particle 1+a ~ 1, to ensure the particle stays slow during the evolution. 
Then the equations of motion for the first two moments of position and momentum are 



d{x) 
dt 

d{p) 
dt 

d(x^) 
dt 

d(f)_ 
dt 

<i{{x,p}) 
dt 



{A, 

m 

-lip) 

{{x,p}) ^ a^ng 

2j[mkBT - {p^)] 

-^-l{{x,p}), 
m 



2/cbT 



rua 



3/2 



where we defined 



Aug ^y2rngkBT 



irm 



(16.51) 
(16.52) 

(16.53) 

(16.54) 
(16.55) 

(16.56) 



In the third equation we used Eq. (16.29) to evaluate g{p) for slow Brownian velocities. 
We keep the position diffusion contribution, to show that it is negligible on time scales 
large compared to 6. These are all linear first order differential equations which are 
easily solved using standard techniques. The second equation is immediately solved by 



{p{t)) = exp(-7t)(p(0)). 



(16.57) 
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This in turn can be substituted into Eq. (16.51) and integration results in 



im) 



1 



(x(0)) + — (1 

7717 



(16.58) 



These are the classical motions for position and momentum of a linearly damped particle. 
We start the second moments with 



m 



A;BT(l-e-2^*)+e-2T'*(p2(o)) 



(16.59) 



which we substitute into Eq. (|16.55|) to find 
2kBT 



{{x,p}{t)) 



7 



(1 - e-^' f + ^^^^ (e-^* - + {{x,pm)e-^\ 



m7 



(16.60) 



which in turn can be substituted in Eq. (16.53) to yield 
{x\t)) . 



{x\0)) + 
ksT 



{{x,p}{0)) 



1 - e'^') + 



rwy 

3 - 4e-^* + + t 



(m7)2 



-7t 



y 



2kBT ^ ^^a^ng f 2kBT\ 



3/2 



m7 



3^ \ rUg J 



(16.61) 



Now we only need to substitute 7 and to use Eq. ( 15.16 ) to see that the term involving 



6 can be neglected compared to the other term in the square brackets. Therefore, the 
unphysical position diffusion is indeed small compared to the real physical position 
diffusion, which occurs on a time scale large compared to 7^^. We also have to check 
how the unphysical position diffusion term behaves on short time scales. Expanding the 



exponentials in Eq. (16.61), and substituting 7, we find for the second line 



no 



Ira, 



]{2kBTfl^ 



(4*3 + tS"^ 



(16.62) 



Hence, for times large compared to the coarse graining time, the effect of the unphysical 
position diffusion is indeed small compared to the real position spreading, as promised 



in section (16.3). Because all our results are only valid on a coarse grained time scale. 



we can rightly claim that position diffusion is an artifact of the approximations used to 
derive the Markovian master equation. 

This is very reassuring. After all, this mysterious position diffusion could only be 



explained by random position jumps, and we already showed in section 12.6 that these 
are not possible in collisional quantum Brownian motion. 
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It is interesting to have a look at the long time limit "ft ^ 1 of these equations 



(x(oo)) = (x(0)) + ^^^ (16.64) 

m7 



(p(oo)) = (16.63) 
{x 

(p'^ioo)) = rnksT (16.65) 

({£,p}(oo)) = (16.66) 

(.V)> = (j(0)> + ^i^^ + f^^-^ + ?^. (16.67) 

We find that the mean energy approaches ksT /2 as required by the equipartition theo- 
rem of thermodynamics. Furthermore, the mean square displacement increases linearly 
in time, as is common in diffusive processes. 

Also interesting is the short time limit 7t ^ 1. While the Brownian particle mean 
squared momentum increases linear in time, the mean squared displacement increases 
according to 2kBTjt^ /{3m). This is typical for the process of momentum diffusion. 

To summarize, we find the expected momentum diffusion. On time scales large 
compared to the damping time 7~^, this results in indirect position diffusion, because 



according to Eq. (16.64), each momentum kick q can be associated with a position 
change of q/imj) in the long term. Direct position diffusion (which is an artifact of our 
approximations) is found to be negligible on time scales large compared to the coarse 
graining time scale h/{kBT). 

16.7 Standard form of QBM master equation 

A common master equation for QBM is of the form 

^ = -{[H,p{t)]-^^[x,{p,p}]-^[x,[x,p]]-^[p,[p,p^ (16.68) 

which we call the standard QBM master equation. To derive such an equation for 
collisional QBM |5H [6OI [59] , the limit of a small collisional momentum change of the 
Brownian particle has to be used. But also for non-collisional models, such as the 
Caldeira Leggett model, this form was put forward [5]. For this equation to be of 
Lindblad form, 

D^,,Dpp ^ {hj/Af (16.69) 

is required. By comparing our equations of motions for the lowest order moments to 
the ones derived from the standard QBM master equation, we can find the parameters 
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7, Dxx, and Dpp corresponding to our theory. We will see that this procedure is only 
possible in the heavy Brownian particle limit, and even then it is not certain whether the 
resulting master equation correctly describes other properties of QBM, like for example 
decoherence. However, for completeness and because of the simple form of the standard 
QBM master equation, we proceed with the task despite these problems. 

The equations of motion for the first and second moments of position and momentum 
found from the standard QBM master equation are 

d(x) (p) 



dt 
d{p) 
dt 
d{x^) 

dt 
d(f)_ 
dt 

d({x,p}) 
dt 



m 



-lip) 
{{x,p}) 



+ 2D^ 



m 



-27 (p^) + 2D^ 



pp 



2(p' 



m 



(16.70) 



These can be compared with the equations of motions, which we derived in sections 16.1 
to 16.5, to find values for the parameters 7, Dpp, and Dxx- Interestingly, we can match 
these parameters only for a heavy Brownian particle (a — )• 0) with small velocities 
<C ^/ksT /rUg. That means, only in such circumstances the standard form of the 
QBM master equation might describe the physics correctly. We find 

A.ng-^2mgkBT 



7 



D. 



pp 



D, 



rnksTj 
0. 



(16.71) 
(16.72) 
(16.73) 



Surely these parameters violate inequality (16.69) and therefore do not result in a com- 



pletely positive master equation. However, had we used Eq. (16.53) for the position 
variance and not neglected the term proportional to (5, we had found 

'2kBT' 



D, 



60F 



rrir, 



6\ 



(16.74) 



Using Eq. (15.18) it is straightforward to show that this parameter fulfills the inequal- 



ity (16.69). Therefore we find a completely positive master equation only if we keep 
contributions from unphysical position jumps. So what's going on here? 

The problem is that QBM is not a strictly Markovian open system. Therefore, one 
needs approximations to produce a Markovian equation of motion. Here, as in most other 
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approaches, that is the introduction of a coarse grained time scale. This approximation 
can result in a non-completly positive master equation even if the approximation is very 
well satisfied in the physical system as in QBM in the high temperature and low density 
limit. 

Another prominent example of this sort of problems is the quantum optical master 
equation [5]. The Born-Markov approximation is usually well satisfied, but it results in 
a master equation which is not completely positive. A second approximation is done, 
called the rotating wave approximation, to cast the master equation into Lindblad form. 

In our approach to QBM the introduction of a coarse grained time scale does actually 



result in a completely positive master equation as is seen in Eq. (15.12). However, it 
introduces non-physical position jumps, and if these are removed by hand, it can not be 
expected that the master equation stays completely positive. 



One should keep in mind that here Eq. (16.68) is postulated, and although its pa- 
rameters are chosen to produce the correct equations of motion for the first two moments 
of position and momentum, it is not necessarily the correct master equation describing 
QBM. For a correct treatment of QBM one has to refer to the more complicated master 



equation (15.12). 



Appendix A 

Here we consider operators first shifted in phase space, and then integrated over phase 
space. For this purpose we will need 

f - fff dx' dp' 

Tr(0) = / dx{x\0\x) = /// dx—^-^—{x\x',py^x',p\0\x) 

dx'dp' , ,| A, y I , ,\ ffdx'dp' , ,| A| , ,\ nr^K\ 
dx——— (x ,p \(J \xlx\ X ,p ) = / / (x ,p \(J \x ,p ). (16.75) 

ZTrh J J 2Trh 

Then the relation 

JJ ^^^(^'P)OD\x,p) = Tr(d)i (16.76) 
is shown by its expectation value on a general Gaussian state \x',p') 

{x,p'\ JJ '^^^b{x,p)db\x,p)\x' ,p') = J J '^^^{x -x,p' -p\d\x -x,p' -p) 

= Tr(d), 



which uniquely determines the operator. Applying Eq. (16.76) to Ox 



Tr(Ox)i = JJ ^^^b{x,p)dD\x,p)D{x,p)xD\x,p) 
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= // -^:^D{x,p)dD''{x,p){x - x), 
leads us to a related formula 

JJ ^^^^^(^^p)OD\x,p) = Tr(d)x-Tr(dx)i. 
In a similar manner, the following relations are found 
JJ '^^pD{x,p)dD\x,p) = Tr(0)p-Tr(d]5)i 

[[ ^^^x^D(x,p)dD\x,p) = Tr(0)x2 -2Tr(Ox)x + Tr(Ox2)i 
JJ 27r^ 

JJ ^^^^p'^D{x,p)dD\x,p) = Tr(6)p2 - 2Tr(dp)p + Tr(Op2)i 

JJ '^^^2xpb{x,p)db\x,p) = Ti{d){x,p} - 2Tr(dx)p - 2Ti{dp)x + Tr(d{£,p})i 

(16.81) 



(16.77) 

(16.78) 
(16.79) 
(16.80) 



Appendix B 

Here we evaluate several integrals over Xg, restricted to the region Ss{xg,pg) defined in 



Eq. (14.9). 
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(16.82) 
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dx dp 



Pg_ 
m„ 



m 



(16.83) 



[dxg(xg + p^] [[ dxdpf{x,p) 

ma m 



dxdp f(x^p) < 
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spg_ 0-Xg Xg lOr m < U 
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= 5 dxdp 



P9__P_ 

mo m 



x + 



2m J 12 \mg m 



fix,p) (16.84) 



Chapter 17 



Decoherence 



Coherences of a quantum state are often described by the off-diagonals of the density 
matrix. This leads to the question: What basis should we use to examine decoherence? 
The two bases which first come to mind are the position basis and the momentum 
basis, and indeed, these are the bases usually used in the literature [21 HJl HU [581 m] - 
Nevertheless, they are not without problems. 

If one uses e.g. the momentum basis, one would be interested in how long a super- 
position of the form {\pa)+ \Pb))/V^ survives, or, equivalently, how fast the coherences 
IPaXPbl /2 decrease. There is a major problem in this sort of question: A momentum 
eigenstate (or any state which is very localized in momentum) is itself a coherent super- 
position of widely separated position eigenstates, and there is no guarantee that these 
position coherences of each individual momentum eigenstate do not vanish on the same 
time scale, or even faster, as the momentum coherences of interest! In other words, 
by the time the coherence between \pa) and \ph) decreases, the Brownian particle might 
significantly change its momentum distribution. The same reasoning applies to using 
the position basis. 

Our reservations are clearly related to the concept of a pointer basis, which we out- 



lined in section 2.2 Could we possibly single out a pointer basis from the measurement 



interpretation of a single collision? In chapter 13 we saw that a colliding gas particle 
\xg,Pg)wg performs a smeared out measurement on the Brownian particle in the basis 
\x,p)wi indicating that these states could be used as pointer basis. However, we also 



know from section 14.1 that we have a choice in the width Wg of the gas particle states. 



Does this mean, that we can use Gaussian states of any width W = y^m/mgWg as 
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pointer basis, as long as Wg satisfies the relation (15.21) 



-r^fr^ « Wg « ^, (17.1) 
yJmgkBT rig 



which we required for our treatment of a single collision? 



As we mentioned in section 2.2, there are still many open questions about the emer- 
gence of a pointer basis from the coupling to an environment. Nevertheless, because of 
the heuristic reasoning in the previous paragraph, and the lack of sensible alternatives, 
we will indeed use the Gaussian states \x,p) for the study of decoherence. In particu- 
lar, we will discuss the vanishing of coherences of states of the form \xa,Pa)+ \xb,Pb), 
commonly referred to as cat states in the literature. 

Rather than using a density operator representation in terms of Gaussian states. 



we found it more graphic to use the Wigner function (see subsection 10.2.1) to dis- 
play coherences of Gaussian states. Decoherence can then be discussed in terms of the 
vanishing of the oscillatory behavior of the Wigner function. 

We wish to mention that a Wigner function description of QMB was already put for- 
ward previously in a very heuristic derivation [63] , as well as in a more precise approach, 
but limited to states of the Brownian particle which are close to a thermal state |59j . 
Both articles are concerned about the general form of a partial differential equation for 
the Wigner function, and do not study decoherence. 

In the following, we will compare the Wigner function without a collision, to the 
Wigner function with a collision. This way, we will obtain the decoherence per collision, 



which can be multiplied by the collision rate Eq. (14.21), R ^ ng\j2kBT / ^irmg, to 
obtain the decoherence rate. 



17.1 Collisional decoherence for general gas particles 

In this chapter, we will frequently encounter density operators of the form 

P = \Xa,Pa){Xa,Pa\ + \Xb, Pb^Xb, Pb\ + Cc"^ \Xb, Pb^Xa, Pa\ + Ce~"^ \Xa, Pa){xb, Pb\ ■ 

(17.2) 

Here, c is bound between zero and one and is a measure for the strength of the coherences 
between \xa,Pa) and \xb,Pb), and (p determines the relative phase between these states. 
In particular, an incoming Brownian particle state \xa,Pa)+ \xb,Pb) before a collision 
corresponds to c = 1 and (p = 0. 
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Let us start by deriving a general analytic expression for the Wigner function of 



the state Eq. (17.2). Because the Wigner function is linear in the density operator, it 



can be calculate term by term, using the position representation of the Gaussian states 



Eq. (12.1), as well as the definition of the Wigner function Eq. (10.13). The result is 



Wp{x',p') 



+ 



+ 



1 

irh 
1 

irh 
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XAVD 
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W\p'-p,f 

w^{p'-pa) 



2h 



PAXD , PA 
h XD 



p' 



PD- 



XA 



where we have defined the average position xa 
{xa — Xb) (analogous for momenta). 



h h 
Xa + Xb)/2 and the difference xd 



17.3) 



We noted in subsection |10. 2. ij that the strength of coherences, indicated by oscilla- 
tory behavior of the Wigner function, do not change due to the unitary free evolution. 
Therefore, we will mostly use the interaction picture for what follow j^^ We can then 



use the general formula Eq. (17.3), with c = 1 and (/? = 0, to plot the Wigner function 
of a state 

According to section 



\xa-,Pa)+ without a collision in figure 17.1 (a). 

a collision with a gas particle state \xg,pg), results in 



12.7 



\Xa/b,Pa/b/ 



\Xalb^Pa/h 

c = exp 



a 



ip = 



(l + a)2 VVF2 ' 
2a{xAPD - xdPa) + (1 - a)pgXD - a(l - a)xgPD 



(17.4) 



where we used Xa = x{xg,Xa), Pa 



given by Eq. (12.73) - (12.76). The Wigner 



function after the collision is plotted in figure 17.1 (b). 

It is quite astonishing, that, despite choosing a gas particle with only four percent 
of the mass of the Brownian particle, and, despite using a superposition of very close 
Gaussian wave functions, almost all coherences are lost after a single collision. If we had 
separated the initial Gaussians only slightly more, or had chosen only a slightly heavier 



gas particle, the coherences would be not visible at all, because c in Eq. ( 17.4) decreases 



^The time evolved Wigner function is obtained by replacing x' by [x' — p't/m) in Eq. ( 17.3 1. Because 



this would unnecessarily complicate matters, we will use the interaction picture Wigner function, i.e. 

without time evolution. 

^We could also say, we study the coUisional process according to figure 12.1 (b). 
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(a) (b) 




Figure 17.1: The Wigner function of the initial coherent superposition \ip)= |a;a,Pa)+|2;6,P6)(a), 
and of the state resulting from a collision with \xg,pg) (b). Parameters are: Xa = 15, X}, = 
0, Pa = 0, = 1.5, W = A, m = 1, = 1, Xg ^ 100, Pg = -1, and a = 0.04. 



exponentially with these parameters. This observation is independent on the initial 
momentum and position of the colliding gas particle, as well as whether the Gaussian 
wave functions are separated predominantly in position or momentum. It is therefore 
fair to say, that, unless the gas particle is much lighter than the Brownian 
particle, the decoherence rate equals the collision rate. This result becomes 
even much more pronounced, if we average over different initial gas particle positions 
and momenta, and we will study this effect in the following section. 



17.2 Collisional decoherence for very light gas particles 



Figure 17.2 shows how the Wigner function of a superposition state changes due to a 
collision with a very light gas particle. The coherences after the collision are still well 
pronounced. The reason is that the measurement, performed by a very light gas particle 
on the Brownian particle, is so imprecise that it can not distinguish between the two 
Gaussian wave functions of the Brownian particle. 

The main effect of a collision with a very light gas particle is a shift of the entire 
Wigner distribution in phase space. That the oscillating part exactly shifts with the 
Gaussians, without obtaining an additional relative phase, might seem surprising, as 
this would certainly not be the case if we simply changed to Ix^/f,, Pa/b) (i-e. 
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the relative height of the oscihations would be different). The reason can be found by 



looking at the argument of the cosine in Eq. (17.3) (and therefore applies also to heavy 
gas particle, where the effect is less easily observed, because of the strong suppression 
of coherences). In particular, the last two terms in the argument account for a phase 
shift according to the change of the average values pA — ^ Pa and xa — ^ xa- The first 
two terms in the argument account for an additional phase shift, but their sum does 
not change in a collision, because {xaPd — PAXd) / {'^K) equals {xaPd — Paxd) / {"^K) + 



with if taken from Eq. (17.4). 

Of course, if the gas particle state \xg^pg) is taken from a thermal gas, we have to 
average over all gas particle momenta pg, weighted by the Maxwell-Boltzmann distribu- 
tion fiTiPg), as well as over all gas particle positions Xg, which can reach the Brownian 
particle in a given time interval (0, t). Because each possible combination of Xg and pg 
results in a different shift of the Wigner function in phase space, it is clear that this 
procedure strongly removes the oscillations. This is the decoherence effect we referred 



to 'phase averaging' in section 12.7 It can still be very strong, even if the measurements 
which the gas particles perform are very weak. 

In the limit of a small mass ratio a, the Brownian particle will be very localized 
compared to the gas particles, in both, position and velocity. Therefore, we do not have 



to use the full rate operator from section 14, but can simply assume that a gas particle 
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collides with the Brownian one during a time interval (0,t), if < —XgiUg/pg < t. 
17.2.1 Position decoherence 

To study position decoherence, we consider an initial Brownian particle state \xa,p) + 



\xh,p). An example of the corresponding Wigner function in plotted in figure 17.3| (a). 



How the Wigner function after a collision with a gas particle state \xg,pg) is obtained, is 



described in section 17.1, but now we also have to average over Xg and pg. To this end, 
we will randomly choose pairs {xg,pg) according to an appropriate distribution, and for 
each we will calculate the Wigner function. Finally, we will average over all of these to 
obtain the average Wigner function Wp^ after one collision. 

First, we need the normalized collision probability distribution for pg 

RiPg) 



C{Pg 



I dpg R{Pg 



exp --^f— , (17.5) 



2mgkBT \ 2mgkBT 



which we obtain from Eq. (14.18) (or Eq. (9.8)) by assuming a slow Brownian particle. 



We will consider the cases ^ separately, starting with pg < 0. We therefore choose 
a random number ui G (0, 1/2), and solve 

fPa ^ ^ 

Ul = dPgC{Pg) 

J —OO 



^ Pg = -2mgkBT\n{2ui) (17.6) 

to obtain a random gas particle momentum < 0. Because of symmetry of C{pg), 
we use Pg = —2mgkBT ln{2ui) for pg > (after choosing a new random number ui). 
For a given gas particle momentum pg < 0, the position Xg has to be in the interval 
Xg G {0,pgt/m). Therefore, we chose a second random number U2 G (0, 1), and use 

Xg = -U2^ (17.7) 
m 

as random gas particle position Xg. 

We note at this point, that the position distribution of the colliding gas particle does 
not only depend on the gas temperature, but also on the time interval (0,t). This is to 
be expected as the initial position of the gas particle is linearly related to the time of 



collision (see figure 12.1), and the time of collision can be anywhere within the interval 
(0, t) of consideration. 
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Figure 17.3: (a): The Wigner function for the Brownian particle before a colhsion. (c), (e), (g): 
The average Wigner function after one collision with a gas particle at temperature T = 0.2, 0.5, 
and 1.5, respectively, (d), (f), (h): The change of the Wigner function due to a collision at 
temperature T = 0.2, 0.5, and 1.5, respectively, (b): The relative change of the Wigner function 
at the origin due to a collision. This serves as a quantitative measure of the decoherence per 
collision. The solid lines are the first order expansion in temperature. Parameters 
20, Xb = -20, Pa ^ 0, pb = 0, W = 4, m = 1, h ^ I, ke ^ I, t = 20, and a = 0.0001. 
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The Wigner functions after a collision, averaged over 200 pairs {xg,pg), are shown 
in figure 17.3| (c), (e), and (g), for the respective temperatures T = 0.2, 0.5, and 1.5. 
The corresponding changes of the Wigner function from the initial one (figure 17. 3| (a)) 



are shown figure 17.3 (d), (f), and (g). We use the maximum peak of the oscillations, 
relative to the maximum peak of the initial state's oscillations, as a quantitative measure 
of coherence. In particular, the relative change of the maximum peak serves as 'deco- 



herence per collision', and is plotted over the temperature T in figure 17.3 (b). There is 
obviously a linear temperature dependence of the 'decoherence per collision', as long as 
the momentum transfer in a collision is small. The small constant decoherence on top 
of the one linear in temperature is due to the very imprecise measurements, performed 
by the gas particles. 

To explain this behavior, let us have a closer look at the oscillating term of the 



Wigner function Eq. (17.3) for the initial state \xa,p)+ \xb,p), 



2c 



exp 



exp 



[X -xa) 



cos 



PXD , P-P 



:i7.; 



with c = 1 and = 0. According to section [l7.1 a collision with a very light (1 + a ~ 1) 
gas particle results in a change of c = ex.pi^—ax'j^/W'^^, which is the responsible for a 
small decoherence independent of the temperature. Next, the position and momentum 
change to Xafi = a^a.fe + 2axg and pa,b = Pa,b + 2pg, respectively. This slightly broadens 



the exponentials in (17.8) when averages over Xg and pg are taken, but the effect is 



too small to be visible in the parameter range of figure 17.3 Further, we showed in 
section 17.1, that [(p — px d / {2h)] does not change due to a collision. Therefore, the only 
significant effect is the change of p to p + 2pg in the cosine, which is a random phase 
shift corresponding to the random momentum of the colliding gas particle. 

As we are interested in the reduction of the peak of the oscillations, we choose p' 
such that cos[ip — px^/i^H) + xd{p — p') /K\ = 1, i.e. tp — pxo/i^h) + xd{p — p')/h = 0. 
The same term after a collision at the very same p' is then 

'2xDPg'' 



cos ^ip + 



-pxD , P + 2pg - p' 

+ XD- 



2h 



cos 
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XDPgV 

h ) 



(17.9) 
(17.10) 



where we expanded the cosine around its maximum. For the 'decoherence per collision', 
we have to average over the squared momenta of the colliding gas particle. Using 



Eq. (17.5), we find 



(^"3)0011 



2m„kBT. 



(17.11) 
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Note that the average squared momentum of a coUiding gas particle is not the average 
squared momentum of a general gas particle, (p^) = rUgkBT, obtained from the Maxwell- 
Boltzmann distribution (faster gas particles are more likely to collide) . We finally deduce 
for the 'decoherence per collision' for position superposition states 



Decoherence 



AmghsT 2 



(17.12) 



Collision It? 
Note that we would also have to add axjj/W'^ to account for the decrease of c, but 



according to relation (15.17), this term is small compared to Eq. (17.12) in the range 



of validity of our work. This establishes, that the leading form of decoherence is phase 
averaging, and not information exchange between Brownian and colliding particles]^ 



The solid lines in figure 17.3 (b) are calculated from this formula, and agree very 
well with the numerical results as long as the temperature is low enough to justify the 



expansion Eq. (17.10). Multiplying the last equation with the collision rate, we find the 



position decoherence rate 



2^/i2 



-X 



D- 



(17.13) 



This very nice result holds when (2x£))^mgkBT <^ fi? is valid, and agrees up to some 
constants with the decoherence rate obtained by Joos and Zeh [Z. Phys. B: Condens. 
Matt. 59, 223 (1985)] 



We can generalize formula Eq. (17.12) to high temperatures and large position sep- 



arations, by directly averaging the cosine of Eq. (|17.9|), where we again use Eq. (|17.5|) 
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(17.14) 
(17.15) 



In the second step, we used integration by parts. We deduce for the 'decoherence per 
collision' of position superposition states 



Decoherence 2x d \/2mgkBT 



Collision 



h 



du e 



sm 



/ 2xDj2mgkBT 

\ 



h 



(17.16) 



The non-zero y-intercept in the inset of figure 17.3 (b) is due to c. But the reduction of coherences 



from this information exchange is only dominating at temperatures so small, that relation ( |15.17[ ) is 
violated and our approach is not valid anymore. Therefore, we can not eliminate the possibility that 
information exchange between the Brownian particle and the gas might not contribute to decoherence 
at all! 



174 



CHAPTER 17. DECOHERENCE 



(a) 



Decoherence 
Collision 



(b) 





Figure 17.4: (a): The 'decoherence per coUision' plotted over the temperature, for parameters 
as in figure [T773| but for larger temperatures. The dots arc from the numerical Wigner function 
at the origin, and the line is from Eq. ( 17.15 ). The small discrepancy is due to the change of the 



Gaussians in Eq. (17.8) 



(b): The Wigner function after a collision at temperature T — 7 shows indeed oscillations in 
opposite phase to the initial Wigner function (see figure 17.3 (a)). This leads to a decoherence 
rate which exceeds the collision rate. 



or for the decoherence rate 

AxDUgkBT 



due sm ^ , u 



(17.17) 



The 'decoherence per collision' is plotted over temperature in figure 17.4 (a) for the 



same parameters as in figure 17.3 (b), but for higher temperatures. It might come as a 
surprise, that the decoherence rate exceeds the collision rate for XD^/^nigkBT /h > 1. 
The reason is that in this regime, the Wigner function after a collision shows oscillations, 
which are out of phase with the oscillations of the initial Wigner function, as shown in 



figure 17. 4| (b). Therefore, if we write down the actual (interaction picture) Wigner 
function after some time t as 



W, 



Pit) 



1 - Rt)Wp,, + RtWp,, 



(17.18) 



where R is the total collision rate Eq. (14.21), and po and pi are the density operators 



corresponding to either no or one collision, respectively, then the oscillations in Wp^ and 
Wp^ interfere destructively. 

We note that it is often stated in the literature that, if the separation X£) of two 
interfering wave packets is larger than the thermal wave length A = h/ y^27rmgkBT of 
the gas, then a colliding gas particle can distinguish between the two interfering wave 
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packets, therefore removing their coherences. In finding that the decoherence rate is 
about the cohision rate if x/j > we confirm the latter part of this statement, but 

we also show that the loss of coherence is by no means due to the information exchange 
between the colliding particles, but due to phase averaging. This in turn results from a 
relative phase, which depends on the momentum of the colliding gas particle. 



17.2.2 Momentum decoherence 



Now, we consider the initial Brownian particle state \x,pa)+ \x,pb), whose Wigner func- 



tion is plotted in figure 17.5 (a). Again, the main source of decoherence will be phase 
averaging. But contrary to the previous subsection, where the relative phase of the two 
Gaussians wave packets after a collision depended on the initial momentum of the collid- 
ing gas particle, for momentum decoherence this phase depends on the initial position 
of the colliding gas particle. The more variation of the initial position of a colliding 
gas particle, the more 'decoherence per collision' we will find. For a given gas particle 
momenta Pg, the gas particle position can be anywhere in the interval {—pgi/nigjO). 
Therefore, the decoherence per collision will not only depend on the temperature T (for 



the distribution of pg), but also on the considered time interval. In figure 17.5 (c) - (h), 
the eff'ects of a collision on the Wigner function is shown for diff^erent temperatures and 
time intervals. The dependence of the 'decoherence per collision' for short times and/or 
low temperatures turns out to be linear in temperature and quadratic in time, as shown 



in figure 17.5 (b). As a result, it is not possible to define a time independent decoherence 



rate for momentum superpositions. We will provide a physical interpretation at the end 
of this subsection, and first give a mathematical explanation of these results. 

For this purpose, we consider again the oscillating term of the Wigner function 



Eq. (17.3) 
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(17.19) 



at its maximum position (f + xpo/C^K) —pd{x — x')/h = 0. As discussed in section 17.1 



a collision does not change the sum ip + xpo/i'^K) in the cosine, and we only have to 
consider the change of x to x + 2axg in the last term in the cosine. Therefore, we find 
the 'decoherence per collision' by averaging over cos{2ax gpo/^)- Here, we need the 
normalized distribution C{xg) of initial positions of colliding gas particle. The phase 



space distribution of gas particles is given by pg{xg,pg) in Eq. (9.6). For a given Xg < 0, 
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a gas particle collides if its velocity is larger than \xg\/t, which leads us to 



C{Xg) oc 



n„ 



y'2'nmgkBT J\xg\mg/t 



dpg exp 



2mgkBT 



(17.20) 



The distribution is normalized either by integration over Xg, or directly by dividing by 



the total collision probability Rt. After substituting u = pg/ y'2mgkBT we find 
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(17.22) 



where we changed the order of integration (see figure 16.1 for detail) in the second line. 
As this functions represents the coherences after one collision, we have to subtract it 
from one to get the 'decoherence per collision' 
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(17.23) 



(17.24) 



The approximation is obtained by expanding the sine to third order. The solid lines in 



figure 17.5 (b) are taken from Eq. (17.23), and agree with the data obtained by averaging 
over Wigner functions with random gas particle positions and momenta]^ 



We recognize from Eq. (17.23) and (17.24), that for short time intervals, it seems 



as if collisions can not reduce the coherence of momentum superpositions states. The 
reason is that momentum decoherence is not a leading order process, but a higher order 
one, resulting from position decoherence due to the position separation XD{t) = tp£)/m, 
which a Brownian particle with momentum separation pu acquires over time. 



The non-zero y-intercept in the inset of figure 17.5 (b) is again due to c (in addition to the spreading 



of the Gaussian in Eq. (17.191). Similar to position decoherence, the reduction of coherences from this 



information exchange is only dominating at temperatures and times so small, that relation (15.151 is 
violated, and our approach is not valid anymore. 
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Figure 17.5: (a): The Wigner function for tlic Brownian particle before a collision, (c), (e), (g): 
The average Wigner function after one collision with a gas particle within a time interval (0,i) 
at temperature T, with < = 20, T = 0.5 in (c); i = 50, T = 0.5 in (e); i = 20, T = 3 in (g). (d), 
(f), (h): The corresponding change of the Wigner function due to a collision, (b): The relative 
change of the Wigner function at the origin due to a collision as function of the time interval for 
T — 0.5 as well as T = 1. Parameters are: Xa = 0, xi, = 0, Pa — 1.2, pb = —1.2, W = 4, m = 
1, /i = 1, fcs = 1, and a = 0.0001. 
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Let us investigate this proposition in detail, by splitting the time interval (0, t) in 



many, but small intervals. Then, according to Eq. (17.24), there should be no momen- 
tum decoherence at all, but the two wave packets will have a time dependent position 
separation X£){t') = t'po/m. Therefore, if a collision happens during the time inter- 
val + dt), the coherence after the collision is given by substituting xd = fpo/m 



into the equation for position decoherence, Eq. (17.14). Averaging over all possible 
+ dt) C (0,t), we find for the coherence after a collision during the time interval 

(0,t) 

' dt / du2ue cos ^ ' 



t .Iq Jq \ mh I 

mh _„2 . / 2t^2mgkBTpD 

' due sm ^ 



ty^2mgkBTpD Jo \ mh 



which is exactly the same as obtained by Eq. (17.22). 

Therefore, we arrive at a very neat interpretation of the decoherence process of 
a superposition of two Gaussians wave packets. The decoherence due to a collision 
depends on the position separation of the two Gaussian wave packets at the time of 
the collision There is no direct decoherence due to the momentum separation p/j of the 
two coherent wave packets. Over time, the momentum separation changes the position 
separation according to X£){t) = X£, + tp£,/m. This leads to an indirect influence of p/j 
on the decoherence rate, which, if there is no initial position separation, is described by 
Eq. pT23l). 



17.3 Conclusions from the study of decoherence 

We have seen in this chapter, that we can disregard the decoherences due to the mea- 
surements, which are performed by colliding gas particle. The reason is, that the mea- 
surement effects are small compared to phase averaging for very light gas particles; and 
that phase averaging destroys all coherences for slightly heavier gas particles, such that 
for both situations, the decoherence due to the measurements is not visible. Therefore, 
the entire dynamics of the Brownian particle due to a collision with a thermal gas par- 
ticle seems to be captured, by changing the momentum in the Wigner function at the 

^As the time of collision we understand here the time, when the mean value of the gas particle wave 
packet crosses over the mean value of the Brownian particle wave packet, such that a collision results in 
a change of pa , but not of xa ■ 
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time of the collision according to 

Wp{x',p') ^ Wp{x',p'-q), (17.25) 

where the momentum change is q = 2(pg — ap')/{l — a). This reminds very much of the 
change of a classical phase space probability distribution due to a collision, as described 
in section 9.3 Noting that the Wigner function Wp{x,p) behaves like a classical phase 



space probability distribution also during free evolution, it follows that it obeys the 



same master equation Eq. (9.12) as the phase space probability distribution p{x,p) 
of a classical particle. If the momentum transfer is small compared to the variations 
of the Wigner function in momentum space, i.e. Wp{x',p') ~ Wp{x',p' — q), then the 



Wigner function also obeys the classical Fokker-Planck equation Eq. (9.25), although 
the condition of small momentum transfer is more stringent than in classical dynamics, 
because of the possibly rapid oscillations of the Wigner functions. 

After the huge effort of finding a very complicated master equation for the density 
operator of the Brownian particle, it finally seems possible that quantum Brownian 
motion could be much easier described, by applying the equations which govern the 
dynamics of the phase space probability distribution of a classical particle, to the Wigner 
function of a quantum particle. More work is needed to rigorously establish the described 
connection between classical and quantum Brownian motion. 

One should mention, that such a description of QBM would include all classical 
features like friction and momentum diffusion, but not the unphysical QPD, which all 
master equations of Lindblad form show. On the one hand, that is a very nice feature, 
but on the other hand, that also means that the dynamics described by such means, can 
not be completely positive. Therefore, we once again observe, that Markovian QBM is 
only an approximation to reality, which either leads to non-positive dynamics, or to a 
fictitious quantum contribution to position diffusion, but can not avoid both. 

It is clear, that the process of a single collision with a gas particle might be quite 
different in three dimensions, including the measurement performed by a colliding gas 
particle. But if it turns out, that the decoherence effects due to the measurement is small 
compared to the phase averaging, as it is in one dimension, then also three dimensional 
quantum Brownian motion might be described by applying the classical equations to 
the Wigner function of the Brownian particle. 

The measurement performed by a colliding particle depends on the interaction po- 
tential. But if the measurement performed by a thermal particle is negligible compared 
by the decoherence effects due to the random momentum transfer, then it might be 
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possible that also scattering of other particles, such as photons, could be described by 
such a Wigner function approach. 



Chapter 18 

General Conclusions 



We first applied the quantum measurement theory approach to the study of the coherent 
transfer of an electron along a rail of quantum dots. If used for the transfer of quantum 
information, we found, that the spin of an electron should be preferred as a qubit, 
rather than its position. Furthermore, we were able to specify the transfer fidelity for a 
given strength of Markovian and certain non-Markovian dephasing sources. This led us 
to conclude, that under certain circumstances, it is preferred to divide a long transfer 
distance into several smaller ones. 

We then used quantum measurement theory, to study collisional quantum Brownian 
motion, a topic which has several controversial results. We first found, contrary to 
common belief, that there is no additional quantum contribution to position diffusion. 
Next we showed, that a colliding particle carries away information of the Brownian 
particle position and momentum. Finally, this information transfer turned out not to 
be the main source of decoherence of superposition states. Instead, phase averaging 
due to a random relative phase, which depends on the momentum of the colliding gas 
particle, was found to quickly and efficiently reduce coherences. 

The first statement of the above paragraph seems to be hard to test experimen- 
tally [58], because the classical position diffusion would quickly dominate anyway. The 
second statement can not be observed in a thermal gas because of the third statement, 
although one might be able to measure the momentum of a particle by a collision with 
a smaller particle, if the latter is in a well prepared initial state. The last statement 
however, could be accessible to observation in e.g. double slit experiments, because it 
predicts, that, for certain experimental parameters, the decoherence rate can exceed the 
collision rate. This effect would be impossible, if information transfer between colliding 
particles were the source of decoherence. 
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Because this effect is predicted to occur at higfier temperatures and witfi larger 
spatial separations of the superposition state, it might be more accessible to experi- 
ments, compared to the quadratic dependence in spatial positions separation for lower 
temperatures and lower spatial separation, which was already predicted previously. 

Our results also let us to conjecture, that comparatively simple equations for the 
classical phase space probability distribution could be applied to the Wigner function 
of the Brownian particle, to successfully describe quantum Brownian motion, possibly 
even for much more general situations as considered in this thesis. 

In conclusion, even in situations where the information transfer between system of 
interest and its environment does not play a significant role, the measurement approach 
was still able to lead us to considerable insight. However, instead of studying the effect 
of an environment on the system directly, but rather to first considers the effect of 
the system on the environment, and then to draw conclusions about the information 
transfer between both systems and eventually about the dynamics of the system itself, 
is somewhat indirect. This can result in considerable work if carried out rigorously, as 
especially the second part of this thesis shows. This possibly limits its value to situations, 
where other approaches failed to deliver conclusive results. 

On the other hand, a measurement interpretation of open quantum system might be 
a powerful tool if used heuristically, as in the first part of the thesis, where we did not 
consider the detailed interaction between electron and quantum point contacts. This 
might be true especially, if a first principle approach does not seem feasible. 
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